ZEITSCHRIFT FÜR ANGEWANDTE 
MATHEMATIK UND MECHANIK 


-INGENIEURWISSENSCHAFTLICHE FORSCHUNGSARBEITEN 
Band 34 August/September Heft 8/9 


Jahrestagung der Gesellschaft für angewandte Mathematik 
und Mechanik 


Die Jahrestagung 1954 der Gesellschaft für Angewandte Mathematik und Mechanik fand in 
der Osterwoche vom 20.—24.4. in München statt. München hatte seine alte Anziehungskraft 
bewährt. Der Besuch der Tagung war groß (etwa 270 Teilnehmer). Entsprechend groß war auch 
die Zahl der Vorträge. Es fanden 6 längere, zusammenfassende und rund 70 Kurzvorträge bzw. V or- 

 führungen wissenschaftlicher Filme statt. Diese große Zahl machte es erforderlich, daß drei V ortrags- 
reihen nebeneinander liefen, eine für angewandte Mathematik, eine für Mechanik und eine dritte 
für Strömungslehre. 

Die Hauptversammlung der Gesellschaft fand am 23.4. vormittags statt. Im vergangenen 
Jahre war eine auffallend große Zahl von Mitgliedern verstorben. Ihrer wurde bei Beginn der Sitzung 
gedacht. Dann berichtete der Vorsitzende Prof. Dr. Walther über das verflossene Jahr. Insbeson - 
dere sprach er über die Fachausschüsse, durch deren Arbeit die Haupttagung von. Einzelberichten 
entlastet werden soll. 

1. Den Fachausschuß für statistische Methoden, den bisher Herr Stange geleitet hatte, gab 

dieser wegen Arbeitsüberlastung ab. Ihn wird Herr Heinhold, München, übernehmen. 

2. Den Fachausschuß für moderne Anwendungen der Funktionentheorie wird weiter Herr 
Ullrich, Gießen, betreuen. 

3. Der Fachausschuß für Rheologie hat im vergangenen Jahr eine Tagung abgehalten, über 
die der Leiter, Herr Sc hultz-Grunow, Aachen, auf der Münchner Tagung zusammen- 
fassend. berichtete. 

4. Der Fachausschuß für Rechenmaschinen unter der Leitung von Herrn Walther, Darm- 
stadt, und der Fachausschuß für Rheologie werden ihre nächste Sondertagung zusammen 
mit den zuständigen Fachausschüssen anderer Organisationen im Frühjahr 1955 veranstalten. 

5. Weiter besteht ein Fachausschuß für Regeltechnik unter Leitung von Herrn Oppelt, Frank- 
furt. Dieser gab ebenfalls auf der Münchner Tagung einen Überblick über Probleme und 
Methoden der mathematischen Behandlung von Regelvorgängen. Im Anschluß daran hielt 
der Ausschuß seine erste Sitzung ab. 

6. Neugebildet wurde ein Ausschuß für Strömungsforschung, dessen Leitung Herr Rvegels, 
Göttingen, übernahm. 

Den Kassenbericht und den Bericht über die Mitgliederzahl, die im Berichtsjahr auf 582 
zugenommen hat, gab der Geschäftsführer Prof. Dr. Görtler. Auf Antrag der Kassenprüfer 
wurde ihm Entlastung erteilt. Satzungsgemäß schieden aus dem Vorstandsrat die Herren Betz, 
Göttingen, und Willers, Dresden, aus. Neugewählt wurde für Herrn Betz der diesjährige 
örtliche Tagungsleiter Herr Heüinhold. Der Unterzeichnete wurde wiedergewählt. Das Amt der 
Kassenprüfer übernahmen die Herren Quade, Hannover, und Cremer, Aachen. Der J ahres- 
beitrag wurde wieder auf 5.— DM festgesetzt. Die nächste Tagung wird in der Zeit vom 31. Mai bis 
4. Juni 1955 ( Pfingstwoche) in Berlin stattfinden. Anschließend an die Hauptversammlung wurde 
die Ausbildung und Stellenvermittlung der Diplom-Mathematiker besprochen. Es wurde beschlossen 
einen Ausschuß für die Angelegenheiten der Diplom-Mathematiker ins Leben zu rufen, mit dessen 
Leitung Herr Schultz-Grunow, Aachen, beauftragt wurde. 

Wie auf allen bisherigen Tagungen war auch auf dieser für Geselligkeit und Gelegenheit zur 
Aussprache gesorgt. Am Dienstag fand ein zwangloser Begrüßungsabend im Regina-Palasthotel 
statt, am Mittwoch abend ein Besuch des Deutschen Theaters. Donnerstag mittag hatte der Ober- 
bürgermeister der Stadt München zu einem Empfang im Ratskeller geladen. Abends versammelte 
man sich im Augustiner-Keller, und am Sonnabend fand ein Ausflug in das Gebiet des Ammer- und 
Starnbergersees statt. Außerdem hatten die Damen der Teilnehmer Gelegenheit, am verschiedenen 
Besichtigungen und Ausflügen teilzunehmen. 

Für die Organisation dieser ausgezeichnet verlaufenen Tagung darf ich im Namen aller Teil- 
nehmer dem örtlichen Tagungsleiter Herrn Heinhold herzlichen Dank aussprechen. 


Willers 
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A. Angewandte Mathematik 


Der Konvergenzgrad bei iterativer Lösung von Gleichungen 
durch inverse Interpolation 
Von K.-H. Bachmann in Dresden 

Gegeben sei eine reelle Funktion /(x), die in einem Intervall 3 genau eine einfache Null- 
stelle & besitze. f(x) sei in stetig und genügend oft stetig differenzierbar. Um & zu bestimmen, 
geht man oft in der Weise vor, daß man f(x) durch eine Ersatzfunktion g(x) ersetzt, deren Null- 
stellen einfach zu bestimmen sind. Ein Iterationsverfahren zur Bestimmung von % erhält man 
dadurch, daß man von n + 1 Anfangsnäherungen 2;<%( = 0, 1,..., n) ausgehend eine Funk- 
tion g,(z) aus den Forderungen \ 


a EEE sy) 


bestimmt. Es möge nun g, (x) in { genau eine einfache Nullstelle x,;ı haben, die als nächste 
Näherung für & betrachtet wırd. Allgemein bestimmt man g,(z) aus 


9,(2:4,—1) = (1) 2: an Basar ie Se ce ER SA RM (2) 
und &,+, als Nullstelle von 9,(x) in %, wieder vorausgesetzt, daß g, in \$ nur eine einfache Null- 
stelle hat. Fordert man noch, daß g, in % differenzierbar sei, erhält man durch Anwendung des - 
Mittelwertsatzes für den Fehler von %,+, 


„ _ MW) 9@) 


a »„—% 7 7 a Si ee ee ehe 3). 
Ent» = Int Tamm e 6) 


Der Nenner verschwindet nicht, da g, in % nur eine einfache Nullstelle hat. 


Im folgenden seien die g, stets Polynome n-ten Grades. Dann ist die Bestimmung von g, 
aus (2) eindeutig möglich und es ergibt sich der Fehler 
n 


u) 
(n Eh 1)! g} En ITen-l (Min 41 sn <Max &r—ı) (4). 


Sind die Voraussetzungen 


lent,| = 


fe+OG) 
(n + 119,8 
ORG Ca NEN. eg 


erfüllt, wobei e eine obere Schranke für die |&| @=0,1,..., n) ist, so konvergiert die Folge {x,} 
gegen &. (V1) und (V2) sind erfüllbar, wenn e genügend klein ist, da dann $ so klein gewählt 
werden kann, daß für jedes mögliche Interpolationspolynom g’ keine Nullstelle in hat und außer- 
dem dann jedes Interpolationspolynom in “% genau eine einfache Nullstelle hat. 

Um die Schnelligkeit der Konvergenz zu untersuchen, werde eine zur Folge {|e,]} majorenne 
Folge {ö,} betrachtet, definiert durch 


SO. (EI fr 


n 
On SO A 0 (=V, RE N) - Te ne En 0) 
i—=0 
durch wiederholtes Einsetzen erhält man aus (5) 


On te I 
Die &, und ß, genügen den Rekursions formeln 


m lt Sa; RUN) VENEN TER 2 re 7), 
Pat» >= 2 Bi-i; ßi == 1 (Ü = 0, 1; eo. n) ar 1 Te (8). 


Maßgebend für die Schnelligkeit der Konvergenz ist das Wachstum der ß,. Definiert man nun 
den Konvergenzgrad g einer Folge {z,} — 2 als obere Grenze aller y, für die 


sup a RAT ee (9), 


!)Vgl.z.B.Bodewig, Z.angew. Math. Mech. 29, S. 44-51. 
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erhält man im vorliegenden Falle als Konvergenzgrad 


g=A = re ee 
Es folgt dies daraus, daß 2 
re el 


ist und die Folge der Exponenten von Q beschränkt bleibt. A, ist die absolut größte Wurzel der 
charakteristischen Gleichung der Differenzengleichung (8), falls diese Wurzel reell und einfach ist. 
Die charakteristische Gleichung lautet 


n 
a a ee, u. 
u=0 
Ihre absolut größte Wurzel ist positiv reell und einfach und genügt der Abschätzung 
n+1 
Een 


Für n = 1 ist das beschriebene Iterationsverfahren die bekannte regula falsi. Hierbei ergibt sich 
der Konvergenzgrad zu 


- 5 (5 + 1) = 1,618. 


Unter entsprechenden Voraussetzungen lassen sich diese Betrachtungen auf die Lösung 
von Gleichungen in allgemeineren Räumen ausdehnen. 


Quadratische Konvergenz der Bernoullischen Methode 
Von F.L. Bauer in München 


Das Bernoullische Iterationsverfahren zur Nullstellenbestimmung von Polynomen 
einschließlich der auf Jacobi zurückgehenden Erweiterung für konjugiert-komplexe Paare 
ist nur linear konvergent. Es läuft bekanntlich auf die Iteration eines Vektors mit der zum 
Polynom gehörigen Frobeniusschen Matrix hinaus. Zunächst wird gezeigt, daß die Iteration 
mitder Transponierten(sie hängt mit der Lagrange schen Variante des Verfahrens zusammen) 
den Vorteil einer einfachen Verprobungsmöglichkeit hat und programmierungstechnisch so 
günstig ist, daß sie sich sogar mit gewöhnlichen Lochkartenanlagen durchführen läßt. Ins- 
besondere aber erlaubt sie eine spezielle Maßnahme, durch die einquadratisch konver- 
gentes Fortschreiten der Iteration erzielt wird. Dazu hat man jeweils n Schritte der gewöhn- 
lichen Iteration (wenn n der Grad des Polynoms ist) durchzuführen und über die Zwischen- 
ergebnisse geeignet zu mitteln. Die Gewichtsfaktoren der Mittelung ergeben sich (und zwar eilt 
das nur für die Iteration mit der Transponierten) nebenbei in unmittelbarer Weise. 


Auch das solcherart abgekürzte Verfahren ist gegen Rundungsfehler selbstkorrigierend. 
Außerdem ist der Iterationsprozeß weiterhin programmierungstechnisch einfach und auch mit 
einer Verprobungsmöglichkeit versehen. Das Verfahren ist daher den Erfordernissen einer 
Durchführung mit programmgesteuerten Rechenanlagen angemessen. Da im wesentlichen nur 
die ersten drei Speziesoperationen vorkommen, ist es auch für Polynome mit komplexen Koeffi- 
zienten geeignet, der Rechenaufwand steigt nicht mehr als unerläßlich an. 


Die Iteration führt nicht nur auf die absolut größte Nullstelle, sondern auch direkt auf 
das Polynom, von dem diese Nullstelle abgetrennt ist. Das Vorkommen mehrerer Nullstellen 
gleichen größten Betrags macht sich bekanntlich durch die Nicht-Konvergenz der Hauptitera- 
tionsreihe bemerkbar; man hat nun mit dieser Reihe einen Prozeß durchzuführen, der formal mit 
der ‚„‚eross-division‘‘ von Routh zusammenfällt, und zwar unter Umständen wiederholt, bis 
eine der so gebildeten Nebenreihen konvergiert. Dann wird nicht nur das (für reelle Polynome 
selbst reelle) Polynom dieser Nullstellen, sondern ebenfalls auch das Polynom der verbleibenden 
Nullstellen approximiert. Von besonderer Bedeutung ist hier natürlich der Regelfall eines ein- 
fachen konjugiert-komplexen Paares. Man wird also stets neben der Hauptiterationsreihe die 
erste Nebeniterationsreihe bilden und die Iteration abbrechen, sobald eine der beiden Reihen 
„steht“. Nur in Ausnahmefällen wird man weitere Nebenreihen bilden müssen. 
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Über die erzwungenen Schwingungen in einem Kondensatorkreise, 
der eine Eisenkernspule enthält 


Von @. Colombo in Padua 


Es ist aus der Erfahrung bekannt, daß, wenn man in einem Kondensatorkreis, der eine 
Eisenkernspule enthält, eine periodische elektromotorische Kraft eingeschaltet, man subhar- 
monische Schwingungen erhalten kann. = 

Diese subharmonischen Schwingungen lassen sich erklären wenn man die starke nicht- 
lineare Abhängigkeit des magnetischen Flusses des Eisenkernes von der Stromstärke :, die in der 
Spule kreist, in Betracht nimmt. Die elektrische Ladung x des Kondensator erfüllt die folgende 


nicht-lineare D. G. 


LEERE + ge =-EN Be, u Zt A 

Wir machen folgende Voraussetzungen: 

a) Die Funktionen &(&), E(t) seien stetig nebst ihren ersten Ableitungen; überdies sei 
2(&) eine gerade, E(t) eine ungerade Funktion; E (t) ist periodisch mit der Periode 2 T. 

b) E, %,, 8,0, 6,x,I, seien positive Konstanten, &,&; gerade Funktionen von &, 71, Ns 
Funktionen von £, welche der Gleichung 3 (t + T) = — n, (£) genügen, ferner setzen wir 
kp &Gtzxa für 2]>15,+e [E+rm, für ö<t<T—8 

2) = . 5 \ 
2 %,+x, für &l<1,—o I-E+rn für T+ö5<t<2T—Ö 
(2). 
%, 0, 6 seien genügend klein, während E, 9,,2,6, R,/,, günstigen qualitativen Bedingungen 
genügen sollen. 

c) X(&) sei monoton für ,+o<|x|<I,+ 0, E (t) sei monoton für O<t<d, T—Ö 
<t<T +6, usw. Unter diesen Voraussetzungen beweise ich, daß eine periodische Lösung 
von kleinster Periode 6 7 existiert und daß diese stabil ist in dem Sinne, daß die Lösungen 
der D. G.!), wenn ihre Anfangsbedingungen genügend nahe den Anfangsbedingungen der 
periodischen Lösung sind, sich nicht von der periodischen Lösung mehr als eine bestimmte Größe 
entfernen. 


E(t) = 


Symmetrisierbare Eigsenwertaufgaben bei gewöhnlichen linearen 
Differentialgleichungen 


Von H.-K. Dettmar 
Max-Planck-Institut für Physik in Göttingen 


Die symmetrisierbaren Eigenwertaufgaben bei gewöhnlichen linearen Differentialgleichungen 
sind eine Verallgemeinerung der selbstadjungverten Aufgaben. Die Selbstadjungiertheit (vgl. [1], 
[2), wird durch die weniger einschränkende Voraussetzung der Symmetrisierbarkeit ersetzt. 
Zugrunde liegt die Eigenwertaufgabe, die aus der Differentialgleichung 


My] = AN WW]au0r Bee 
und den Randbedingungen 


U. =0.w= 1,00 2020: (2) 
besteht. Dabei ist 


Miy= 22h) y9@, NY) = 942) ya) 
und E nn. 
m—1 
U.lyl =. [xy a) + Bun y®(b)]. 
Die f, und g, sindina <& <b gegebene, reelle stetige Funktionen mit f,, +0 und mE0. Es 


ist0 <n < m. Die «&,;, und ß,; sind reelle Koeffizienten. Die Gleich 
unabhängig vorausgesetzt. r ungen (2) werden linear 


1) Die Arbeit wird in Rend. Sem. Mat. Univ. Padova Vol. XXIII. veröffentlicht. 


a Ba Zn a San Da Ua Tau ı AM len ana Ale us nu ul u Inu 25 
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1. Zusammenstellung der verwendeten Begriffe 


Mit 3 wird der reelle lineare Funktionenraum, dessen Elementeina <xz <b m-mal stetig 
differenzierbar sind und die Randbedingungen 2) erfüllen, bezeichnet. Eine Funktion ve 
heißt Vergleichsfunktion. 


Ferner werden lineare homogene Differential-Integral-Ausdrücke ®[y] betrachtet, die von 
der Form 


Z,p0(@) y®(&) (0 <m <m), 
fan & pr1(e 2) ya) O<m<m—),!......® 


T %, 702 
fd, [ds > pro 2.) da) VO <m<m—i) 
€ Cz | 
ae arlISiWE 


sind. Dabei sind die Funktionen 94, Pr» -..ina <x <bin allen Argumenten reelle stetige 
Funktionen mit p,,e =0. Für die ©, gilta <c, <b. Die höchste in ®[y] auftretende Ablei- 


tungsordnung z nennen wir die Ordnung des Ausdrucks ®[y]. 


Unter Verzicht auf größte Allgemeinheit werden die symmetrisverbaren Eigenwertaufgaben 
wie folgt definiert (zum Begriff der Symmetrisierbarkeit vgl. auch [3], [4], [5])): Die Eigenwert- 
aufgabe (1), (2) heißt symmetrisverbar, wenn ein linearer homogener Differential-Integral-Aus- 
druck Ö© [y] existiert, der ein Ausdruck der Gestalt (3) oder auch eine Summe von endlich vielen 
derartigen Ausdrücken ist, und wenn mit zwei beliebigen Vergleichsfunktionen u, v für die Formen 


{u,v} = for M |[v] dx und <u, v) = [ Ofu] No] dx 


die Symmetriebeziehungen _ 
{u, vo} = {p, u} und ZU, v> = vu) 


bestehen. Den Operator ®, der den Ausdruck © [y] erzeugt, nennen wir dann einen Vergleichs- 
operator der Eigenwertaufgabe. Sind überdies die beiden Formen positiv definit in ® ({u, u} > 0 
und <w, uy > 0 füru=0€%), so heißt die Eigenwertaufgabe volldefinit symmetrisierbar und 
Ö ein volldefiniter Vergleichsoperator. 

Jede selbstadjungierte Eigenwertaufgabe ist symmetrisierbar (®Ö= 1). Darüber hinaus gibt es symme- 
trisierbare Aufgaben, die nicht selbstadjungiert sind (s. Beispiel 2). 

Für viele Betrachtungen ist es vorteilhaft, neben ® den umfassenderen Funktionenraum ®y >% heran- 
zuziehen, dessen Elemente im Grundintervall n-mal [n = Max (x, n)] stetig differenzierbar sind und so viele der 
Randbedingungen (2) erfüllen, daß< u, vyin Qy symmetrisch und positiv definitist. (In®yhat<u,v)alle Eigen- 
schaften eines skalaren Produktes und kann zur Festlegung einer Metrik dienen.) i 


2. Sätze über die volldefiniten symmetrisierbaren Eigenwertaufgaben 

Die Eigenwertaufgabe wird im folgenden durchweg volldefinit symmetrisierbar vorausgesetzt 
mit ® als volldefinitem Vergleichsoperator. Es gelten u.a. folgende Sätze: Die Eigenwertauf- 
gabe besitzt abzählbar unendlich viele Eigenwerte ohne Häufungspunkt im Endlichen. Alle 
Eigenwerte sind positiv und lassen sich in einer Folge 0 <A, <A, <... anordnen (jeder Eigen- 
wert ist entsprechend seiner Vielfachheit gezählt). Zu dieser Folge gibt es einen Satz von zuge- 
hörigen ne Yı» Ya» «+»... Die Eigenfunktionen können orthonormiert angenommen 
werden, d.h. <%;, ya? = Öir. Ist R[u] = {w, ul/<u, u) der „mit ® gebildete Rayleighsche 
Quotient‘, so gilt 


M=MinRlul 2.2222 nn0 


wo u alle Funktionen = 0 aus ® durchläuft, die den Gleichungen <u, y >= 0(o = 1,...,k—]) 
genügen. y; ist Lösung dieser Variationsaufgabe. 
Für die Greensche Resoivente @(x, &,%) der Eigenwertaufgabe gilt die Symmetriebeziehung 
DRS = Dale]: See rs) 


in den beiden Dreiecksbereichen a <r <E<bunda<£i<sr<o (das x bzw. &am Operator © zeigt an, für 
welches Argument ® gebraucht wird). 


Ferner gelten folgende Relationen: Für u E®x die Ungleichung 


IK U, YED- KBieisisterlischerUnnleichüne)ege 2 an a we (6), 
Kt 
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und für ueQ®x und vER die Gleichung 


<Uu,v) = 2a OR Dre en Bra alt note (7). 


Genügt der Vergleichsoperator gewissen zusätzlichen Bedingungen, so ist die iterierte Funktion 


b 
q(@)= S@(®, 8,0) Ne [u(S)] ds 
a 
einer Vergleichsfunktion u(x) entwickelbar: 
© 1 
(2) = 2 oryn(@), wo = 9, yRr = 7 KU, YR2. 
k=1 k 


Die Reihe konvergiert mitsamt den ersten m Ableitungen absolut und gleichmäßig ina <x<<b. (Weitere Ein- 
zelheiten und Sätze finden sich in [6)). 


3. Verfahren zur praktischen Lösung 


Die bekannten Verfahren zur genäherten Lösung von Eigenwertaufgaben (Verfahren der 
schrittweisen Näherungen, Ritzsches Verfahren und andere) sowie die bekannten Einschlie- 
Bungssätze (vgl. etwa [2]) lassen sich auf die symmetrisierbaren Eigenwertaufgaben anwenden und 
näher analysieren. Es sei dies andeutungsweise am Beispiel des Verfahrens der schrittweisen 
Näherungen gezeigt. 


Die Vorschrift dieses Verfahrens lautet: Von einer Funktion g=E0E€ 2x ausgehend werden Funktionen 
4» 92 -- EV bestimmt als Lösungen der Randwertaufgaben 


M [4] —-tNil= N lg); U,[gi) = 0 Ver e)e :- u er (8). 
Dabei ist t reell und =Ax (t= Iterationszentrum). Mit den g; werden die Schwarzschen Konstanten 
a, = (gi, »—i) und mit diesen de Sch warzschen Quotienten w+1= a,/a,+1 gebildet. Die Schwarz- 
schen Konstanten und damit auch die Quotienten hängen von dem gewählten Vergleichsoperator ®ab. Hierin 
kommt die Verallgemeinerung gegenüber der Anwendung des Verfahrens bei selbstadjungierten Aufgaben zum 
Ausdruck. 
Genau wie bei den selbstadjungierten Aufgaben läßt sich zeigen, daß lim 1, = Ag —t = A, gilt, wo A, der 
v>0 
zu tnächste Eigenwert mit<g,,%s> = Dist. Konvergenzbeschleunigung kann erzielt werden, wenn von Schritt zu 
Schritt das Iterationszentrum t verbessert wird, d.h. näher zum gesuchten Eigenwert A, gelegt wird. 

Überdies gilt der Einschließungssatz von Kr yloff — Bogoliubov [7]fürsymmetrisierbare Eigen- 
wertaufgaben. Der Satz besagt, daß zwischen u, — us (u — 4,) und Bst Vn, (u} — 4,) mindestens ein Eigen- 
wert4-=/r,- liegt. (Zum Beweis: Mit (99, yk> = 6% wird (gi, YkD = Cr/A und mittels (6) und (7) zeigt man 

ı 


Baltı — Be) = er (a,u3— 2a,4, + Q,) == (or /Ar)? (iz An)? &(en/in)® > (n,—7,)?, wenn A, der zu u, +1 


nächstgelegene Eigenwert ist). 

Zur Durchführung eines Ritzschen Verfahrens können die Extremaleigenschaften des mit ® gebildeten 
Rayleighschen Quotienten R [u] oder auch die des mit ® gebildeten Tem pleschen Quotienten 7; [w] 
= iM, (u, — T)/(n, — r) (u = gu) herangezogen werden. Das Verfahren gestaltet sich genauso wie bei den selbst- 
adjungierten Problemen. Die dort geltenden Sätze behalten bei den symmetrisierbaren Aufgaben ihre Gültigkeit. 
Der - m plesche Quotient kann ebenso wie dort zur Einschließung von Eigenwerten herangezogen werden 
(vgl. [8]). 

Zur Illustration seien zwei Beispiele vorgeführt. 

Beispiel]. 
MY 20ya—Ayt eye El) 
Der erste Eigenwert soll näherungsweise bestimmt werden (vgl. [2], S. 157 und S. 199). 
Diese Aufgabe besitzt u.a. die volldefiniten Vergleichsoperatoren 


2 Rt 0 
Era ze 
e 3 $ ’ +1 
(Ein weiterer Vergleichsoperator ist z. B. durch ®[y]= — [ p(x,&) y(&)dE, wo p(x,&) die Greensche 
—1 


Funktion zum Differentialausdruck Z[y]= y’’ und den Randbedingungen y(+1) = Oist, definiert). 
Zur Durchführung des Verfahrens der schrittweisen Näherungen werden zwei Funktionen besti 
die mit 6= Oin der Beziehung (8) stehen. Wählt man als Ausgangsfunktion q,(@) = 61 — 1 N 


dann wird g,(2) = 644, — 790 x2 + 155 x — 10 26 +8. 


In der folgenden Tabelle sind die Werte von 4», die sich mit g,, g, bei Verwend = 
ergeben, gegenübergestellt. ; 2 ma 


1 wer Fe ee 
Be 10,574 | 10,576 | 105731408 
fs >, 10,5752 | 10,5733 | 10,5731 20 
KB 10,3896 | 10,5186 |  10,5559 14 


ah 


AS Zum al nah a a ad a u A a 
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In der letzten Zeile der Tabellesind dienach Kr yloff-Bogoliubov berechneten unteren Schran- 
ken angegeben; u, und w,liefern obere Schranken. Die exakten Eigenwerte der Aufgabe sind gegeben durch 
Ar = (20 4(k n/2)*)/(k n/2)? (k= 1,2,...). Der erste Eigenwert ist A, — 10,5730 96. 

Beispiel 2. 

a a a 9] 

y—u= ri deyn)>t.: 
Hier ist ö= N ein volldefiniter Vergleichsoperator. Zur Bestimmung des ersten Eigenwertes wurde von 
90= x (l + x)? EB ausgehend durch numerische Integration die erste Iterierte g EB (M [q,]= N [go] ‚ Iterations- 
zentrum i = 0) berechnet. Als Resultat ergabsich: u,—= 24,828, u,—= 24,665 und die Schranken (Kr yloff- 


Bogoliubov) 
22,660 <A, < 24,665. 


Zu besseren Schranken gelangt man, wenn man anstelle des Iterationszentrums {= 0 das Zentrum t,— R[g)] 
= 25,125 benutzt und eine Funktion gf€®, die mit q,in der Beziehung M [g*] —t, N [g%] = N [gu] steht, be- 
stimmt. Die Rechnunglieferte: „,—= — 1,180, #,—= — 1,169, undnach Kryloff-Bogoliubov ergeben 


sich die Schranken — 1,289 <A, < — 1,049 bzw. 
23,836 <A, < 24,076. 


(Die Rechnungen wurden auf der Göttinger elektronischen Rechenmaschine @1 durchgeführt.) 


Zusammenfassend läßt sich sagen: Durch die Betrachtung der symmetrisierbaren Eigenwertaufgaben wird 
eine über die selbstadjungierten Probleme hinausgehende, umfassendere Klasse von Eigenwertproblemen der 
praktischen Behandlung zugänglich gemacht, falls sich für diese Vergleichsoperatoren angeben lassen. Für die 
Behandlung selbstadjungierter Aufgaben ergibt sich insofern ein neuer Gesichtspunkt, als man in den Fällen, in 
denen außer dem Einheitsoperator (® = 1) weitere Vergleichsoperatoren angegeben werden können, Freiheitin der 
Wahl des Operators für die numerische Rechnung besitzt. 
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Funktionentheoretische Methoden bei Übertrasungs- und 
Regelungsproblemen 


Von Johannes Dörr in Darmstadt *) 


Die Eigenschaften eines geplanten elektrischen Übertragungssystems lassen sich in groben 
Zügen durch Wahl des Betrages der Übertragungsfunktion |Z(w)| festlegen. Dabei muß Z (w) 
eine gebrochen rationale Funktion sein, die nicht mehr Nullstellen als Pole hat und deren Pole 
alle in der oberen Halbebene liegen. Außerdem muß für die reelle Achse gelten 


IZ(o)| = |Z(— w)|, arcZ(w) = — arcZ(— wo). 
Nur unter diesen Voraussetzungen ist das System mit endlich vielen Schaltelementen realisierbar 


und besitzt stabile Übertragungseigenschaften. 
Im Einzelnen lassen sich die Übertragungseigenschaften erst beurteilen, wenn man zu typi- 
schen Ursachen «(t) die zugehörigen Wirkungen U(t) berechnet. Gibt man sich die Ursache als 


Fourierintegral E 
ul) = [ Aw) e*'do . ee TE 


so gewinnt man für die Wirkung ebenfalls eine Integraldarstellung 


+” 

Es Ss ee FAKE ee rec ea 0 

Dieser Zusammenhang stellt zugleich die mathematische Definition der Übertragungsfunktion 
dar. Die vorher angegebenen Bedingungen für Z (w) folgen aus physikalischen Betrachtungen. 

Wie findet man nun eine solche regulär analytische Funktion Z(w), deren Betrag auf der 

reellen Achse einen vorgegebenen Verlauf hat oder mit diesem vorgegebenen Verlauf möglichst 


*) Institut für Praktische Mathematik der Technischen Hochschule Darmstadt. 
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gut übereinstimmt? Zur Lösung dieses sogenannten Realisierungsproblems wird hier am Beispiel 

des elektrischen Tiefpasses ein funktionentheoretisches Verfahren angegeben. 
Als Einführung wird zunächst dieButterworth-Approximation funktionentheoretisch 

hergeleitet. Durch 
1 1 d . 9 "für jo]>1 

Im 2nV yizer tele 

wird eine Funktion beschrieben, die im Inneren des Einheitskreises, also auch auf dem Stück der 


reellen Achse — 1 <o < +1 konstant, außerhalb des Einheitskreises aber Null ist. Durch 
Übertragung ins Finite gewinnen wir daraus gebrochen rationale Funktionen der Form 


1 X ee" 1 
!x(®) = Bee = Io: EM N RETTET (2a) 
mit R 
7 
Yp = @r—1)5% NT ab 
Man kann (2a) auch auf die Form bringen - 
In(o) = Ho a NE EREIRN DICH 


Für die reelle Achse hat fy(w) bei genügend großem N den für einen Tiefpaß gewünschten Verlauf, 
aber /v(w) hat auch in der unteren Halbebene Pole. Lassen wir in der Produktdarstellung (2ec) 
alle Faktoren fort, deren Pole in der unteren Halbebene liegen, so entsteht, wie man aus den 
Symmetrieeigenschaften von (2c) leicht erkennt, eine Funktion, deren Betrag auf der reellen 
Achse ebenfalls den für einen Tiefpaß gewünschten Verlauf hat. Das heißt, 
N 
); 
Bx(o) = II 
il 
ist bei genügend großem N als Übertragungsfunktion eines Tiefpasses geeignet. 

Diese Übertragungsfunktionen haben jedoch die unangenehme Eigenschaft, daß arc By(w) 
stark vom Approximationsgrad N abhängt. Die Phase arcZ(®) beeinflußt aber ebenfalls die 
Übertragungseigenschaften. Man muß deshalb bei der Realisierung auch den Phasenverlauf in 
gewissen Grenzen beherrschen. 

Um zu Übertragungsfunktionen zu gelangen, die sowohl in ihrem Betrag als auch in ihrer 
Phase festgelegt sind, gehen wir von folgenden speziellen Funktionen aus: 


2.) = | ei u) mine bo e 


Eine Untersuchung von |Z,(o)| und arc Z,(o) für die reelle Achse zeigt, daß diese Funktionen als 
Übertragungsfunktionen für Tiefpässe hervorragend geeignet sind!). Zur Realisierung muß man 
aber die Z,(w) durch gebrochen rationale Funktionen approximieren. Zur Bewältigung dieser 
Aufgabe suchen wir nach dem Vorbild des vorhergehenden Beispiels nach einer Integraldar- 
stellung 


@— ei® 


1.) 
Zu) = N rd er ß 
= (a (3b) 
Wenn eine solche Integraldarstellung existiert, so liefert die Übertragung ins Finite sofort eine 
Approximation durch gebrochen rationale Funktionen. 
Wir gehen aus von der Cauchy schen Integralformel in der speziellen Form 


1 Pen 1 
Fo) = 5, 8 Fl | a) lie ER EERTFT Ea: 


Dabei sollen die Punkte v — z und v = O im Innern des vom Weg s umschlossenen Gebietes liegen, 
und der Weg s soll im abgeschlossenen Bereich eines durch die Punkte + laufenden Kreises 
liegen. Auf s und in dem von s umschlossenen Gebiet soll F(v) regulär analytisch sein. 

Mit der Transformation 


= lot io, v= pl me ee 
3 F(2 (o)) = fo) = g(@) + h(o) VI oe RE (4c) 


ie) A.W altherund J.Dörr, Typische Einschwingvorgänge zu fundamental bert ionen 
A.E.Ü. 7, 379386 (1953). ne © Ühersca gunestupl siegen 


a u a A u Lat an Hm a LE aD Lu Lu. 
‘ 
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entsteht nach einigen Umformungen aus (4a) die Darstellung 


1 A | 1 Ay0® Br 
ee I ee . (da). 
Dabei geht der geschlossene Weg s in den geschlossenen Weg o über. Der Punkt o liegt aber hier 
außerhalb des von o berandeten Gebietes. Daraus folgt wegen der Zerlegung (4c): 


Aa) 


Ko) = 358 Se, 


mE fe a AR nn oh), 


Wählen wir als Weg s einen Kreis, der durch die Punkte 2 = 7 i läuft, so läßt sich dafür auch 
schreiben 


+ +1 


4 9a) 1 (aayı— | 
110) —— = | yı a dA — j ee et a (de). 
Der Integrationsweg ist dann ein beliebiges von = — 1 nach = + i geführtes Kreisbogenstück. 


Man kann aber den Weg s allgemeiner so wählen, daß sich in der Darstellung (5c) ein doppel- 
punktfreier Weg ergibt, der beliebig von —1 nach +1 läuft. 
In (5c) haben wir die gewünschte Integraldarstellung gefunden. Mit F(z) = z", also 


Io) = (1 — © — io)” = (—üR [ To) + iU,(o) Y1 — o2] 
ergibt sich 
+1 


a. 1-22 
Neger ur | 0,0) 11 —A: di 
7 A—0® 
—ı 
Die U,„(A) sind dabei die Tschebyscheffschen Polynome zweiter Art. Denken wir uns in 
der letzten Integraldarstellung den Integrationsweg im Halbkreis von —1 nach +1 durch die 
obere Halbebene geführt, so liefert die Übertragung ins Finite die realisierbare Approximation 


1 i > BR sin @,U„(e!®) ei 3/2 9, 
Ze a ee 


Rn, 
v=iI © a 


ER: eb) 
mit s 
It c 
Bon (2v—1). 


Diese Reihen ergeben auch bei relativ wenig Summengliedern N recht gute Approximationen. 
Wenn deshalb N nicht extrem klein gewählt wird, darf man bei der Berechnung der Wirkungen 
nach Integral (1b) als Übertragungsfunktion den Ausdruck (3a) an Stelle des unhandlicheren 
Summenausdrucks (6) einsetzen?). Dabei treten, wenn man als Ursachen Sinusquadratfunktionen 
wählt, Integrale über Produkte von Besselfunktionen mit Exponentialfunktionen auf ?). Nur 
für die Realisierung der Übertragungsfunktion in Form eines Netzwerkes ist in diesem Fall 
die Approximation unbedingt erforderlich. 


Neue Verfahren zur direkten Lösung des allgemeinen 
Matrizeneigenwertproblemes 


Von Sigurd Falk in Braunschweig 


Die Umwandlung des allgemeinen Matrizeneigenwertproblemes Ar = kBr der Ordnung n 
in das spezielle Problem B-! Ar = kr erfordert rund 4/3 n? Multiplikationen, die anschließende 
Transformation von Hessenberg (oder eine ihr gleichwertige andere) etwa n? Mult., so daß 
zur Aufstellung des charakteristischen Polynomes 9 (k) = |X — k ®| rund 7/3 n? Mult. erforder- 
lich sind. Es werden zwei Verfahren entwickelt, die nur 6/3 n? Mult. benötigen und sich stark 
schematisieren lassen. Voraussetzungen werden nicht gemacht; X und ® dürfen beide singulär 
sein. 

2) Siehe !), 8.380. 


3) J. Dörr, Unbestimmte Integrale über Produkte von Besselfunktionen mit Exponentialfunktionen, 
Wiss. Zeitschrift Techn. Hochsch. Dresden 2, 353—354 (1952/53). 
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Verfahren 1: Verallgemeinerung der Methode von Hessenberg. Die: gegeben 
Matrix 8 wird mittels Zeilenkombinationen in eine untere Dreiecksmatrix B übergeführt, dabei 


geht X in Nüber. X und.-® lassen sich nun in einem geschlossenen Rechengang in 9 (Hessen- 
bergmatrix) und & (Einheitsmatrix) transformieren: wir multiplizieren die Gleichung 


Ir = kB ET ER (1) 
von links mit der nichtsingulären Matrix & und führen neue Vektoren ein: a 
“ | ey Rn ee See 
Dann wird aus (1) - £ 
SÜRH Eh LONn (3) 
Jetzt soll 3 > & 
LAR=H oder NE REH ee (4) 
und ’ R 
LEN=E oder BREITE aa TEA (5) 
sein, daraus folgt nach Elimination von I!: 
ÄRSEERE SORTE (6). 


In dieser Matrizengleichung sind X und B gegeben, R und 9 gesucht. Nun zeigt sich, daß R, 
somit auch BR eine untere Dreiecksmatrix ist, demnach hat Gl. (6) die Form: 


Hierin ist © das Produkt von B und R, & einerseits das Produkt von X und R, anderer- 
seits das Produkt von ® und ®, also it A\R = BNS, wie Gl. (6) es verlangt. Die erste Spalte 
von ® wählt man nun beliebig, (am einfachsten einen Einheitsvektor); dann lassen sich der Reihe 
nach die ersten Spalten von ®, ® und ® eindeutig berechnen, außerdem die zweite Spalte 
von @. Damit liegen dann ebenso alle zweiten Spalten von ®, ® und ® und die dritte Spalte 
von ® fest, und so fort, bis sämtliche Matrizen ausgefüllt sind. Dabei ist jedes gesuchte Element 
einzige Unbekannte in einer linearen Gleichung. In Ausnahmefällen müssen in ® und ® 
Spalten vertauscht werden, auch kann 9 zerfallen. Zum Schluß wird das Polynom |$ — k €] 
rekursiv aufgestellt, der Eigenvektor ), zum Eigenwert k, berechnet und nach (2) zurücktrans- 
formiert. 


Falls ® singulär ist, verläuft die Umwandlung in die untere Dreiecksmatrix ® in etwas 
abgeänderter Weise. 


Verfahren 2: Es wird nicht A und ® je für sich, sondern die Polynommatrix ® (k) = 
— kB + X insgesamt transformiert, wobei lediglich Zeilen, nicht Spalten von ®(k) zu kombi- 
nieren sind. Mit Hilfe des Begriffes der einer beliebigen Polynommatrix zugeordneten ‚Koeffi- 
zientenmatrix“ lassen sich die von den Matrizen mit konstanten Elementen geläufigen Rechen- 
verfahren auf Polynommatrizen übertragen. Zum Schluß bleibt im wesentlichen eine zwei- 
reihige Matrix übrig, deren Elemente Polynome in k sind, und zwar im allgemeinen vom Grade 

n 1 —1 
n/2 bei geradem n und vom Grade 2 5 bzw. m bei ungeradem n. Auch Aufgaben der 
At (PUALHRBH+ELT=0, AU+HKB +) = 0 usw. lassen sich ohne weiteres behandeln. 
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Sämtliche Eigenvektoren sind bei berechnetem Eigenwert unabhängig voneinander leicht zu 
gewinnen. 


Hauptveröffentlichung mit mehreren Zahlenbeispielen: „Abhandlungen der Braunschwei- 
gischen Wissenschaftlichen Gesellschaft, Band VI, 1954“. Verlag Friedr. Vieweg und Sohn, 
Braunschweig 1954. Dort ist das Verfahren 1 mit II, Verfahren 2 mit III bezeichnet. 


Eine Methode zur Lösung von Anfangswertproblemen bei gewöhnlichen 
und partiellen linearen Differentialgleichungen zweiter Ordnung 
Von Helmut Grunsky in Mainz 


a) Das Wesen der Methode wird bei einer gewöhnlichen linearen Differentialgleichung 
1. Ordnung am deutlichsten und der äußerst. einfache und übersichtliche Weg zu bekannten 
Resultaten kann als Vorbild für die höheren Fälle dienen. 


Gegeben sei die Differentialgleichung y’ + a, (x) y = 0; gesucht ist die Lösung y(a) mit 
Ya) = Yo. 
Wenn 1, (£) eine Stammfunktion von a, (£) bezeichnet, so ist: 


zT z T T 
YWE,=[yOE=-—- [uayE-— hy + fhyie=—lhye —[haryde. 
% as do ur 
Setzt man nun |, a, = q,, so läßt sich die partielle Integration iterieren und man erhält schließlich: 


[y]z, m (2 ı) y\ = An+ı Y dE ’ = 4,— l,_ı 4,» 
= % = 


Die 1,(£) können noch in verschiedener Weise normiert werden. Wählt man stets 1,(x) = 0, so 
strebt das Restintegral gegen 0 mit n— oo und man erhält 


Ye) = y(a)- exp -/ a, (&) ae), 


wobei exp {: - } als Abkürzung für die e-Reihe steht. — Daß auch eine andere Normierung von 
Interesse sein kann, sieht man an dem speziellen Fall a,(«) = 1. Wählt man hier die |, als 
dieBernoullischen Polynome für das Intervall x, <& <x, so erhält man mith= 2 — 
die bekannte Beziehung ‘ 


oo B, oo B, 
Zurem)o-(Zir)um. 
wo die B, die Bernoullischen Zahlen sind. 
b) Liegt für die Differentialgleichung 2. Ordnung 
EAN NEE En Ve re (1) 


das Anfangswertproblem y(2,) = Yo, Y'(%) = Y, Vor, so gehe man von der Formel aus: 


* 


ya)=y+ty (2) +[@-8 y(öde. 


Ersetzt man hier y’’(£) vermöge (1) durch y’ und y und formt die dann auftretenden Integrale 
durch partielle Integration um, indem man wieder auf y’’ hinstrebt, so erhält man mit den 
Bezeichnungen: 


AHd=al) =R— Hal), BO=bH) = R—H5)), 
y(z) = Yo + Y0° (8 — %) — [Bı y]z, + [3 — 4ı) 7], hr ge: —B)Y de. 


Verfährt man mit dem Restintegral ebenso, so kommt nach n-maliger Iteration: 


voa=n+w en [Zain], + [ (An— B,) y” de 
@ Be nanBene as 022: 
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Bei der noch freien Normierung der A,(£) und B,(£) hat man darauf zu achten, daß y'(x) nicht 
18 . * 
auftreten darf, d.h. es ist zu fordern 3 (B, (2) — A,(«)) =0. Natürlich wird man i. a. schon 
‚—1 


das einzelne Glied dieser Summe gleich 0 wählen. Es ist leicht zu sehen, daß der Prozeß sicher 
dann konvergent ausfällt, wenn man ferner B,(x) = 0 setzt. Hat man es mit dem Fall a (x) = 0 
zu tun, wenn man ferner B,(x) = 0 setzt. Hat man es mit dem Fall a(x) = 0 zu tun, so sind 
alle A, = 0 und für B, hat man die günstige Abschätzung: 


M?- |® 3 ger 
Be Ta 
wo M eine Schranke für |b(&£)| in x, SE < x bedeutet. 


In gewissen Fällen kann eine Modifikation des Verfahrens zweckmäßig sein; liegt z. B. die 
Besselsche Differentialgleichung vor 5; 


zyt Hy FREI NER ERNEREHRIENE (2) 
so gehe man aus von 


[evr@a = Ey @ — vo, -— [WO+EyO)@. 


Nun strebt man bei der partiellen Integration immer auf £y’’(£) hin und ersetzt diesen Ausdruck 
vermöge (2) durch — (y’ + & y). Man erreicht dadurch, daß als bekannte Funktionen, mit deren 
Hilfe die Lösung ausgedrückt wird, nur Polynome auftreten. 


c) Liegt eine hyperbolische partielle Differentialgleichung 2. Ordnung in der Normalform 
vor: 


uy +) +) tra, )u—0, 
so wollen wir uns (der Einfachheit halber) das charakteristische Anfangswertproblem stellen: 
u=ia(2)‘ Ey) w=Bly)- für ge. 
Mit R sei das Rechteck ©, SE <x, y <n Sy bezeichnet, mit R sein Rand. Wir gehen aus 


von der Formel: 
[uln = LJ Uny(&,n)d&dn 
in der [u] die Wechselsumme der Funktionswerte von u in den Ecken von ® bezeichnet. Nun 


folgt durch partielle Integration, wenn /,(£, n) eine Funktion mit 


en) un: +) — al) 
bedeutet: 


lula = — I (a u: + 51%, + cu) dE dn 
Damit in | nur Bekanntes auftritt, ist /,(&, 7) so zu normieren, daß 
R 


h:&YJ)=b(& 4%); nam) =-a@n). 
Damit ist l, bis auf eine additive Konstante festgelegt. Das erste Integral der rechten Seite ist 


— [hu]a, also teils durch Bekanntes, teils durch die gesuchte Größe u(x, y) ausgedrückt. 
Iteration führt schließlich zu 


a Y 
(21), - [228 w)aoa + (24 m)BoDan + [| Imundsan 
Te = Yo BE R 

mit 
In;en = An;g + Dn;n — On» In;&(&, )=bn(E N) In;n(& n) =a,(2,n); el 5 
N An ae Ines B, =bn — Imst nr. 


Falls a, = b, = 0 ist und alle I, so normiert werden, daß l„(2, %y) = 0 (in welchem Falle übrigens 
das Verfahren mit dem Picardschen Iterationsverfahren zusammenfällt) findet man 


ME 
n!? 2 


In &; n)| s 


wo M eine Schranke für |e| in A bedeutet. 
Auch hier können gewisse Modifikationen wie bei b) zweckmäßig sein. 
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Zur Nomographierbarkeit höherer Funktionen 
Von W. Händler in Hamburg 


Als neuere Entwicklungsrichtung in der Nomographie kann die systematische Unter- 
suchung der Nomographierbarkeit von bestimmten Funktionenklassen gelten. Funktionen 
komplexer Veränderlicher F(w,z) =0 wurden insbesondere von I. A. Willner (z.B. 
Mathematische Zeitschrift Maremarngecrnä Cöopnus T27(69) N.1 1950 S. 3—46 (russ.)) 
untersucht. Dabei soll die Lösung als eine Fluchtgerade in einem Nomogramm mit 4 Funk- 
- tionsleitern dargestellt werden. (a,b,9,gq; wenn z=a-+tib und w=p-üig). Läßt man 
zu, daß bis zu zwei der Funktionsleitern krummlinig sein können, so ist als Bedingung für die 
Nomographierbarkeit notwendig und hinreichend, daß die folgenden Differentialparameter mit 


de 
de 

Be N 
AH 1 --lkum), ne). 


konstant und reell sind. Sämtliche nomographierbaren F unktionen können im übrigen aus der 
Integration einer der beiden folgenden Differentialgleichungen gewonnen werden: 


WW I WAT, WE (2), 
Mi Eee EN URL ale 8 Rn En (3% 


Lösungen von (2) bzw. (3) sind neben einer Reihe bekannter elementarer Funktionen die ellip- 
tischen Integrale 
ve) 


w—Ww=M ee ; wor aM y2reellis 2 2 (4) 
— k? sin 
bzw. 
Nk v(e—-%) de 7 
j Ü 

ee BE RLE SR 2 gr Kr 

u—mw=ı ne: | ee ‚„ worinN, y?reell und + IM (5). 
r } 


Von diesen beiden Grundformen gelangt man auf dem Wege über bekannte Beziehungen z.B. 
zwischen den elliptischen Funktionen Jakobis, Weierstraß zueiner Reihe von anderen 
Nomogrammdarstellungen. Für einen festen Modul % = konstant erhält man zunächst ein 
Nomogramm mit 2 geradlinigen Leitern für « und 5 (nach Voraussetzung) und mit 2 krummlinigen 
Leitern für p und q in Form eines Kegelschnittes. Variiert man.den Modul k, so erhält man ein 
ganzes Kegelschnittbüschel und kann mit einem einzigen Nomogramm den ganzen Wertevorrat 
einer elliptischen Funktion, eines elliptischen Integrals, bzw. des Integrals oder Logarithmus einer 
elliptischen Funktion in Abhängigkeit von Modul %k übersehen. 


Der Entwurf eines Nomogramms führt über die sogenannte ‚„kanonische Darstellung‘‘, 
nach der alle Größen in reeller Form als Gleichungssystem geschrieben werden. Die „kanonische 
Darstellung‘* lautet etwa für (4) mit reellem y: 


2 
+1) 0 [2ya—a)l + BE Co [27 — dp] + NE 0 
3 
an en see., fe 
PrPr : b Pr Br 
Darin sind 

b = on Ber; = on Re ;2. = dn|Re" 7; a], hr 
Bi = 8n re. Ba = en 1 m 3 e|, ß; = dn Im? 2 e| 


(sn,en,dn Jacobische Elliptische Funktionen). 
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Zum Zwecke der geometrischen Interpretation bildet man aus (6) eine 4-reihige Matrix vom | 
Range 2: 


er 01 03 — 03 \ | 
i Pu Fu Be | 
3 3 3 | 
Kun, 1—k': pi | 
EN ee | 
2. Ps 2 Ps | 
1 0 
0 I &f2y6—) 


Die Struktur eines Nomogramms liegt damit bis auf Kollineationen fest. Das Nomogramm hat 
im einfachsten Fall die Gestalt 


[ 2 k2a?o2 Ne 

al —az 0a — 

ke? p% 2 B,2 i 
rag EN SR 7 

1— kB 1— kB .. (9). 
! _ i 
cos2y(a—a,) 
0 ! 1 


| &f2y0—b,) 
Die Leitern der Veränderlichen 2, bzw. a und b, sind geradlinig, wie vorausgesetzt und aus der 
Matrix zu ersehen ist (3. und 4. Zeile von (9)). Die Leitern für die Veränderliche w, bzw. p und g, 
liegen auf Kegelschnitten (1. und 2. Zeile von (9)). Verändert man den Wert des Moduls %, so 
erhält man, wie bereits erwähnt, ein ganzes Kegelschnittbüschel. 


Die Punktdarstellung einer Matrix nach Konrad Friedrich und 
sein Gitterdeterminanten-Verfahren; Anwendungen und Methodisches 


Von Werner Jenne in Frankfurt am Main 


Es ist der Zweck der folgenden Zeilen 


1. eine Analogie hervorzuheben, die meiner Ansicht nach besteht zwischen dem Hornerschen 
Verfahren zur Berechnung des numerischen Wertes eines Polynoms einer Veränderlichen 
einerseitsunddem Gitterdeterminanten-Verfahren von Konrad Friedrich(1)(2) 
andererseits; 

. die von K. Friedrich 1930 angegebene ‚„Punktdarstellung eines Normalsystems‘‘ (3) 
auf den allgemeinsten Fall der Punktdarstellung einer Matrix oder einer Determinante aus- 
zudehnen, der denkbar ist; 

3 die methodischen Eigenschaften dieser Darstellung kurz zu schildern; 

4. die Anwendungen und Anwendungsmöglichkeiten in Theorie und Praxis anzudeuten. 


D 


1. Analogie Horner-Friedrich 

Bekamntlich schuf Horner sein so überaus einfaches Verfahren zur Berechnung des 
numerischen Wertes eines beliebigen Polynoms einer Veränderlichen deshalb, weil bei der sich 
direkt an die übliche Schreibweise der Polynome anschließenden Berechnungsart die numerische 
Berechnung der auftretenden hohen Potenzen der Veränderlichen zu mühsam ist. 

Ähnlich liegt aber der Fall, wenn der numerische Wert einer Determinante berechnet 
werden soll, bei der die Anzahl der Zeilen gleich 5 oder größer ist. Es gibt bekanntlich zwei immer 
anwendbare klassische Verfahren zur Berechnung des numerischen Wertes einer Determinante: 
erstens durch Summation zahlreicher Produkte aus n Faktoren (n = Reihenanzahl), zweitens 
Entwickelung nach den Elementen einer beliebigen Zeile oder Spalte. Beide Verfahren werden mit 
wachsender Zeilenanzahl der zu berechnenden Determinante rasch unübersichtlich und mühsam 
in der Handhabung. Das Friedrichsche Gitterdeterminanten-Verfahren vermeidet ver- 
schiedene Nachteile, die den beiden klassischen Berechnungs-Verfahren anhaften, durch Ein- 
führung einer bezifferten ebenen Abbildung der Determinante, anhand derer nach .einfachen 
Vorschriften gerechnet wird, die Friedrich aus der Theorie der Determinanten abgeleitet 
hat, veröffentlicht zuerst in (1), später ausführlicher in (2). Allerdings zeigt sich das Gitter- 
determinanten-Verfahren im allgemeinen nur dann den klassischen Verfahren überlegen, wenn 
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die zu berechnende Determinante zahlreiche verschwindende Elemente aufweist (Gitterdeter- 
minante). Ich möchte an dieser Stelle der heute merkwürdigerweise häufig anzutreffenden An- 
sicht gegenübertreten, daß nur solche Verfahren zur Auflösung linearer Gleichungssysteme von 
Interesse seien, die bei jedem Gleichungssystem anwendbar sind, ohne Rücksicht auf besondere 
Eigenschaften der Matrix des Systems. Ganz im Gegenteil, neben diesen etwa 20 Verfahren 
brauchen wir auch solche, welche die Möglichkeit bieten, besondere Eigenschaften der Matrix 
des Systems zu berücksichtigen. Unter diesen ist wohl das Friedrich sche z 
minanten-Verfahren ohne Zweifel das allgemeinste, welches bisher bekannt wurde. 


2. Die Punktdarstellung einer Matrix oder einer Determinante 

Abbildung einer beliebigen Matrix in der Zeichenebene nach folgender Vorschrift: die Ele- 
mente der Hauptdiagonale werden bezifferten Punkten der Ebene zugeordnet, die gemischten 
Elemente dagegen ebenfalls bezifferten Kurvenstücken, die bei beliebigem Verlauf, beliebiger 
Länge und mit beliebiger Krümmung das zugehörige Punktepaar verbinden. Meist können hier- 
für Geradenstücke gewählt werden. Verschwindende Elemente: Kein Abbild. Unsymme- 
trische Matrizen oder Determinanten: Doppelbezifferung der Kurvenstücke, je nach dem 
Durchlaufungssinn. Nichtquadratische Matrizen: freiendende Kurvenstücke. Symmetrische 
Matrizen oder Determinanten: Durchlaufungssinn der Kurvenstücke kann unberücksichtigt 
bleiben (Fortfall der Doppelbezifferung). Konstanter Wert aller von Null verschiedenen gemisch- 
ten Elemente: Ersparnis der Kurvenbezifferung (häufiger Fall auf verschiedenen Anwendungs- 
gebieten). Die Elemente der Hauptdiagonale haben einen festen Wert: Ersparnis der Punkt- 


— a2 : 
a] -— 
92 


q 
1 
Bild 1. Ausschnitt aus der allgemeinsten Punktdarstellung eines linearen Gleichungssystems 


bezifferung, daher nun rein geometrische Abbildung, frei von jeder Metrik, an kein Koordinaten- 
system gebunden, Krümmung der Kurvenstücke gleichgültig, entscheidend sind nur Lage- 
beziehungen und Symmetrieeigenschaften der Bildfiguren. Zum Beziffern der Bilder können 
auch Matrizen, Tensoren, Gebilde der logistischen Mathematik und dergleichen verwendet 
werden. Vgl. Bild1. 


3. Die methodischen Eigenschaften der Punktdarstellung 
In methodischer Hinsicht besteht wohl am meisten Verwandtschaft zu den bezifferten 
Kräfteplänen, die bei der Berechnung statisch unbestimmter Tragwerke Verwendung finden, 
z.B. zum Cremonaschen Kräfteplan. 
Es wäre noch festzustellen, wie weit die naturgemäß auftretenden topologischen Probleme 
aus der Theorie der Gelenkwerke (Graphen) bekannt sind oder nicht. 


4. Anwendungen und Anwendungsmöglichkeiten 

a) Geodätische Ausgleichungsrechnung. Ohne Zweifel lag bei den An- 
wendungen das Hauptgewicht bisher auf dem Gebiete der geodätischen Ausgleichungsrechnung. 
Und hier wiederum bieten die, durch Arbeiten von mir nicht unwesentlich weitergetriebenen 
Friedrichschen Verfahren da besondere Vorteile gegenüber den klassischen Methoden, wo 
das Boltzsche Entwickelungs-Verfahren oder sein Substitutions-Verfahren zur Anwendung 
gelangen. Dies ist überwiegend der Fallin den Ländern deutschen, italienischen und französischen 
Sprachgebiets; dementsprechend sehen wir in diesen Ländern eine durch zahlreiche Abhandlungen 
belegte lebhafte Beschäftigung mit den Friedrichschen Verfahren, vgl. die Literaturüber- 
sicht. Ein besonders tiefes Erfassen und Weiterführen der Ideen von Boltz, Friedrich 
und der meinigen ist in Frankreich festzustellen, und zwar speziell durch die Arbeiten von 
M. Dupuy (4) (5). In der letzterschienenen Abhandlung (5) beschäftigt sich M. Dupuy 
eingehend mit den von H. Boltz (6) und mir (7) veröffentlichten „Kehrmatrizen-Tafeln“ — 
ein, in französischer Sprache, von ihm geprägter Ausdruck — und macht interessante und zweifel- 
los sehr brauchbare Vorschläge für die weitere Ausgestaltung dieses neuartigen Rechenhilfs- 
mittels. Eine Stellungnahme von mir zu der genannten Abhandlung vonM. Dupuy ist für 
das Septemberheft 1954 des „Bulletin Geodesique“, des Organs der Assocat ion Inter- 
nationale de Ge&od&sie, Erscheinungsort Paris, zum Abdruck in englischer 
Sprache !) angenommen. Ferner habe ich einen Arbeitsvorschlag ausgearbeitet zur schnellen und 

L) Weil einerseits im englischen Sprachgebiet die Friedrichschen Arbeiten noch wenig bekannt zu 
sein scheinen, andererseits dem englischen Mathematiker J. J. Sylvester grundlegende Arbeiten über Kon- 
tinuanten zu verdanken sind. 
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rationellen Durchführung der Berechnung von geodätischen Kehrmatrizen-Tafeln mit Hilfe von 
Lochkarten; dieser Vorschlag wird, wie ich hoffe, in einer deutschen geodätischen Zeitschrift 
. erscheinen können. ALT 

b) Andere Gebiete. In seiner grundlegenden Arbeit (3), erschienen 1930, nennt 
Friedrich bereits als Anwendungsgebiete seiner Methoden: Technische Mechanik — 
Dimensionierung der statisch vielfach unbestimmten Tragwerke, und Elektrotechnik — Berech- 
nung elektrischer Leitungsnetze. 

Hierzu erwähnt er, daß sich im letzteren Falle die Punktdarstellung meist bereits mit dem 
Netzbilde decke. Verschiedene Beispiele für diese Anwendungen führt Friedrich a.a.0. 
an und durch. 

Als weitere Anwendungsgebiete, die sich bis auf III. erst nach dem Tode von K. Fried- 
rich (1944) ergeben haben, sind zu nennen: 

I. Reduzibilitäts-Eigenschaften spezieller, von Nörlund definierter Tschebyscheff- 
scher Polynome, wobei insbesondere meine Ergebnisse bezüglich der Summenspaltung und 
Differenzenspaltung bei symmetrischen Punktdarstellungen vorteilhaft zur Anwendung 
gelangen (1); . 

II. Säkulargleichung und Eigenwertprobleme (Arbeiten im Institut für Praktische Mathematik, 
Darmstadt); 
III. Randwertprobleme bei partiellen Differentialgleichungen — Poissonsche Gleichung, 

Laplacesche Gleichung, vgl. z.B. (5); 

IV. Gewisse Probleme der Theoretischen Chemie. 


Schlußbemerkung 


Daß Punktdarstellung und Gitterdeterminanten-Verfahren von K. Friedrich wert- 
volle Hilfsmittel darstellen bei der Inversion von Matrizen, speziell von Gittermatrizen, 
ergibt sich aus der angezogenen Literatur und den obigen Ausführungen. Für weniger wahr- 
scheinlich halte ich es, daß die Punktdarstellung auch bei der Multiplikation von 
Matrizen Vorteile zu bieten vermag; immerhin wäre wohl auch diese Frage zu prüfen. 
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Die Kernstrahlen in der Darstellenden Geometrie 
Von H. Lenz in München 


Die Aufgabe, aus zwei gegebenen linearen Bildern eines Gegenstandes ein drittes zu kon- 
struieren, läßt sich nach G. Hauck mit Hilfe der Kernstrahlen lösen. Die Übertragung von 
Strahlen in projektiven Büscheln, auf die es dabei ankommt, kann nach S. Finsterwalder 
mit dem Papierstreifen ausgeführt werden; ebenso die u. U. vorher nötige Auffindung dreier 
Paare entsprechender Kernstrahlen zu je zweien der drei Bilder. 


dr Beispiel: In Grund- und Aufriß sind Gegenstand und Projektionsrichtung s 
gegeben. Gesucht ist der Riß auf eine zu s senkrechte Ebene. Nach der Nyströmschen 
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Formel cos & = cos 9 siny (Bezeichnungen s. Bild 1) findet man den Winkel &£ und damit die 
Kernstrahlen zu Aufriß und drittem Riß wie folgt: 


Fälle von einem Punkt Z der 2’’-Achse das Lot Z8 auf s’’! Trage 0’S =O’T auf s’ ab! 
Zeichne den Kreis X um 0” 
mit dem Abstand (7, y‘) als 
Radius! Die Tangente t von 
Z aus an K schließt dann 
mit der y’’-Achse den Win- 
kel & ein. Die &-Achse ist 
also das Lot zu t durch O. 
Für beliebige Punkte P ist 
dann Abstand (Pr, si) = 
Abstand (P,x) und Ab- 

stand (P’, s’) = Abstand 
(P, 2). Zur Abstandsmes- 
sung und -übertragung ver- 
wende man, soweit möglich, 
den Stechzirkel. 


2. Beispiel (all- 
gemeine Axonometrie): Ge- 
geben sind Grund- und Auf- 
riß, sowie die Achsenbilder 
des Schrägbildes mit Ein- 
heitspunkten A, B, ©. Man 
findet den Grundriß s’ (und 
analog den Aufriß s’’) mit un 
dem Papierstreifen (in Bild 2 
nur für s’ eingezeichnet) auf 
Grund der Doppelverhältnis- 
gleichung 


[4 


(USER TR) 

bzw. 

[(y" u 2° 0) = (99 Bl. 
Die Übertragung der Kern- 
strahlen durch gegebene 
Punkte P erfolgt dann etwa 
mit Hilfe geeigneter schräg 
eingezeichneter Skalen m, n. 
In’Bila2 ist 40 =Ab- 
stand (A, 2) und O"R=Ab- 
stand (0, &)] 


3. Beispiel: Ge- 
geben sind Grund- und 
Aufriß. Gesucht ist ein _per- 
spektivisches Bild, in dem 
der Einheitswürfel vorge- 
geben ist. Die Kernstrahlen 
im dritten Bild gehen durch Bulasa 
die Fluchtpunkte der z-Achse 
bzw. der x-Achse. Man kann mit dem Papierstreifen je drei dieser Kernstrahlen (etwa durch 
drei Ecken des Einheitsquadrats der «-y- bzw. y-z-Ebene) in den Grund- bzw. Aufriß über- 
tragen und dann weitere Kernstrahlen, die durch gegebene Punkte gehen, aus den gegebenen 
Rissen entnehmen und in das perspektivische Bild übertragen. 
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Bemerkungen zur günstigsten projektiven Abbildung von Skalen 
und Nomogrammen 


Von Rudolf Ludwig in Braunschweig 
Eine im Bereich x, < x < x, stetige, eindeutige und monoton zunehmende F unktion F(x) 


läßt sich durch 
Fa) — Fi) 


ee (1) 
Fe) TG) As 
als geradlinige Funktionsleiter der Länge 1 darstellen. Durch eine projektive Abbildung 
va Men 2. Am Be EEE (2) 


mit 


T (9), en > und A = (ar), (di, =, 2), A=+0 
und mit der speziellen Matrix 


u-(°,) a (3) 


T 


erreicht man, daß die Randpunkte der Skala festbleiben und daß für 7 > 0 der Durchlaufungs- 
sinn der Skala erhalten bleibt. 


Die allgemeine Projektivität (2) mit 3 wesentlichen Konstanten läßt sich zerspalten in 
eine Affinität und die spezielle Projektivität (3). Zwischen 7 und den a;; besteht dann die Be- 
ziehung 


Der Teilstrichabstand ist für kleine Schrittweiten proportional f(x). Daher wird als 
Teilstrichkennzahl T definiert. 


TMHDIE= FM nn ee (5). 
Für die durch (3) transformierte Skala ist 
TE) ER an SE Ro: 


1+@— Dial 


T als Funktion von x und x ist deutbar als Fläche über dem Bereich Blmn ser su 0<7<el, 
sie hat folgende Eigenschaften: 


a) Das Produkt der Teilstrichabstände. an den Skalenenden ist unabhängig vom Ver- 
zerrungsparameter T: 


Ike): Tife)] = Ra): Pla) = const 2... ee 


b) T(&, 7) hat in® stets mindestens eine stationäre Stelle &,, zu. Aus den beiden 
Bedingungen 97/0x = 0T/ör = 0 folgt nach Elimination von 7 


en A War u 


und 


= — 1, 4,.12,% VIERTE HH 

(zo) N 
Denkt man sich beide Seiten von (8) über f(x) aufgetragen, so sieht man leicht, daß unter den 
angegebenen Bedingungen beide Kurven sich innerhalb ® wenigstens einmal schneiden müssen. 


Die Entscheidung über die Art der stationären Stelle liefern die zweiten Ableitungen und 
die Diskriminante. Man rechnet aus: 


3 ) f Fi E 3 1 4 ff f? 
RES To) ee) ’ N, T,) sum Ta) = gm) =, 
2 1 vr jrtı 2 e 
Det) = Ta In — Ta = —- 641 — pi (ey |, Sale Pi: m (11). 
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Daraus folgt für die Stelle x, 7, : 


Sara 0... DD, 


JE 
Spap 


# 5; iR 
oder 2 —3f"” 5 
a Er 
Ist D <. 0, so hat die Fläche einen Sattelpunkt. Des letzte Fall würde ein Maximum bedeuten. 
Da T am linken Rand proportional 7, am rechten Rand aber umgekehrt proportional r ist, sieht 
man geometrisch leicht ein, daß es wenigstens einen Sattelpunkt geben muß. 

Die Funktionsleiter, die dm Parameterwert z, des Sattelpunktes ent- 
spricht, soll als de günstigste Abbildung der Skala betrachtet werden. Für r- 
Werte nahe Null und für sehr große Werte von r wird nach (7) die Schwankung von T beliebig 
groß. Bei konstant gehaltenem 7, hat 7 bei x, ein Maximum oder- Minimum; bei konstant 
gehaltenem x, hat aber T stets ein Maximum bei z,, d.h. der Teilstrichabstand ist dort größer 
als bei anderen Verzerrungen. 

Ein Fluchtliniennomogramm mit den drei geradlinigen Skalen 


Done made Brei aa 


D2=205: 


b;ı €iı 
v =f;(t)b; +c; mitden Vektoren 6, =|b;, |, Tat PR NL 2) 
0 1 


wird durch die ebene Kollineation 
eur mit A =ila;;), (ze, I U (13), 
transformiert. Durch Parallelprojektion einer geraden Skala auf eine Koordinatenachse erhält 


man eine ähnliche Skala, d.h. aber, daß die Teilstrichkennzahl dabei erhalten bleibt. Die zwei- 
reihige Matrix W für die Projektivität zwischen der Skala r; (12) und der Projektion des kolli- 


nearen Bildes z. B. auf die x-Achse ist 
i _ (1 bi + An bie» AıCiı F ia Cia + us 
A, = a ee nn 
Az, biı + Ay bias Azı Ciı + Agg Cig + ss 
Der zu dieser Projektivität gehörende Verzerrungsparameter 7; ist nach (4): 


eye ER Es UL EN SER 
Azı Ci + gg Cig + Qsg 


oder auch 
le —1) 1 —bilaı + lu — Das —bielae + u Das =0...... (16). 
Gibt man also z.B. nach (9) die günstigsten 7; für zwei Skalen vor, so kann man aus zwei Gln. (16) 
das Verhältnis az} : dag : Ga, berechnen. Der Verzerrungsparameter der dritten Skala ist dadurch 
bestimmt, denn aus drei Gln. (16) folgt: 
(H— Deu — du (4 — 1) 2 — bie n—1 
(1, — 1) Ca — bar (— 1) 09a — ba a EA 
(73 — 1) 05 — b5ı (73 — 1) 35 — by er 
Eine Kollineation zur günstigsten Abbildung eines Nomogrammes mit geradlinigen Leitern kann 


man daher bestimmen durch Original und Bild von drei Punkten und durch zwei optimale Ver- 
zerrungsparameter. 


Verallgemeinerte Multipolfelder 
Von R. Lüst, A. Schlüter und E.-Trefftz (vorgetragen von R. Lüst) 
Max-Planck-Institut für Physik in Göttingen 


Zur Untersuchung von zylindersymmetrischen Vektordifferentialgleichungen, insbesondere 
um siein ein System gewöhnlicher Differentialgleichungen zu überführen, die nur noch vom Mit- 
telpunktsabstand r abhängen, werden skalare und vektorielle verallgemeinerte Multipolfelder 
definiert. Denn häufig ist es bei solchen Problemen sehr unzweckmäßig, ein spezielles Koordi- 
natensystem einzuführen und die Vektordifferentialgleichungen in Komponenten aufzuspalten, 
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5 4 da dabei leicht die Übersichtlichkeit über die Bedeutung und Größe der einzelnen Glieder, so vie 


die geltenden Identitäten — z. B. div rot = 0 oder rot grad — 0 — verlorengeht. Insbesondere = Ei 


r wird durch die Einführung spezieller Koordinaten die exakte oder näherungsweise Überführung 
.. Ta der partiellen Differentialgleichungen in ein endliches System gewöhnlicher Differentialglei- 
Er chungen oft praktisch undurchführbar. Ir TEL EERE PL 

Zu Problemen dieser Art wird man z. B. geführt, wenn man die Rotation kosmischer Gas- 
massen unter dem Einfluß turbulenter Reibung und bei Gegenwart von Magnetfeldern unter- 


, suchen will. a 

e> Jedes beliebige zylindersymmetrische Vektorfeld läßt sich zerlegen in ein Feld DO + ®, das 
eier in den Meridianebenen verläuft, undin ein Feld T, dessen Feldlinien parallel zur Aquatorebene 
5 E; verlaufen. Dabei bedeutet Q den rotationsfreien und ® den divergenzfreien Anteil in der Meri- 
A: dianebene, während T stets divergenzfrei ist. Diese Vektorfelder DO, ® und T werden nun in 
1002», folgender Weise entwickelt. a 

er An Stelle des Potentials der Punktladung 1/r, aus dem sich die gewöhnlichen Multipolfelder 
= - "herleiten lassen, wird eine beliebige kugelsymmetrische ‚„Stammfunktion“ g(r) eingeführt. Es sei 
Br. nun ein skalares Potential nullter Ordnung definiert durch die Gleichung: 


Sit=9(). 


Die höheren skalaren Potentiale der verallgemeinerten Multipolfelder werden analog zu denen der 
= % gewöhnlichen Multipolfelder definiert: 


= Snyılg} = (m grad S,{g}); n=0,1,.... 


Hierin ist m der Einheitsvektor in Richtung der Symmetrieachse. Durch Bildung des Gradienten 
von S„{g} wird dann ein „‚Quellfeld‘‘ Q,{g} defniert: 


Q„fg} = grad 8,{g}; n=(0,1,.... 


Jede Feldlinie dieser Felder liegt ganz in einer Meridianebene und ein solches Feld ist natürlich 
wirbelfrei. Aus diesem Quellfeld wird ein ‚„‚toroidales‘‘ Vektörfeld T,{g} gebildet: 


T.{g} == [On-ı{gt m]; la oc 


Die Feldlinien dieser Felder verlaufen parallel zur Äquatorebene und außerdem ist ein solches 
Feld quellenfrei, wie durch direktes Ausdifferenzieren leicht bestätigt werden kann. Durch Ro- 
tationsbildung gelangt man weiter zu einem „poloidalen‘‘ Vektorfeld: 


Pig = rot Sig) meine: 
das quellenfrei ist und wie Q,{g} ganz in der Meridianebene liegt. 
Für die so eingeführten Felder gelten folgende einfachen Differentiationsregeln: 


div D,fg}= S,{Ag}; rot Bulg}= Tut Ag}; 
AS„{g} = S„{Ag} R A D.fg} — D{Ag} ; 
AT, {g} = In{4g}; AR.{} = Pntdg} . 


Der A-Operator (Aa = grad diva —- rotrot a) kann also bei allen Multipolfeldern direkt auf die 
Stammfunktion abgewälzt werden. 

Für die skalaren Potentiale S,{g} kann man nun Reihenentwicklungen nach (mr) = z be- 
kommen, in denen nur die Ableitungen der Stammfunktionen g(r) auftreten, (die Symmetrie- 
achse sei die z-Achse) und weiterhin lassen sich mit Hilfe dieser Reihenentwicklungen auch die 
Vektorfelder D,, T, und ®, darstellen. Außerdem sind Reihenentwicklungen für die Produkte 
zweier skalarer Potentiale sowie auch für die Produkte zweier Felder möglich. Schließlich kann 
man auch Reihenentwicklungen für eine vorgegebene zylindersymmetrische Funktion fr, 2) bzw. 
für ein vorgegebenes zylindersymmetrisches Vektorfeld angeben!). 

Als Beispiel zur Anwendung der verallgemeinerten Multipolfelder sei auf eine Lösung der 
hydrodynamischen Gleichungen verwiesen, bei der die Stromlinien d überall parallel zu den 
Wirbellinien rot vd verlaufen?). Dieses Problem hat seine Analogie bei der Untersuchung kraft- 
freier Magnetfelder, bei denen das Magnetfeld überall zum Strom parallel ist). In der Hydro- 
dynamik magnetisierter Materie sind gerade solche Magnetfelder von Bedeutung. 


it) SiehehierzuR. Lüst, A. Schlüter und E. Trefftz,J.of Rat. Mech. and Analys. 2,485 (1953 
?2) R. Lüst und A. Schlüter, Z.angew. Math. Mech. (im Druck). IE: h 
®) R. Lüst und A. Schlüter, Z.f. Astrophys: 34, 263 (1954). 
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Zur Interpolationstheorie der reellen Funktionen 
’ gar Von W. Quade in Hannover 

F Die Lag ra ngesche Interpolationsformel und damit auch die Newtonsche weist 
zwei Defekte auf, die den Wert dieser Formel beeinträchtigen. Zum einen strebt der Lagrange- 
sche Interpolationsausdruck bei zunehmender Ordinatenzahl im allgemeinen nicht gegen die zu 


 Interpolierende Funktion, zum anderen bewirkt eine geringe Änderung einer der gegebenen Ordi- 


naten eine unverhältnismäßig große Änderung der zwischen den gegebenen liegenden Ordinaten. 
Es wird ein Interpolationsverfahren vorgeschlagen, das von den beiden soeben erwähnten 
Mängeln frei ist. Das Verfahren benutzt Ausdrücke, die sich additiv aus einem Polynom und 
einem endlichen trigonometrischen Ausdruck zusammensetzen. Es werden Fehlerabschätzungen 
angegeben, die erkennen lassen, in welcher Weise der Interpolationsausdruck gegen die zu inter- 
polierende Funktion strebt. 
Eine ausführliche Darstellung wird demnächst in dieser Zeitschrift erscheinen. 


Die Orthogonalmethode zur Lösung von halbhomogenen | 
partiellen Randwertaufgaben *) 


Von Sven T. A. Ödman in Djursholm, Schweden 


Es sei ein halbhomogenes, lineares Randwertproblem 2 m-ter Ordnung 
L(w) =f(z,y); U;(w) =0; =1,2,--:-4m) 


gegeben. Gesucht ist eine eingliedrige Näherungslösung y = c - u(2; A) -v (y; u), die die Rand- 
bedingungen wenigstens angenähert erfüllt und die drei unbekannte Parameter ce, A und u enthält. 

Mit einem Ansatz y=X (x): Y(y) und unter Vernachlässigung des Restgliedes wird die 
partielle Differentialgleichung mit zwei gewöhnlichen ergänzt und läßt sich dann symbolisch in 
entweder der Form ZL, (X) =f, (x) oder der Form Z,(Y) = f,(y) schreiben. Die Koeffizienten 
sind dabei veränderlich, d. h. abhängig von y bzw. x. 

Um die drei Gleichungen gleichzeitig zu lösen wird X=c, :u(x) und Y= c, : v(y) gesetzt, 
WO Cı'Cg —=c. Gemäß der Fehlerproduktmethode, die eine Variante der Fehlerquadratmethode 
ist, wird verlangt, daß das mittlere Fehlerprodukt der zwei Fehlerfunktionen jeder Gleichung 
möglichst klein wird. Man setzt also 


N [L(euv) — f|(ceuv— w) de dx = min 
u H—Hlau— Nds min 
IE N—h] (6, v — Y)dy = min 
Wenn sämtliche, zugehörige Differentialgleichungen selbstadjungiert sind, dann läßt sich 
zeigen, daß die drei Minimalforderungen in folgender Form umgeschrieben werden können 


SJ [L(euv) —2fleuvdedy = min 
J [Z, (ce, u) — 2 fı] cı u dx —= min 
J IL; (ev) —2 f,] cv. dy = min 
Nach einer Ableitung mit Bezug auf die respektiven Konstanten und nach Wiedereinfüh- 
rung der ursprünglichen Bezeichnungen ergibt sich 


JEW) —fMwuvdady — 0 
]EL@) — fl uda =0 ee en 2 a 1 RE EN 
TEm—irdy  =0 | 

das als 1. Orthogonalformel bezeichnet wird. 

Nach der nahe verwandten Ritz schen Methode wird bekanntlich ein Minimum eines ein- 
zigen Integrales verlangt und dies mit Bezug auf jede ausführbare Variation des Integrales. Hier 
ist dagegen gefordert worden, daß drei Integrale ihren Minimalwert erreichen, und dies für eine 
Variation nur mit Bezug auf die amplitudenbestimmenden Konstanten — dieses Verhältnis führt 
zu einer erheblichen Verminderung der numerischen Rechenarbeit.— Die Methoden können das- 
selbe Resultat ergeben wenn erstens alle drei Differentialgleichungen selbstadjungiert sind, und 


*) Die Methode ist in folgender Arbeit eingehender beschrieben worden und mit Beispielen versehen: 8. T.A. 


Ödman, Studies of BoundaryValue Problems, Part I: Minimum Methods for Solving Semi-Homogeneous Boun- 
dary Value Problems, Svenska Forskningsinstitutet för Gement och Betong, Handlinger Nr. 20, Stockholm 1953. 
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Nomogramme zur Ermittlung der mag 

in stromdurchflossenen dünnen ‚Zylind 

Von Elenore Schwarz in Brin “ 

Es werden Nomogramme beschrieben, die die Radial- und Axialkomponent der m: 

schen Feldstärke in stromdurchflossenen dünnen Zylinderspulen in der gleichen uerm 

gestatten wie das 1952 von A. Weigan d angegebene Verfahren. Dieses er | Geg: nsatz. 
zu älteren Methoden die Bestimmung der Feldstärkekomponenten nicht nur in Punkten der 

Achse sondern in allen Punkten des Feldes. Das Verfahren führt auf die Gleichungen ; : 


on FC) AARON 
: UNO sr le ee een ea Re, 5 
wen MR (1 +2) K(k}) N 
ee Va Ya +? Hr +2): 
er: +3, [Al PD + Aus, PB ee en nn, 
Br er a " Ar Be B | B 2 
= e . ki "AdArnIpl—z: =/SIU TAT TI re ae ee er ee BR IR = 
Se Ar en | 
E kz FÜHRT Sin? 2 Tr ns ae ER SER TR ; (4), “ 
$: ; a) 
: r Te ae, ke Free ee Ad ee Ms 0% 
x een 
k yid+2 
i Be N I ER 
> Ne 


Dabei sind y, r, z die unabhängigen Veränderlichen und 47 bzw. Hf die gesuchten Komponenten 

der magnetischen Feldstärke. C (ky), K (ky) sind vollständige elliptische Integrale zweiter Gat- 

tung und A, (a,, ß/) elliptische Integrale dritter Gattung. : 
Aus den Gleichungen (1) bis (6) werden neue Zusammenhänge abgeleitet, die eine nomo- 

graphische Vertafelung zulassen und auf zwei einfach zu handhabende Nomogramme führen. 
Siehe Deutsche Elektrotechnik, Jhg. 8, Heft 3, 1954, S. 83—87. 


Untersuchung spezieller Funktionen der Potentialtheorie. 
mit Hilfe der konformen Abbildung rn 


Von Friedemann Stallmann in Gießen 


In der Theorie der konformen Abbildungen von Kreisbogenpolygonen wird gezeigt, daß 
der Quotient zweier linear unabhängiger Lösungen einer linearen Differentialgleichung zweiter 
Ordnung mit reellen rationalen Koeffizienten eine Halbebene auf ein Gebiet abbildet, das von 
endlich vielen Kreisbögen und evtl. endlich vielen logarithmischen Windungspunkten berandet 
ist; ein solches Gebiet sei als Kreisbogenpolygon bezeichnet. Ist die Differentialgleichung von 
der Fuchsschen Klasse, hat sie also nur Stellen der Bestimmtheit, so erhalten wir gewöhn- 
liche Kreisbogenpolygone ohne logarithmische Windungspunkte. Kreisbogenpolygone mit 
logarithmischen Windungspunkten — wir wollen sie außergewöhnliche Kreisbogenpolygone 
nennen — gehören zu Differentialgleichungen mit Unbestimmtheitsstellen. Wie der Verfasser 
zeigen konnte !), läßt sich zu einer gegebenen Differentialgleichung das zugehörige Kreisbogen- 
polygon durch asymptotische Integration ohne großen Rechenaufwand näherungsweise bestimmen. 


!) F.Stallmann, Konf. Abb. v. Kreisbogenpolygonen I. Math. Zeitschr. 60, 187—212 (1954). 


en 


und 


> 


heran werden, um die Differentialgleichungen zu untersuch 
" Potentialgleichung Au = 0 bzw. der Wellengleichung Au+ R®u=0 
inaten auftreten. Danach gehört zu jeder derartigen Differentialgleichung 
lygon, und die grundlegenden Eigenschaften der Lösungen drücken sich in 
»sonderheiten dieser Polygone aus. Insbesondere können auch die Randbedingungen 
uftretenden Eigenwertprobleme leicht als geometrische Bedingungen für die zu- 
n Kreisbogenpolygone formuliert werden ®). Die Diskussion der Differentialgleichungen 
Lösungsfunktionen wird damit zurückgeführt auf anschaulich geometrische Betrach- 
tungen und man gelangt zu sehr übersichtlichen Ergebnissen. Auch für numerische Unter- 


‚suchungen sind diese Betrachtungen — in Verbindung mit den oben erwähnten Methoden der 


asymptotischen Integration — recht geeignet. Eingehendere Untersuchungen dieser Art sind 
vom Verfasser durchgeführt für Lam&sche, Mathieusche und Sphäroidfunktionen; ihre 
Ergebnisse sollen demnächst veröffentlicht werden. 3 


| Bestimmung der Eigenwerte einer elliptischen Differentialgleichung 


mit der Randbedingung der Periodizität (Schrödingergleichung 
in Metallen) : 
Von E. Trefftz 


Gesucht sind die möglichen Energieniveaus eines Elektrons im Innern eines Metalls. Sie 
bestimmen sich als Eigenwerte der Differentialgleichung 


— 5 Ay + Ve)y=Ey. 
Dabei ist V(x) das effektive ‚‚Potential‘‘ des Kristalls, in dem sich das betrachtete Elektron 
bewegt. Es ist eine dreifach periodische Funktion des Ortsvektors rt, 
Yke+9)= VW a „Gittervektor“ . 
Für die ‚‚Wellenfunktion‘ y folgt daraus, daß 
y(t-+a)=eite y(r) (Blochsche Bedingung) !) 


Bild 1. Elementarzelle für raumzentriertes Gitter Bild 2. Elementarzelle für flächenzentriertes Bild 3. Elementarsektor im f-Raum 
Gitter (Erste Brillouinzone zuBild1) 


Bild 1 und 2 entnommen aus Handb. d. Phys. 24,2 S. 402 


ist, mit willkürlich wählbarem Wellenvektor f. Die Aufgabe besteht darin, eine Beziehung 


- zwischen f und E herzustellen, E = E(f), durch die dann die Verteilung der Elektronen in der 


Energieskala bestimmt ist. Es genügt y(t) in einer „Elementarzelle“ (Periodizitätsbereich) zu 
betrachten. Diese legt man zweckmäßig so, daß das Atom (es sei nur 1 Atom je Periodizitäts- 
bereich vorhanden) im Mittelpunkt der Zelle (Koordinatenursprung) sitzt. Dann ist für kleine t 
das Potential V(t) kugelsymmetrisch und wird wie 1/|t| unendlich. Die Berandung der Zelle 
ist gegeben durch Flächen, die die Verbindungslinie benachbarter Atome in der Mitte senkrecht 
durchschneiden. Auch in bezug auf fist Y und # periodisch 


E(t+2rb)= Ei) 


2) F. Stallmann, Konf. Abb. gew. Kreisbogenvierecke als Eigenwertproblem. Math. Zeitschr. 59, 


211—230 (1953). 
ı) Sommerfeldt-Bethe, Hab. d. Physik 24,2 p. 368 ff. 


hre 
t= = an) hat auch Punkt H. 


angesehen werden. h 
Superposition von einer vonr— 
einer in r— 0 a a 


Daß beide Anteile zusammen Kr Schrö din; 
friedigen müssen, gibt eine lineare Beziehung z 


Yout = = Yreg - di yieg: 
Eine zweite Beziehung ist dadurch gegeben, daß PEN der reg [ 
Anteil zusammensetzt aus den auslaufenden Wellen 1 
Zentren 


Proglt) — = = ete Yya(t + A) = yo. 


entsprechend der Blo c # schen. Bedingung. Die Kompatibilität “ 
beider Beziehungen liefert den gewünschten Zusammenhang zwi- 

schen Eundf. Es wird also in dieser Methode die Erfüllung der 
Differentialgleichung und die der Randbedingung in gleicher 


S ir un H Weise behandelt. Zur praktischen Durchführung wird y in ein u 
ee: Bild 4. Energiein Abhängigkeit orthogonales Funktionssystem entwickelt. Die Winkelabhängig- 5 
E 2 N see keit ist durch Kugelflächenfunktionen (Ordnung I) gegeben, die ae 
E n y 


d; 


t=(0,0,2) = (y,Y, 2) Yın <®) 
Bild 5 \ 


nach ihrer Symmetrie der kubischen Dreh-Spiegelungsgruppe gegenüber geordnet sind ( eubic 
harmonies‘ %)). s 
Die Radiusabhängigkeit von y ist dann durch halbzahlige Zylinderfunktionen gegeben, 


regulärer Anteil iılku r) = ; Fr 1 0 
-Aykor 2 
auslaufender Anteil h,(k,r) = IE ur A, ı (kor) 
or = 
; sin i 
mit = und he) = = 


?2) Bouckaert-Smoluchowski,-Wj 
; gner, Ps Rev. 50, 58 (1936). 
®) Korringa, Physica 13, 392 (1947); vgl. auch 
gl. auc Härri son,Ph 
) v.d.Lage- Bethe, Phys. Rev. 71, 612 (1947). ae Ba 


—— 


N =. = Ye Y ee j 7a ‚ ka + . h 5 7 ge j 
e x kois en N r die Energie bezogen auf das (annähernd konstante) Potential V,im Randgebiet, a 
m. Damit werden /’und $ Matrizen und die Kompatibilitätsbedingung lautet Bi: 
en Tan Dem a1 00, BE 


00.0 Die Hauptarbeit der praktischen Berechnung liegt in der Ermittlung von T. Im alge- 
Bi meinen hat man bei Entwicklung bis zur 3. Ordnung eine hermitesche Matrix 16. Grades, d.h. x 
136 Koeffizienten, die ihrerseits Funktionen von f und k, sind. Liegt fin einer Symmetrieebene 
R oder gar -achse, so reduziert sich die Zahl der Koeffizienten entsprechend dem Zerfall der Wellen- 
funktionen in Symmetrieklassen. Da /'nur von der Symmetrie des Gitters, nicht vom speziellen 
Ansatz V(t) abhängt, können die Werte für sämtliche Kristalle gleichen Gittertyps verwendet ” 
werden. S hingegen ist wesentlich abhängig von V(t). In erster Näherung kann jedoch V kuge- 
symmetrisch angenommen werden und die Berandung der Elementarzelle durch eine Kugel a. 
ersetzt werden (z. B. gleichen Rauminhalts). Dann ist S = $(E) eine Diagonalmatrix, die nur . 
4 . soviel verschiedene Elemente enthält als /-Werte auftreten, bei Entwicklung bis zur 3. Ordnung Be 
E“ also 4. Diese müssen als Funktionen von Z durch Integration der Schrö dinger gleichung ge- wi: 
funden werden. Bild 4 zeigt die Ergebnisse für Niobium-Metall (Z=41)für den Fal {=(0,0,%5) 
(0 <x <1). Die Buchstaben a, ..., &, &’ bezeichnen die Symmetrieklassen der zugehörigen re 
Wellenfunktionen. y ist zweifach, ö, e, e’ dreifach entartet 2). Für x +0 oder 1 hat man ent- DRG: 
sprechend der geringeren Symmetrie des Wellenvektors weniger Symmetrieklassen (A,,... 4;) 2), 7 A 
die sich als Linearkombinationen der Klassen a, ... e’ darstellen lassen . A, ist zweifach entartet. Be. 
Für benachbarte f-Vektoren würde die Energie aufspalten. Außerdem verschwinden dann die a 
Überkreuzungen: (schematische Zeichnung) B- 


t= (0,0, 2) T= (Yı, Y %) Yo<2). 


Die Energie E=1 entspricht der Energie eines freien Elektrons mit der Wellenzahl |f|=1 
(in unserm, dem Gitter angepaßten Maßsystem). Die endgültige Bestimmung der Elektronen- 
dichte als Funktion von E kann erst vorgenommen werden, wenn im f-Raum soviel Punkte 
berechnet sind, daß dazwischen interpoliert werden kann. e: 

Eine ausführliche Darstellung wird in den Phys. Rev. erscheinen. 
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Über die glatte Ausgleichung von Wertefolgen 
Von Fr. Wecken in Haltingen 


Es seien n reelle oder komplexe Zahlen &,(» = 1, -,n) gegeben als Näherungswerte bei 
t =»für eine Funktion f (f), die für 1<t<n stetige (k— 1)te und stückweise stetige k-te 
Ableitung hat. Erfüllt g(t) diese Differenzierbarkeitsbedingungen, so verwenden wir 


ar a A ch 


n n 
Ag) = 3 |g@) — z,|? und B(g) = [ |g® (t)| dt als Maß für die Anpassung von g(t) an die x, 
1 1 


bzw. für die Glätte von g(t). Wir schreiben {z,,  - ,2,} = {%,}ıs»sn) =t. Bei festem rt, k und 
A > 0 stellen wir das Problem: 


ADAABÜü)EemIn RT RBTRRETE: (1). 


Ohne Beschränkung der Allgemeinheit kann dabei gE®P(t) = 0 fürv<t<v H+1O<r<n) 
gefordert werden. A(g) und B(g) hängen dann nur von den nk Werten gb) O<x<K, 
0<»v<n) bzw. den k Vektoren g® = {g® (v)}o<,<m) ab. Die zu (1) gehörigen Normal- 
gleichungen haben eine eindeutige Lösung 9” = f® von der Form 


ee De en 


die Matrizen M und D, des Grades n hängen von k, nicht aber von / und den x, ab. Aus den f(”) 
ergibt sich f(f) eindeutig durch Polynominterpolation in jedem Intervallv<t <v+1(0<v<n); 
dieses f(£) ist die Lösung von (1). — Bei festem r und k hängen /(t), A(f) und B(f) analytisch 
von A ab. Die Grenzfälle A> 0, A> » schließen sich stetig an und erfüllen die Forderung 
B(g) = min = B,bei A(g) = (reine Interpolation) bzw. A (g) = min = 4, bei B = 0 (Aus- 
gleichspolynom (k — 1)ten Grades). — Das gewonnene Verfahren ist jedoch schwerfällig; insbe- 
sondere stört die Abhängigkeit aller Größen von n. 


Für en AnenillicheN J 
wesentlich handlicher. Ist eine 
gr EANZUEB auf — o<v = ©. — 


n z t - vg Fi en re so 
chungen mit eindeutiger Auflösung (2). Dabei sind 

Q{2,} = {z,}1} erklärten Operators Q mit von k ab IR: 
ae (nach, I: J. Schoenberg » a Sr Math. 4, ” ü 


Istr € a en - € ER T— fE Bnazio, Sun € ee ae ea 
zu fordern. — f(f) entsteht wie oben durch Polynominterpolation. Die Zuordn 
für} = 0 eine gewöhnliche, für A > 0 eine glättende Interpolationsformel nachSchoenb 
Br ist, ebenso wie @ und die Differentiatoren Dx, exakt für Polynome vom Grad 2k — 1. A 

E, nicht allgemein ausführbar. Für die „harmonischen Schwingungen“ we dee ee 
x (—n<u<sa)git Gy, =y (u) Y, mit der oharakteristischen Funktion y (u) = r (et); st 
0<yW<ifüru+0,ywW)=1—Au#+..-füru 0. Bei Anwendung auf Linearkom- 
binationen der 4, wirkt @ wie ein phasentreues Inearas Übertragungsglied mit dem Dämpfungs- Ss 
faktor y (u). Für k < 6 läßt sich auf Grund der expliziten Gestalt von @ noch aussagen: & (u) 
Di - fällt monoton mit wachsendem ||; mit wachsendem % nähert sich (u) einem Verlauf, der einem 
Dr - Tiefpaß entspricht, wobei die Grenzfrequenz mit wachsendem A abnimmt. 


BE Der Wert von %k ist in der Regel aus der physikalischen Situation eindeutig bestimmbar. 
“2 Fr Bezüglich A jedoch scheint eine willkürfreie allgemeine Festlegung kaum möglich. Für spezielle 
Ben Fälle legt die folgende Überlegung eine bestimmte, nur auf £ und k bezugnehmende Wahl von A 
- Be, nahe. — Sind x sowie g® (0 <x< k) Linearkombinationen der 4, so lassen sich A (g) und x 
2 B (g) durch Mittelbildungen (wie bei fastperiodischen Funktionen) erklären, und f (f) ist wieder Be 
e-; die Lösung von (1). Werden A(f) und B(f) als Funktionen von A bei festem x und k betrach- ’ 
N. tet, und ist s die Anzahl verschiedener Frequenzen |u| +0, die in der Darstellung von £ vor- 
> kommen, so gilt: für s=0 ist A = B()=0 für jedes 4. Für s=1 ist YJA(N/A- en 
2 + yYB(f)/B,— 1 für jedes 4. Für s> 1 hat YA()JAx + VB(f)/B,—= min genau eine Lösung 
a = 4*. Für s = 2 ist A* nur von den beiden Frequenzen, nicht von den Amplituden abhängig, 
’ . und @ (mit A = A*) bewirkt bei nicht zu kleinem k nahezu eine Reproduktion der langsameren 

und Vernichtung der schnelleren Schwingung. Für s> 2 ist A* auch von den Amplituden ab- 

hängig, jedoch gilt Ähnliches wie bei s— 2, wenn die Frequenzen sich auf zwei schmale, deut- 

lich getrenn te Bereiche verteilen. 
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Referate der Vorträge: 


R. Albrecht, München: Ein Iterationsverfahren zur ee Abbildung 
zweifach zusammenhängender Gebiete auf einen Kreisring; 


E.Bukovics, Wien: Zur Abschätzung des Ba beim Er ui zZ 
und beim Verfahren von Blaess. 


L. Holzer, Rostock: Laplace-Transformation der Besselschen Funktionen. 


G. Lyra, Göttingen: Über eichinvariante nun der Räume konstanter 
Krümmung. 


L. Ochini, Sao Paolo: Determinazione degle zeri di un polinomio a coeffizienti 
reali. 
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 Programmorganisation | bei kleinen Rechenautomaten mit 
R innerem Programm 


N Ra Von H.Pösch und Th. Fromme i in Weil/Rhein. 


man die Programmierung auf einem ‘modernen Rechenautomaten studieren will, 
braucht man ein sprachliches Mittel, um die Arbeitsweise der Maschine darzustellen. 


Schon beim Bau einer Rechenmaschine müssen sich Mathematiker und Elektroniker 


darüber verständigen, was die Maschine leisten soll. Verständigen heißt in diesem Fall die 
gleichen Bezeichnungen und den gleichen Code für. die einzelnen Funktionen der Maschine 
benutzen. Wir haben uns zunächst nur zum Studium der Programmierung und zum Unterricht 
in Programmierungsfragen das Modell einer Maschine überlegt, deren Befehlscode so einfach ist, 
daß man ihn schnell erlernen kann, die aber alles das leistet, was an nn Operationen 
in modernen Rechenautomaten auftritt. 
Zur Ausführung der arithmetischen Operationen muß eine Maschine mindestens folgende 

Bauteile enthalten: 

T..einen Speicher 
 2.ein Addierwerk für die Addition mit Resultatsregister (Akkumulator) 

3. ein Komplementwerk, das zu jeder Zahl die Komplementzahl bildet und 

4. ein Register, das eine Stellenverschiebung der Zahlen vornimmt. (Shiftregister). 
Diese arithmetischen Operationen und Zahlentransporte müssen von einem Leitwerk nach An- 
weisung des jeweils im Befehlsregister b gespeicherten Befehls gesteuert werden. Die 
Befehle für die Zahlenoperationen drücken wir durch ein Operationszeichen und eine Speicher- 
adresse aus. Sie bedeuten: 


Am: (Akku. + <m) = Akku 
Sm: <Akku) — <m) => Akku 
im. (Akku) — m 


<) bedeutet: „Inhalt von“, m heißt Speichernummer und —> „Transport nach“. 
Das Befehlsregister kann seine Befehlstexte von einem Lochband, einer Lochkarte etc. erhalten, 
Wir sprechen in diesem Fall von einer Maschine mit äußerem Programm. 

Die volle Flexibilität der großen Universalmaschinen erreicht man erst mit gespeichertem 
Programm. Die Modellmaschine soll auch mit innerem Programm arbeiten, muß daher einen 
Befehl enthalten, der den Transport vom Speicher zum Befehlsregister besorgt. Wir vercoden 
ihn mit 

Em: <m>—»br 


Der E-Befehl transportiert nicht nur nach b, sondern löst nach Rückführung des arithmetischen 
Befehls den Befehl E(m + 1) aus. Eine solche Maschine könnte zwar alle Befehle (auch die 
bedingten) durchführen, wollte man sie aber realisieren, so ergäben sich zu lange Operationszeiten. 

Da man in der Praxis doch noch einige Bauteile hinzufügen muß, wie zum Beispiel ein 
Bedingungs-Flip-Flop, haben wir sogleich bedingte Befehle vorgesehen. 

Unser Maschinencode weicht also von dem einer theoretischen Minimalmaschine bereits ab. 
Für die Bedingung verwenden wir die Zeichen P und Q. Dabei bedeutet ein P vor irgend einem 
Befehl (z. B. PAm): der Befehl ist nur durchzuführen, wenn der Inhalt eines besonderen Registers 
des Kennregisters, positiv ist. Q reagiert entsprechend auf negatives Kennregister. 

In der Sprache dieser so erweiterten Minimalmaschine wollen wir zunächst einige Programm- 
typen besprechen. 

Als Rufprogramme bezeichnen wir solche Befehlsfolgen, in denen die Maschine veranlaßt 
wird, nach Ausführung der Operationsbefehle wieder an die Stelle zurückzukehren, von der aus 
das Programm gerufen wurde. Schreiben wir vor einen senkrechten Strich die Nummern der 
Speicher und daneben jeweils den Speicherinhalt, so erhalten wir das folgende Beispiel für ein 
Rufprogramm. 


Hauptprogramm Rufprogramm 
400 A 402 er 2 2274300 T320 
401 sr 301 


402 E 403 302 


en ee 
e za, . ... 
az (E 403) 


wu 


in. 
B. Rechenmaschinen 


grammieren nicht mehr zu beachten. Alle Speicherinhalte werden durch ein Lochband einge- 
geben und in jeder Maschine dient ein spezielles Programm, das wir Leseprogramm nennen, dazu, 
die vercodeten Lochgruppen in Maschinentexte zu übersetzen. Dieses Leseprogramm kann man 
benutzen, um eine Vielzahl von Aufgaben zu erfüllen. Im ersten Beispiel wählen wir die soge- 
nannten symbolischen Befehle. Der Mathematiker schreibt auf der Schreibmaschine z. B. x A600, 


sin A 700 mit der Bedeutung: multipliziere den Inhalt von 600 (z. B. y) mit dem letzten Rechen- 


resultat, berechne den Sinus von x, wobei x in 700 steht. Das Leseprogramm übersetzt solche 
symbolischen Befehle durch geeignete Transportbefehle und setzt statt des Mal-Zeichens, bzw. 
statt des Wortes „sin‘‘ einen Sprungbefehl in das Multiplikations-Rufprogramm (zZ. B. E 300), 
bzw. in das Sinus-Rufprogramm (z. B. E 780). Man kann im Leseprogramm, das natürlich die 
Befehle nacheinander in die fortlaufenden Speichernummern n, n + 1: : - bringt, auch Umrech- 
nungen der Adressen der Befehle durchführen lassen. Auch die sofortige Ausführung gewisser 
Befehle kann das Leseprogramm veranlassen. Solche Befehle heißen in der EDSAC Controlcom- 
bination. Zu diesen gehört auch die Notierung der ersten Einschreibadresse n in das Lesepro- 
gramm. Zahlen muß das Leseprogramm in das Dualystem übersetzen unter Berücksichtigung 
des fließenden Kommas. 

Um so vielfältige Aufgaben zu lösen, bedarf es einer besonderen Organisation des Lesepro- 
gramms. Sehr starke Hilfsmittel sind Stellwerk und Weichen. Unter einem Stellwerk verstehen 
wir einen Programmteil, der auf Grund gewisser Informationen (in diesem Fall Lochkombina- 
tionen) einen oder mehrere Befehle, insbesondere einen Sprungbefehl eines anderen Programm- 
teils errechnet und dort einschreibt. Einen solchen variablen Befehlstext nennen wir eine Weiche. 

Diese Beispiele mögen genügen, um zu zeigen, wie bereits an diesem einfachen Modell kom- 
plizierte Programmierungsfragen diskutiert werden können. Darüber hinaus kann dies theore- 
tische Modell Auskunft darüber geben, welche Rufprogramme häufig gebraucht werden, viel Zeit 
erfordern und daher in einer größeren Maschine durch einen elektronischen Automaten ersetzt 
werden sollten. Wenn man viele Programme auf der einfachen Maschine programmiert hat, kann 
man sogar die Vorteile einzelner elektronischer Automaten gegeneinander abwägen. Auch kann 
man die Frage diskutieren, wie weit es möglich ist, den einzelnen Rufprogrammen ein spezielles 
Unterbefehlszentrum zu geben und sie in Mikroprogrammierung schnell durchführen zu lassen. 

Man kann auch daran denken, eine solche einfache Maschine zu bauen. Natürlich muß man 
der Rechengeschwindigkeit halber noch einige Elektronik hinzufügen. Eine solche Maschine 
wäre für Institute, die über geringe Geldmittel verfügen und schon mit einer kleinen Rechenge- 
schwindigkeit auskommen, interessant. Es ist auch in Erwägung zu ziehen, ob neben einer großen 
Maschine eine kleine folgenden Zwecken dienen kann: 


1. Der Übersetzung langer und komplizierter Programmfolgen aus der dem Mathematiker ge- 
läufigen Formelsprache (Superprogramm). 


2. Zur Erprobung der für die große Maschine aufgestellten Programme. 


Man kann der kleinen Maschine ein spezielles Ausführungsprogramm geben, das alle 
Schritte der großen Maschine nachahmt. Die kleine Maschine würde dann gestatten, Pro- 
grammkontrollen ohne Unterbrechung der Arbeiten der großen Rechenmaschine durchzu- 
führen und insbesondere dann Zeit und Geld sparen, wenn auf der großen Maschine viele ver- 
schiedene Programme von vielen Mathematikern mit unterschiedlicher Maschinenerfahrung 
durchgeführt werden sollen. 


Der dreiwertige Aussagenkalkül der theoretischen Logik und seine 
Anwendung zur Beschreibung von Schaltungen, die aus Elementen 
mit zwei stabilen Zuständen bestehen 


Von Hans Rohleder in Dresden 


Zur Beschreibung von elektrischen Schaltungen wird normalerweise der Aussagenkalkül 
der theoretischen Logik verwendet. Dabei ist es üblich, neben den logischen Grundverknüpfungen 
wie Negation, Disjunktion und Konjunktion, das Gleichheitszeichen zu verwenden. Ist z.B. 
ein Relais X über den Kontakt A mit der Spannungsquelle verbunden (Bild 1a), so beschreibt 
man den Sachverhalt durch die Gleichung 


XeA... ee} 


Als Rufprogramme benötigen wir zunächst die Operation in fließendem Komma (4, x > 
aber auch komplizierte Operationen (e*, sinx) können so organisiert sein. Sind die Rufpro- 
gramme in der Maschine fest notiert, so braucht der Mathematiker ihre Organisation beim Pro- 
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5 dürch‘ Msiheh, Wert von “X erfüllt werden. Diese Schwierigkeit kann Pe werden, wenn. 
man den zweiwertigen Aussagenkalkül zu einem dreiwertigen verallgemeinert. 


Die im folgenden betrachteten Aussagen sollen immer drei Wahrheitswerte haben en 
Der Wahrheitswert einer Aussage A, dermit W(A) bezeichnet werden soll, ist eine Zahl und zwar 
entweder 0 oder 1/2 oder 1. Zu den aussen Werten O0 und 1 kommt also als dritter ve 2: 
(unbestimmt) hinzu. 

Das Relais X in Bild 1b zieht an ve fällt gleich. ‚darauf wieder en usw. Schaltet man 
einen Kontakt von X in Reihe mit einem geschlossenen Kontakt eines anderen Relais, so erhält 
man die gleiche Ausgangsspannung, als wäre nur der Kontakt von X vorhanden. Berücksichtigt 
man alle entsprechenden Ergebnisse und verlangt, daß die Konjunktion auch hier die Reihen- 
schaltung beschreiben soll, so erhält man: 


W(A& B) =Min. wa, a ze, 


Entsprechend 
;WlAN-B) = Max: [WA W (De re ee) 
und : 
W(4)=1—W(A) un A: 5), 


Der Kalkül ist mit diesen Operationen noch nicht funktionell vollständig. Deshalb werden 
zwei weitere Operationen eingeführt: 


„ [WM für W(4) #1/2 | 
W (4) -| 0 für W(4)=1/2 ’ ©) 
WA Tue WA) 12 


Eine dieser beiden Operationen ist an und für sich entbehrlich; sie stehen in einem ähnlichen 
Verhältnis zueinander wie Konjunktion und Disjunktion. 


Bildi Bild 2 Bild 3 


Die Einführung der beiden zuletzt genannten Operationen in dieser Form wird durch den 
Erfolg gerechtfertigt. Der zweiwertige Aussagenkalkül, wie er z.B. im ersten Kapitel von 
D. Hilbert, W. Ackermann „Grundzüge der theoretischen Logik“ ?) dargestellt wird, 
kann nämlich bis auf den Teil, der sich mit der Axiomatisierung beschäftigt, fast vollständig 
übertragen werden. 

Die Analogien gehen so weit, daß selbst das Dualitätsprinzip beinahe wörtlich übertragen 
werden kann. Es lautet: Von einer in Normalform vorliegenden Funktion erhält man die Nega- 
tion, indem man an allen Stellen die Konjunktion und Disjunktion sowie die Operationen ° und ” 
vertauscht und alle Grundaussagen durch ihre Negation ersetzt. 

Die Angabe näherer Einzelheiten über den dreiwertigen Kalkül sowie auch die Erklärung 
des im Dualitätsprinzip verwendeten Begriffs Normalform würde hier zu weit führen. 


1) Das Konjunktionszeichen & wird im folgenden fortgelassen. 
2) Springer-Verlag: Berlin 1949. 
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- zeichneten Wahrheitswert nichts zu bemerken. 


‘werden. Es scheint deshalb wünschenswert, auch bei der 
_ zeichneten Wahrheitswert auszukommen. 


wird®). Bei allen Betrachtungen, die fü 


vollkommen gleichberechtigt und zueinander dual. 
hungen unabhängig davon, ob ein oder zwei ausge 


Im folgenden soll die Verwendung des dreiwertigen Kalküls zur x Beschrei 


tungen betrachtet werden. 
Jede Schaltung wird durch ein GSChInae SER beschrieben. = Big 2 2 


System 


r ER 
Z=AyBYa u DQsDXF RR ae et. 


wenn X wird angezogen, wenn entweder A geschlossen ist oder B und Y, und Y wird angezogen, 
denn ( geschlossen ist oder X und D. | 

In Bild 3 sind dieselben Symbole wie in „Synthesis of Electronic Computing and Control 
Cireuits‘“) verwendet. Das Potential an einigen Punkten, bei denen der Potentialverlauf 
nicht sofort klar ist, wird als Unbekannte des Systems eingeführt. Man erhält so das Gleichungs- 
system: 


x -AvB(OsD)Y(OvD) 
Y=dvP(BvDvX{ayD)| no Se ze 


ZIXHBYIR 
Allgemein erhält man ein System von der Gestalt: 
KT (4,B:. 0,2, 142 
4 RE EP =, FEN 
Zeh (ABER CNN 82) 
Dabei wirken die „Elemente“ A, B,...C von außen in die Schaltung hinein, die „Elemente“ 
X,Y,...Z werden erst durch das System definiert. 

Als nächstes werden die Lösungsverhältnisse eines Gleichungssystems der Form (10) be- 
trachtet. Es kann gezeigt werden, daß jedes System, welches eine Schaltung beschreibt, minde- 
stens eine Lösung hat. In der Bedingung ‚‚welches eine Schaltung beschreibt‘ liegt eine Ein- 
schränkung der Funktionen f, 9, .... h. Diese Einschränkung muß natürlich exakt formuliert 
werden, indem die Klasse der zulässigen Funktionen genau angegeben wird. Das würde jedoch zu 
weit führen. 

Im allgemeinen haben die betrachteten Systeme nicht nur eine sondern mehrere verschie- 
dene Lösungen. Eine dieser Lösungen ist besonders ausgezeichnet. Sie soll Hauptlösung genannt 
werden. 

Die Hauptlösung ist die Lösung, welche am häufigsten den Wert 1/2 für die abhängigen 
Variablen liefert. Faßt man eine spezielle der abhängigen Variablen ins Auge, so wird für ge- 
wisse Wertekombinationen der unabhängigen Variablen vom Gleichungssystem für die be- 
trachtete abhängige Variable der Wert O0 oder 1 vorgeschrieben werden. Die Hauptlösung liefert 
für alle Wertekombinationen, für die nicht unbedingt der Wert 0 oder 1 vorgeschrieben ist, den 
Wert 1/2. Sie.ergibt also auch dann den Wert 1/2, wenn dasbetrachtete Relais Haltespannung hat. 

Die Hauptlösung kann in jedem speziellen Fall rechnerisch ermittelt werden. Das Verfahren 
soll nicht in allen Einzelheiten angegeben werden. Es besteht im wesentlichen aus einer Iteration 
des Systems. Der Iterationsprozeß ist immer nach endlich vielen Schritten beendet. Diese Tat- 
sache hat ihre Ursache darin, daß es bei einer festen Zahl von Variablen nur endlich viele Funk- 
tionen von diesen Variablen gibt. 

Zum Abschluß noch einige Anwendungsmöglichkeiten. 

Mit der Ermittlung der Hauptlösung ist ein Fortschritt bei der Analyse von Schaltungen 
erzielt. Im Beispiel von Bild 2 würde man normalerweise bei den Gleichungen (8) stehenbleiben. 
Durch Einsetzen der zweiten Gleichung in die erste erhält man: 


X = AvBOVBDRT. oe 
Hier ist auch das indirekte Anziehen von X über BC zu erkennen. 


®)Man vgl.z.B.: J.B.Rosser,A.R.Turquette ,,‚Many-valued Logics‘ North-Holland Publishing 
Company, Amsterdam 1952. 


4) Harvard University Press; Cambridge, Massachusetts 1951. 
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> Dieser era ist lach praktisch nicht immer ausreichend; denn bei vielen BI; 
EN Schaltungen erhält man für alle abhängigen Variablen in jedem Fallden Wert1/2. SE ke 
> . Bei einer noch tiefer gehenden Analyse dürfte man nicht nur den Zustand einer Schokuiie 
= bei einer bestimmten Wertekombination der unabhängigen Variablen ins Auge fassen, sondern 
man müßte auch das Verhalten der Schaltung beim Übergang von einer Wertekombination zur Ber 
A anderen betrachten. Dazu ist es jedoch notwendig, bestimmte Annahmen über die verwendeten Fed 
R Schaltelemente zu machen. Bei Relais muß z.B. die Anzugs- und Abfallverzögerung bekannt Bar 
s sein. Esist klar, daß der Aussagenkalkül dazu nicht ausreichend ist. ir 
x - Außerdem wird durch den dreiwertigen Kalkül die Aufmerksamkeit auf eine Klasse von Be 
j Röhrenschaltungen gelenkt, die bisher nicht beachtet worden ist. Von der Schaltung Bild 3 wird, 
wenndAs)K Bere Aussagen sind, der logische Ausdruck =B 


—[AvAvB(OvD)vGRODI[GVEVF(BVDWABED ....... 0 


gebildet. Die Pr hat also mit einem Flip-Flop nichts gemeinsam. Der Ausdruck könnte 
auch ohne Rückkopplung gebildet werden. Es wären dann aber mehr Röhren erforderlich. a 

Schließlich kann eine rein rechnerische Methode zur Aufstellung von Kontaktbrücken- "7 
schaltungen angegeben werden. Die bisher bekannten Methoden?) greifen neben der Rechnung = 
auf Schaltskizzen zurück. 


E Projekt einer elektronischen Rechenmaschine an der E.T.H. Bi: 


(ERMETH) Fi 

E Von A. P. Speiser in Zürich (Inst. f. angew. Math. d. E.T.H. Zürich) | 2 
4 Allgemeines E 

Die Darstellung der 16stelligen Dezimalzahlen erfolgt in Tetraden nach Aiken, indem E b% 

j eder Dezimalziffer vier Dualziffern mit den Gewichten 2*, 4, 2, 1 zugeordnet sind. Ein Wort hat BR; 

4 also 64 Dualziffern. Die Tetraden werden parallel geführt, die Dezimalziffern in Serie. Die EN 

ı Folgefrequenz der Impulse beträgt (außer im Speicherwerk) 32 kHz. Somit benötigt die Über- ee 

- tragung eines Wortes 4 Leitungen und dauert 0,5 ms. N = 
E- Wir unterscheiden 3 Klassen von Wörtern, nämlich: Tr 

1. Zahlen in halblogarithmischer Form a: 10%. Dabei ist der Exponent a durch eine 3stellige Be 

Dezimalzahl dargestellt. Die Mantisse a ist ein Dezimalbruch mit 11 wesentlichen Stellen und 5 ä 


einem Vorzeichen. 

2. Zahlen mit festem Komma. Wir haben für -Zahlen dieser Art infolge des Wegfallens des Expo- 
nenten 14 Dezimalen zur Verfügung. 

3. Befehlspaare. Ein Befehl umfaßt 7 Dezimalen, nämlich 2 Operationsziffern, eine Indexziffer !) 


und 4 Adressenziffern (der Speicher enthält 10 000 Wörter). 


Außerdem ist jedem Wort noch eine 16te Dezimale zugefügt, welche das sog. Q-Zeichen!) 
enthält und im übrigen als Kontrolle dient. 


Befehlsgebung 

Die beschriebene Maschine ist eine Einadreßmaschine, d.h. jeder Befehl enthält eine Adresse 
und eine Operation. Die Adresse bezeichnet den zur Operation gehörigen zweiten Operanden, 
während das Resultat der vorangehenden Operation als 1. Operand dient. Bei gewissen Befehlen 
(Operationen mit nur einem Operanden) fällt der 2. Operand weg, die Adresse ist in einem solchen 
Fall irrelevant. 

Es ist zu beachten, daß jeder Befehl außer der Adresse (die von 0 bis 9999 laufen kann) und 
den Operationsziffern noch eine sog. Indexziffer enthält. Wenn diese Indexziffer 0 ist, so ge- 
schieht nichts Besonderes. Ist sie aber eine der Ziffern 1 bis 9, so wird die Adresse dieses Befehls 
vor seiner Ausführung noch um die im Indexregister ] stehende 4stellige Zahl erhöht. 


Rechenwerk und Leitwerk: Schaltungen 


Der Wert 0 wird durch eine Spannung von —29 V, der Wert 1 durch —7 V dargestellt. Die 
Verknüpfung von Variablen geschieht ausschließlich mittels Dioden-Netzwerken, 
deren Eingänge durch Kathodenverstärker gespeist sind. Zur Bildung der Negation und als 
Zwischenverstärker dienen Inverter (zweistufige Triodenverstärker), die das Signal auf die 


5) H. Piesch, „Über die Vereinfachung von allgemeinen Schaltungen“. Archiv für Elektrotechnik, 
XXXIII (1939), S. 733—746 und M. A. Gawrilow, „Relaisschalttechnik für Stark- und Schwachstrom-An- 
lagen “ (aus dem Russischen), Berlin 1953. 

1) Diese Begriffe sind in „Automatische Rechenplanfertigung‘“ (Mitteilungen aus dem Institut für ange- 
wandte Mathematik Nr. 3) von H. Rutishauser erklärt. 


 Normalpegel verstärken und von Störungen befreien; sie geben sowohl die 


_ gegengekoppelt, so daß ein Röhrenabbrand von über 50% ihre Funktion nicht beei 


auch deren Negation ab und werden durch Kathodenverstärker gefolgt. Die Inverte, 


Kathodenverstärker, Diodennetzwerk und Inverter sind durchwegs gleichstro: ge 


und bilden zusammen eine „statische Schaltung“. An ihren Eingängen stehen Flip-Flops, die BT 
den Anoden die normierten Pegel von —29 V resp. —7 V abgeben. Die Flip-Flops sind, wie die 


Inverter, zur Erhöhung der Betriebssicherheit gegengekoppelt. An den Ausgängen der statischen 
Schaltung stehen ebenfalls Flip-Flops, welche unter Zwischenschaltung von Uhrimpuls-Toren ge- 
schaltet werden. Die schaltenden Uhrimpulse haben eine Dauer von 1 us und eine Wiederholungs- 


_ frequenz von 32 kHz. Vor ihrem Eintreffen müssen alle Ausgleichsvorgänge beendet sein. 


Verzögerungen’erfolgen durch Ketten von Flip-Flops mit zwischengeschalteten Uhrimpuls- 
Toren. 
Verwendete Röhren und Dioden: Ausschließlich E 90 CC und OA 55. 


= 
Reehenwerk und Leitwerk: Organisation | 

Das Rechenwerk besitzt 3 regenerierende Register, welche Sektoren auf der Magnettrommei 
sind; ihre Kapazität beträgt je 1 Wort. In ihnen werden, mit Hilfe des Addierwerkes, die Rechen- 
operationen ausgeführt. Bei Operationen mit gleitendem Komma werden Exponent und Mantisse 
in den gleichen Registern und im gleichen Addierwerk verarbeitet; ihre Trennung geschieht ledig- 
lich durch Zeitstreuung. 

Die Rechenzeiten betragen (ohne Suchzeiten): 


Addition 5ms 
Multiplikation 16 ms 
Division 30 ms. 


Die mittlere Suchzeit für ein Wort (Zahl oder Doppelbefehl) dauert 5 ms, und die mittlere 
Zeit pro Operation in einem typischen Problem beträgt, einschließlich aller Suchzeiten, 16 ms. 


Das Magnettrommel-Speicherwerk 

Das Speicherwerk besteht aus einer Magnettrommel von 25 cm Durchmesser und 50 cm 
Länge; sie rotiert mit 6000 U/min. Die Aufzeichnung erfolgt ebenfalls tetradisch, so daß 4 Magnet- 
köpfe einen dezimalen Kanal bilden. Das Speicherwerk enthält 50 Kanäle (200 Köpfe); in jedem 
Kanal befinden sich 200 Wörter. Die Kapazität erreicht also 10 000 Wörter. Die axiale Breite 
eines Kanals beträgt 2 mm, die Impulsdichte etwa 4 pro mm und die Impulsfrequenz 320 kHz. 
Um beim Ablesen und Speichern eine Anpassung an die 32 kHz des Rechenwerkes zu erreichen, 
werden die Impulse so geschachtelt, daß jeder zehnte Impuls zum gleichen Wort gehört. 

Geschrieben wird im „return-to-zero‘‘-Verfahren, d.h. für die Ziffer O0 erhält der Schreib- 
kopf einen negativen, für 1 einen positiven Stromimpuls von 1 us Dauer. 


Eingang und Ausgang 

Als Zwischenträger für Eingangs- und Ausgangswerte dienen Lochkarten. Die ver- 
wendeten Karten sollen wahlweise auf der ERMETH und auf Lochkarten-Normalmaschinen ver- 
arbeitet werden; daher wird auf den Karten die Darstellung der Dezimalen mit 6 Löchern ver- 
wendet. Auf einer Karte werden 5 Wörter gelocht. 

Die ERMETH hat zwei Kartenstationen; jede kann nach Wahl als Locher oder als Abtaster 
verwendet werden. Die Stationen verarbeiten pro Sekunde 2 Karten. Total können also pro 
Sekunde 20 Wörter gelocht oder abgetastet werden. Die Umschaltungen zwischen den Feldern 
einer Karte und den Stationen erfolgen mit elektromechanischen Koordinatenwählern. Die durch 
die Maschine zum Lochen abgegebenen Wörter werden zunächst in Gruppen von fünf durch Koor- 
dinatenwähler gespeichert, worauf die Karte gesamthaft gelocht wird. 

Endgültige Resultate können auch einer Schreibmaschine odereinem Tabula- 
t or zugeführt werden. Die Schreibmaschine druckt in zwei Sekunden ein Wort; auf einem Steck- 
brett können 20 verschiedene Druckschemata vorgewählt und alsdann beliebig programmiert 
werden. Der Tabulator druckt 2 Wörter pro Sekunde. 


Der Kartenlocher 
Die Vorbereitung und die Kontrolle der Lochkarten durch das Personal geschieht auf dem 
Kartenlocher, der von der Maschine vollständig unabHängig ist und der zwei eigene Karten- 
stationen sowie einen eigenen Tabulator besitzt. 
Weitere Beschreibungen der ERMETH werden erscheinen in-der Zeitschr. f. angew. Math. 
u. Phys. sowie als Mitteilungen aus dem Institut für angewandte Mathematik. 
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En Obere und untere Schranken für den Drillwiderstand 
2 Von E. R. Berger i in Wien. 


nn 


aufgefaßt werden; behandelt man diese nach den direkten Methoden, erhält man obere und 


untere Schranken. Ausgehend vom Spannungshügel, bzw. der Torsionsfunktion T, erhält man 


die untere Schranke: 


Jen „aT—o)dr =4V—G, [Tar=V, gedt=v, [war=G. 


Die een T muß am Kubentand Null kn im en stetig sein. Integrationsbereich 5 


ist der ganze Bereich innerhalb der äußeren Kontur, über Durchbrüchen 9’ muß 7 = const sein. 
Nach der Methode von Trefftz kann das Problem so umgeformt werden, daß man die 


obere Schranke erhält: 


1=jv dF, Das grad. 
Die ech funklion T muß innerhalb der äußeren Kontur, also auch über Durchbrüchen, 
der DGl: AT =—2 genügen. Integrationsbereich ist nur der Querschnitt selbst, ohne die 
Durchbrüche. Wirbelbehaftete Felder d wurden hier willkürlich vom Vergleich ausgeschlossen, 
obwohl ihre Verwendung manchmal raschere Konvergenz liefert }). 

Zur Berechnung dieser oberen Schranke verwendet man zweckmäßig die komplexe Koor- 
dinatez= x +iy. Die Vergleichsfunktion 7 wird angesetzt in der Form 


ie ze 
T=-R|-% +10). 
Dem Vektor v entspricht die komplexe Zahl 


v=—-2z+f@ E99 


Ansatz von f(z) als Taylorreihe ergibt für das Funktional einen quadratischen Ausdruck in den 
Koeffizienten c„; die Extremalforderung liefert ein unendliches lineares Gleichungssystem mit 
Hermite scher Matrix, dessen Koeffizienten die Form haben 


J 2 zi dF = Mys+t, s—t° 


Diese Ausdrücke können als ‚komplexe Momente‘ des Bereiches gedeutet werden. Ihre Ordnung 
ist s + 2 ; sie lassen sich durch die gewöhnlichen Momente derselben Ordnung linear darstellen. 
Das Gleichungssystem wird orthogonalisiert und man erhält den Extremalwert als unendliche 
Reihe, deren Teilsummen eine monoton fallende Folge von oberen Schranken bilden. Die erste 
und die zweite Näherung ergeben die bekannten Schranken °) 


J=J, PLN, 


(Jy = Bes Trägheitsmoment um den Schwerpunkt; J}, Js, = Hauptträgheitsmomente.) 
Zur Bestimmung der unteren Schranke kann man den Spannungshügel zunächst der Höhe 
nach affin verzerren. Das ergibt allgemein den einfacheren Ausdruck !) 


J, >A Ve. 


Für sternförmige Bereiche kann man den Spannungshügel als Kegelfläche ansetzen, besser als 
Kegelstumpf, am besten ist eine Parabel als Meridian, mit dem Scheitel in der Kegelachse. Diese 
Fläche sei als ‚„„Parabelkonus‘“ bezeichnet. Man erhält die Schranke 


ds 
nor dt, 
z p 


wobei p der Abstand der Kurventangente vom (willkürlich gewählten) Ursprung ist. 

Beim Dreieck und Viereck kann man T in bekannter Weise als Produkt ansetzen ?). 
Streifen zwischen parallelen oder leicht geneigten Geraden überdacht man unter Erfüllung der 
DGl mit einem parabolischen Zylinder bzw. einem hyperbolischen Paraboloid ®). Am Ende muß 


1)0. Weber, Z.angew. Math. Mech. 11 (1931), S. 244—45. 
2)W.Hoffer berth, Stahlbau 17 (1944), S. 12—16. 

®)E.Ni coli, Z.angew. Math. Mech. 4 (1924), S. 181—82. 
4)6.Weber,Z.angew. Math. Mech. 2 (1922), S. 299—302. 
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# » Dollieed J, kann bekanntlich 1) 2) als Extremalwert zweier nn. 


_ richtig wiedergeben. Die Schranke Jıy einerseits, u 


‘schied ein. (Jr beim langgestreckten gleichschenkligen Di 


ein geeigneter Abschluß angesetzt werden: 
konus. a Be 

Anwendung auf einige einfache Figuren (R 
liges Dreieck) zeigt, daß die obere Schranke J,, un 


abgeschlossen mit dem Produktansatz, schließen den 


schnitte werden besser angenähert als langgestreckte. t Se 


a Er n j = > G 4 wer 

Ein neues graphisches Verfahren zur Schwungradermittlung 
bei Kolbenmaschinen _ u 

Von K. Desoyer in Wien !) 


Für die Ermittlung des Schwungrades für ein System Kolbenkraftmaschine-Arbeits- 


maschine, bei dessen stationärem Lauf ein gegebener Ungleichförmigkeitsgrad nicht überschritten 
werden darf, wurde ein neues graphisches Verfahren gesucht. Es sollte möglichst geringen 
Zeichenaufwand erfordern, jedoch keine wesentlichen Näherungsannahmen verwenden, die die 
Zuverlässigkeit der Resultate beeinträchtigen könnten. 

Eine Kolbenkraftmaschine sei mit einer beliebigen Arbeitsmaschine zwangläufig gekoppelt. 
Wir untersuchen den stationären Lauf dieses Systems. Die gesamte kinetische Energie T läßt 
sich in der Form anschreiben I: 


TB) 
a En) EA ee 0) 


(9, bzw. d,,; Drehwinkel derMotorwelle bzw. der Welle der Arbeitsmaschine, E, bzw. E,, Getriebe- 
funktionen von Motor bzw. Arbeitsmaschine). Statt der beiden Drehwinkel führen wir eine 
Winkelvariable p so ein, daß sie während der vollen Systemperiode um 2x wächst. Wir stellen 
E, und E,, als Funktionen von 9 dar und teilen beide in je zwei Anteile: 


E,@) =9, +48), En) =O9n+AEn)....... (2) 


9, bzw. 9 ,, die zu ermittelnden Trägheitsmomente der zusätzlichen Schwungmassen am Motor 
bzw. an der Arbeitsmaschine, AE,(p) und AE,,(p) sind nunmehr als bekannt vorauszusetzen). 
Aus (1) und (2) erhält man: 


T-[O+ ABS. ..... ee 
mit 
od, = ABA) + AED)... 
Übersetzungsverhältnis u = d/dır- 
Bei einer Verdrehung des Systems wird von diesem die Arbeit geleistet: 
dA=MI) HI, —- WON) Here. Fa t0)% 


Dabei ist M(d,) das Moment an der Motorwelle, das man erhält, wenn Den die auf die Kolben 


wirkenden Gaskräfte, Reibungen und evtl. Schwerkräfte auf das Wellenmittel reduziert („effek- 


tives Drehmoment“), W(9,,) das analog gebildete Moment für die Welle der Arbeitsmaschine. 
Durch Einführen der Winkelkoordinate @ geht (5) über in 


d4 =k, | — I Wo] dp = 9) dr REN MG 


(k)... Anzahl der Drehungen der Motorwelle in der vollen Systemperiode i 
Ungleichförmigkeitsmoment‘‘). si R ee 


Der Arbeitssatz für das System lautet in differentieller Form 
w* 
alte + AE(p)] 3) = dA: = Od (2): 


1) Der vorliegende Bericht bezieht sich auf eine Arbeit, die i i i i i 

i i ; gemeinsam mit G. Heinrich, W - 
geführt wurde: Ein neues graphisches Verfahren zur Schwungradermittlung bei Kolbenmaschinen: An, 
u. Wärmewirtschaft, 8. Jg., Heft 9, 1953, S. 241— 247 Springer-Verlag, Wien. £ 
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Wir bezeichnen die Werte von 9, die den Extremwerten von w und 7 entsprechen in folgender 
Weise: 


Be OT Oma SIE ee 
Die Integration von (7) zwischen g, und 9, ergibt 
Omaz min fe 
[0 + AEG [OH Od... 8). 
Pı 
Für Oin und Opa, erhält man Näherungswerte durch den üblichen Ansatz 
0) 6) 
Omin = Om (1 >>> =) 3 er — Om (1 + a DEI RE TE ARD (9) 


mit der mittleren Winkelgeschwindigkeit 


T 
1 
m | o(t) dt 
N ; 
und dem zulässigen Ungleichförmigkeitsgrad 


se Omas — min 
Om 
(* bedeutet hier und im Folgenden Näherungswerte). 

Zur Ermittlung der Stellen 9, und 9, suchen wir vorerst die Stellen 9, und o, auf. Dort hat, 
wie (7) zeigt, die Funktion Q(p) Nullstellen. ‚Wegen der Veränderlichkeit von E, und E,, werden 
die Stellen 9, und 9, bzw. 9, und 9, nicht zusammenfallen, aber benachbart sein. Wir führen die 
Differentiation in (7) aus 
dAE(p) w: do 

a, 5 + [@ are 


AN 
2 LE TE Darstellung des Zeichenaufwandes 


für die praktische Durchführung der Schwung- 
N radermittlung 


» _ 1 /29@) 
ei fat) N UN 
| dp rag 
-+ 


Zur Ermittlung der Stellen 9] bzw. »% (Näherungswerte für @, bzw. @,) hat man also nur in den 
an die Stellen 9, bzw. $, anschließenden kurzen Intervallen, in denen Q(p) und d AE/dg gleiches 


200) 
Vorzeichen besitzen, die Kurvenstücke Ka zu zeichnen und mit den Werten on 
+ 


dAE(p) 
dp 


bzw. &%., zum Schnitt zu bringen. Im übrigen Bereich braucht man die Getriebefunktionen 
nicht zu ermitteln, ein Umstand, der große Arbeitsersparnis bedeutet. Die Fläche unter der 
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wird. Im Bedarfsfalle können alle mit Stern versehene 
beit!) angegebenen Verfahren kontrolliert und durch Ite 
Die angegebene Methode ist auch für den Fall verwendb 9 gegeben R 
Das Bild zeigt die Durchführung des Verfahrens für ein X stem Kolbenme ’ 
kompressor (Viertakt-Einzylindermotor 1500 Umdr./min, einfachw rkende zylinder-Kolben- 
kompressor, Übersetzungsverhältnis u= 1/2, zuläss. Ungl. Grad ö = 0,05 leichzeit: Far 


AN 


den nötigen Zeichenaufwand dar. Die zu planimetrierende Fläche ist schraf 1 bi 
Eine parallel zu obigem Verfahren ausgearbeitete rein rechnerische Methode 
erwähnt 3). x 


Zur Stabilität von Reglern mit Nachlaufzeit 
Von Wolfgang Hahn in Braunschweig 
Es sei die Gleichung des Regelsystems auf die Gestalt!) 
YW@+2KY$ + YW—alt—N=0...........( 


gebracht. Partiallösungen von (1) sind die Funktionen e%!, wobei A; eine Wurzel der transzendenten 
Gleichung e 
AQ) EI ak ol ae re ee (2) 


n Y 2 
ist. In der einschlägigen Literatur wird als notwendige und hinreichende Bedingung für die Stabi- 
lität des durch eine Gleichung vom Typ (1) beschriebenen Systems angesehen, daß alle A, negative 


Realteile haben, ohne daß dies Kriterium bisher mathematisch exakt begründet worden ist. Zum - 


Beweis muß man folgendes zeigen: 1. Jede Lösung von (1) läßt sich in eine Reihe der Gestalt 
3 0; e*i entwickeln. 2. Die Reihen konvergieren auf der ganzen positiven t-Achse gleichmäßig. 
(Da bei speziellen Anfangsbedingungen einzelne der Koeffizienten C, verschwinden können, ist 
die fragliche Bedingung streng genommen nur hinreichend.) 


Für die Standardlösung der Regeltechnik, die durch die Anfangswerte YW)=0(<O0), 
-Y (0) = 1 gekennzeichnet ist, erhält man aus (1) durch Laplace-Transformation #7 
x y() Al) =1 


und daraus durch die komplexe Umkehrformel und Berechnung des Integrals mit Hilfe des Re- 
siduensatzes zunächst formal 


YO) > Kan en EN 8). 


(Dabei ist von etwaigen Mehrfachnullstellen, die höchstens in endlicher Anzahl auftreten können, 
abgesehen.) Drückt man die Exponentialgrößen e*: mittels (2) durch die Potenzen von 4; aus 


und beachtet, daß A() eine ganze Funktion von der Ordnung 1 ist, so zeigt sich, daß die Reihe -» 


(3) und ihre zwei ersten formalen Ableitungen absolut konvergieren. Daraus folgt, daß (3) eine 
Lösung ist; daß sie die „richtigen“ Anfangswerte hat, muß man durch direkte Berechnung zeigen. 
Wegen der absoluten Konvergenz der Reihe ist die geforderte gleichmäßige Konvergenz ge- 
sichert. — Aus der Standardlösung läßt sich in üblicher Weise durch Differentiation bzw. Faltung 
die Lösung der inhomogenen Gleichung mit beliebigen Anfangswerten gewinnen. 


Der Beweis läßt sich in gleicher Weise für Gleichungen von höherer Ordnung führen, sofern 
die höchste auftretende Ableitung Y® (t) nicht mit Verzögerungen behaftet ist. Tritt ein Glied 
Y®=D (€ — 7) auf, so muß man in (3) die konjungiert komplexen Glieder zusammenfassen. 


?)R.v.Mises, Die Ermittlung der Schwungmassen im Schubkurbelgetri 
eo, g gmassen i chubkurbelgetriebe. Z. Österr. Ing. u. Arch. 

®)K.Desoyer u. A. Slibar, Die rechnerische Ermittlung des Ungleichförmigkei i - 
maschiren. Österr. Ing.-Arch. VII, S. 100 (1953). 5 a 


(1947) Han Collatz, Zur Stabilität von Regelungen mit Nachlaufzeit. Z. angew. Math. Mech. 25/27, 60—62 
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Einen starren Rotor, d.h. einen Rotor der weit unterhalb seiner ersten kritischen Dreivahı 
läuft, kann man bekanntlich immer dadurch auswuchten, daß man in zwei beliebigen Ebenen 


‚senkrecht zur Achse zusätzliche Massen anbringt bis die Drehachse zu einer Hauptträgheitsachse 


des Rotors wird. 

Ist der Rotor elastisch, so genügt dies nicht; denn durch die unbekannte Unwucht und die 
zusätzlichen Massen können Biegemomente entstehen, die den Rotor so stark verbiegen, daß er 
bei voller Drehzahl ganz aus dem SIERT kommt, trotzdem er bei niedriger Drehzahl voll- 
kommen ausgewuchtet war. 

Mit den Orthogonalitätsbedingungen, welche die elastischen Linien bei den verschiedenen 
kritischen Drehzahlen miteinander Diner: kann ein besseres Verfahren angegeben werden. 


Es sei: 


 Ym(x) Elastische Dede Rotors bei der kritischen Winkelgeschwindigkeit as 
Y,(2) Elastische Linie des Rotors bei der kritischen Winkelgeschwindigkeit , , 
m(x) Rotormasse je Längeneinheit 


so gilt die Bedingung: 
A i 1 

(om — ©) | Ym‘ Yn‘ m da = 0 
ö 


und daraus folgt für das Integral: 


>) mM=N 
IE ale Io m<n' 
0 


Das Integral läßt sich wie folgt deuten, wenn man esz. B. mit co, multipliziert: 
Die Arbeit, welche die Fliehkräfte om y%, m auf die Durchbiegung Y., ausüben, ist stets Null. 


Daraus folgt, daß eine zusätzliche Massenmomentverteilung, welche mit y, m proportional 
ist ausschließlich auf die kritische Drehzahl », einen Einfluß haben kann, die übrigen kritischen 
Drehzahlen jedoch überhaupt nicht beeinflußt. 

Daraus ergibt sich das folgende Auswuchtverfahren für elastische Rotoren. Man bringt den 
Rotor bei der ersten kritischen Drehzahl durch eine zusätzliche Massenmomentverteilung, welche 
zu Y,-m proportional ist, ins Gleichgewicht. Nachdem dies erreicht ist, geht man auf die zweite 
kritische Drehzahl und wuchtet den Rotor durch eine zu %,: m proportionale zusätzliche Massen- 
momentverteilung aus. Durch diese zweite Massenmomentverteilung kann das Gleichgewicht 
bei der ersten kritischen Drehzahl nicht gestört werden, wie die Orthogonalitätsbedingung es 
zeigt. Man geht dann auf die dritte kritische Drehzahl usw. 

. Da die Durchbiegung des Rotors bei beliebiger Drehzahl als Summe der kritischen Durch- 
biegungen — mit von der Drehzahl abhängigen Konstanten multipliziert — dargestellt werden 
kann, so wird ein Rotor, der bei sämtlichen kritischen Drehzahlen ausgewuchtet ist, auch bei be- 


 liebiger Drehzahl ausgewuchtet sein. 


Praktisch genügt es, den Rotor bei denjenigen kritischen Drehzahlen auszuwuchten, die 
unterhalb der zulässigen Höchstdrehzahl liegen. Eventuell kann man bei der Höchstdrehzahl 
noch mit derjenigen Massenverteilung, welche der nächsthöheren kritischen Drehzabl entspricht, 
nachwuchten. 

Da praktisch nur wenige kritische Drehzahlen unterhalb der Höchstdrehzahl liegen, ist es 
streng genommen nicht nötig, die theoretisch geforderten Massenmomentverteilungen einzuhalten. 

Es genügt pro kritische Drehzahl grundsätzlich eine Auswuchtebene; bei zwei kritischen 
Drehzahlen unterhalb der Höchstdrehzahl, und Nachwuchten bei der Höchstdrehzahl also drei 
Ebenen, die an sich beliebig gewählt werden können, praktisch jedoch so weis auseinander liegen 
sollten wie nur möglich. 


Die Durchbiegungen der drei kritischen Formen in den drei Ebenen seien 


Yı1 Yıa Yıs 
Yaı Ya2 Yaz 
Yzı Y32 Ya3 


der erste Index bezeichnen die kritische Drehzahl der zweite die Auswuchtebene. 
Die 3 x 3 Massenmomente, die für das Auswuchten bei den drei kritischen Drehzahlen nötig 
sind, ergeben sich dann aus den folgenden Gleichungen: 
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Für die erste kritische Drehzahl: 


M;ı Yu + Mi Yı2 + M;s Yıs = 4 
M,ı Yı + Mıs Ye2 + Mıs Yes = 0 
M,ı Ysı + Mıs Ys2 + Mıs Ya = 9: 


Für die zweite kritische Drehzahl: 


My Yın + Ma Yız + Mas Yıs = 0 

M;ı Yaı + Ma Ya2 + Mas Y 25 = b 

M;ı Ysı + Mae Yss + Mas Ya; = 0. 
Für die dritte kritische Drehzahl: 


M;ı Yıı + Mag Yız + Mas Yıs = 0 
M;, Yaı 1: M;; Ya3 1 M;; Yaz = 0 
M;ı Yaı + Mag Ya + Mas Yas = ©» 


Die Konstanten a, b und c sind an sich beliebig, und können z. B. durch 1 ersetzt werden. 

Es ist vermutlich am zweckmäßigsten, in frei beweglichen, federnd aufgehängten Lagern 
auszuwuchten. Denn in solchen Lagern entsprechen die beiden ersten kritischen Drehzahlen dem 
üblichen ‚statischen‘ und ‚dynamischen‘ Auswuchten. Diese beiden Auswuchtungen sind aber 
unbedingt notwendig, bevor man es überhaupt wagen kann einen größeren Rotor auf volle Dreh- 
zahl zu bringen. Auf freien Lagern schließt sich das ‚‚elastische‘‘ Auswuchten dann zwanglos an; 
beim elastischen Auswuchten auf festen Lagern ist dies nicht der Fall. 

Die elastischen Linien der kritischen Drehzahlen kann man sich heute durch ein elektrisches 
Analogie-Gerät oder mit Hilfe einer elektrischen Rechenmaschine verhältnismäßig einfach ver- 
schaffen. 


Gesetzmäßigkeiten und Grenzwerte von Drehschwingunssfrequenzen 
Von H. Neuber in Dresden 


Während für die kritischen Drehzahlen biegesteifer Wellen mit Schwungmassen einfache 
Beziehungen für obere und untere Schranken bekannt sind (u. a. dieFormelvon Dunkerley), 
kennt man bisher solche Grenzwerte für das Torsionsschwingungsproblem noch nicht. Für über- 
schlägige Rechnungen mit wenig Zeitaufwand, sowie für Kontrollmöglichkeiten bei der genauen 
Berechnung von Drehschwingungssystemen mit beliebig vielen Massen ist aber die Kenntnis der 
Frequenzschranken von großem Vorteil. Durch Einführung der elastischen Torsionsmomente 


n=2 n=3 n=4 n=5 n=6 n=2Z n=6 


wet 
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Bild 1. Die Frequenzdeterminante bei Verwendung der elastischen Torsionsmomente als Unbekannte 


En 
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anstelle der Drehwinkel als Unbekannte wird eine Differential-Differenzgleichung aufgestellt, 
deren Lösung eine besonders übersichtliche Determinante liefert (Bild 1). Die weitere Unter- 
suchung dieser Determinante, deren Verschwinden die Frequenzgleichung gibt, läßt anschauliche 
Zusammenhänge des betrachteten Schwingungssystems mit einer Reihe von Ersatzsystemen 
erkennen, welche aus dem Hauptsystem durch teilweise Erstarrung (Bild 2) oder teilweise 
Wegnahme von Steifigkeiten (Bild 3) herstellbar sind. Daraus ergibt sich die Möglichkeit, die 
Koeffizienten der Gleichung (n— 1)-ten Grades für das Frequenzquadrat durch die Eigen- 
frequenzen der Ersatzsysteme darzustellen. 
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(3-Massen- System) 


Bild 2. Zurückführung der Koeffizienten der Frequenzgleichung auf Eigenfrequenzen von Ersatzsystemen, die 
durch teilweise Erstarrung aus dem Hauptsystem hervorgehen 
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Bild 3. Inverse Auflösung der Determinante und Darstellung der zugehörigen Koeffizienten durch Eigen- 
frequenzen von Ersatzsystemen, die durch teilweise Wegnahme von Steifigkeiten entstehen 


So läßt sich eine Reihe von Sätzen für den Zusammenhang der Eigenfrequenzen des Haupt- 
systems mit denen der Ersatzsysteme aufstellen. 


Insbesondere liefert die Beziehung 


1 1 
ah in] nnd) 
einen unteren Grenzwert für die Frequenz der Grundschwingung: 
ze 
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_ Für das Torsionsschwingungsproblem wi 


Sachverhalt der Dunkerleyschen Fı 
Andererseits besteht die Beziehung, 


welche einen oberen Grenzwert für die höchste 


Schwingungsberechnung von Stabwerken a 
Von E. Pestel in Hannover ee 


Obwohl es naheliegt, den Grundgedanken von Holzer zur Berechnung der Torsions- 
schwingungen von Maschinenwellen auch auf die Ermittlung der Längs- und Biegeschwingungen 
von Stäben anzuwenden, hat erst vor. 10 Jahren Myklestad für den Fall reiner Biege- 
schwingungen eine dem Holzer-Verfahren verwandte Rechenmethode entwickelt. Im fol- 
genden soll nun die Vorgehensart von Holzer auf allgemeine Stabwerke mit diskreter Massen- 
und Massenträgheitsmomentenbelegung angewendet werden. Dabei braucht der Schwerpunkt 
der Massen nicht auf der elastischen Achse der Federstäbe zu liegen, die selbst als masselos an- , 
genommen werden und gleichzeitig Längs- und Querkräfte sowie Biege- und Torsionsmomente E 
aufnehmen können. Wenn die Eigenschaften der Massen und Federn bekannt sind und aus dem 
interessierenden Frequenzbereich eine Frequenz vorgegeben wird, ist es möglich, aus der Kenntnis 
der 6 Auslenkungen und der entsprechenden 6 Schnittgrößen, sagen wir unmittelbar rechts von 
der Masse r die entsprechenden 12 Größen unmittelbar rechts von der Masse r + 1 zu berechnen 
(vgl. Bild 1); denn jede der 12 Größen an der Stelle r + 1 läßt sich dann als Linearkombination 


Bild 1 


der 12 an der Stelle r gefundenen Größen ausdrücken. In übersichtlicher Weise kann man somit 
die Übertragung der 6 Auslenkungs- und der 6 Schnittgrößen von einer Stelle r zum folgenden 
Punkt r + 1 als Produkt einer quadratischen Matrix und einer Spaltenmatrix darstellen. 

Wir können also bei einem freien mit diskreten Massen belegten Stabzug (vgl. Bild 1), die 
aus 12 Elementen bestehende Spaltenmatrix ıt, durch folgende Matrizenmultiplikation bestim- 


men. 
1, Se N ee . EL bg 
Darin stelle die quadratischen Matrizen X, die Feder- und Masseeigenschaften jedes einzelnen 
Abschnittesr dar. Da an jedem Ende eines solchen räumlich beanspruchten Stabzuges i.a. 6 in 
den 12 Größen homogene Randbedingungen gegeben sind, können von den 12 Größen, die 
unmittelbar rechts von der äußersten linken Masse vorhanden sind, 6 durch die übrigen 6 aus- 
gedrückt werden. Aufgrund der Tatsache, daß sich am rechten Ende des Stabzuges die vor- 
erwähnten 12 Größen als Linearkombinationen der 6 am Anfang frei wählbaren Größen ergeben, 
liefern die hier zu erfüllenden 6 Randbedingungen ein System von 6 homogenen linearen Glei- 
chungen mit den 6 an der Stelle 1 frei gewählten Größen als Unbekannten, die nur dann von 
Null verschieden sein können, wenn die Nennerdeterminante dieses Gleichungssystemes ver- 
schwindet. Der Wert der Nennerdeterminante ist aber eine Funktion der eingangs der Rechnung 
gewählten Frequenz. Nur, wenn zufällig eine Eigenfrequenz in die Rechnung eingeführt wird, 
ergibt sich die Nennerdeterminante zu Null. Man berechnet nun in dem interessierenden Frequenz- 
bereich die Nennerdeterminante für einige Frequenzen und findet durch Interpolation die 
gesuchten Nulldurchgänge. 
Bei verzweigten Stabwerken sind die Verhältnisse nicht so einfach wie bei dem Stabzug 
des Bildes1. Bei solchen Stabwerken erweist es sich als notwendig, zunächst das Stabwerk in ein- 
zelne freie oder federgefesselte Stabzüge zu zerschneiden. Die Schnittgrößen, welche die Wirkung 


!) Eine ausführlichere Darlegung wird demnächst in der Zeitschrift „Ing.-Archiv“ veröffentlicht, 
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der abgeschnittenen auf den jeweils betrachteten Stabzug wiedergeben, sind dann als Unbe- 
kannte einzuführen. Die Forderung nach dem kinematischen Zusammenhalt des Stabwerkes 
und die Betrachtung des Axioms ‚„‚actio gleich reactio“ liefern die zusätzlichen Gleichungen für 
die Bestimmung dieser Unbekannten. Es wurde im Vortrag gezeigt, daß die gesamte Schwin- 
gungsrechnung durch Verwendung geeigneter Matrizenschemata einheitlich programmiert wer- 
den kann, so daß mit der Durchführung der Zahlenrechnung entweder ungelernte Kräfte oder 
Rechenautomaten betraut werden können. Das Verfahren bietet ferner die Möglichkeit, in die 
Rechnung Parameter einzuführen. Auf diese Weise können z.B. die Möglichkeiten der Ver- 
lagerung von Eigenfrequenzen durch Änderung der elastischen Eigenschaften oder (und) der 
Massenbelegung im Hinblick auf eine gegebenenfalls notwendig werdende Verstimmung des 
Stabwerkes beurteilt werden. - 

Eine Variante des vorgetragenen Verfahrens erlaubt, mit einer einzigen Durchführung der 
Matrizenproduktbildung entsprechend (1) die dann allerdings sehr umfangreiche Determinante 
für die Berechnung des gesamten Frequenzspektrums des Stabwerkes aufzustellen. 

Schließlich ist das Verfahren auch geeignet, erzwungene ungedämpfte und gedämpfte 
Schwingungen von Stabwerken zu berechnen. Eine ausführliche Darstellung mit Beispielen ist 
in der Arbeit des Verfassers „Ein allgemeines Verfahren zur Berechnung freier und erzwungener 
Schwingungen von Stabwerken‘“ in den Ende Sommer 1954 erscheinenden Abhandlungen der 
Wissenschaftlichen Gesellschaft der Technischen Hochschule Braunschweig gegeben. 


Bemerkungen zu einer Arbeit von Caratheodory 
| Von R. Gran Olsson in Trondheim 


In einer Arbeit über die Bewegung eines Schlittens als ein nicht-holonomes System macht 
C.Carath&odory die Annahme, daß eine Kraft R in Querrichtung auf den Schlitten 
wirkt, d.h. senkrecht zur Bewegungsrichtung [1]. Die Bewegung möge in der Horizontalebene 
stattfinden, und die Reibung an den Kufen soll durch den Zug des Pferdes in entgegengesetzter 
Richtung dauernd aufgehoben sein. 

Es dürfte indessen mit den tatsächlichen Verhältnissen besser übereinstimmen, wenn ein 
Moment statt einer Kraft am Schlitten wirkend angenommen wird. Es sei F die Kraft, die an den 
Kufen des Schlittens in Richtung der Bewegung von gleicher Größe wie die Reibung wirkt. Die 
Kufen mögen von verschiedener Rauhigkeit angenommen sein, so daß die Reibungskräfte F, und 
F, an den Kufen verschieden sind und die Resultante (F, + F,) eine zur Symmetrieachse des 
Schlittens exzentrische Lage hat, während der Zug des Pferdes F längs der Symmetrieachse 
wirken möge. Wenn nun #=F, + F, ist, wirkt keine exzentrische Kraft sondern ein Moment 
M=F-e am Schlitten, wo e die Exzentrizität der Resultante (F, + F,) in bezug auf die 
Symmetrieachse des Schlittens bedeutet. 

Dies Moment mag konstant angenommen werden, wonach die Geschwindigkeiten des 
Schlittenschwerpunktes sich durch einfache Integrationen ergeben. Der Hodograph der Be- 
wegung ist durch Ellipsen (ähnlich wie bei der Planetenbewegung) gegeben. Die Schwerpunkts- 
bahnen ergeben sich durch eine weitere Integration und zwar als Fresnelsche Integrale, die 
bekanntlich tabuliert sind [2]. 

Eine andere Annahme für das Moment ist die der Proportionalität zwischen Moment und 
Geschwindigkeit in Richtung der Bahntangente des Schlittens. Dies bedeutet wiederum Pro- 
portionalität zwischen der Winkelbeschleunigung und der Geschwindigkeit. In diesem Falle ist 
eine Integration der Bewegungsgleichungen möglich, und zwar ergibt sich die Zeit durch ein 
elliptisches Integral, wo im Integranden die Winkelgeschwindigkeit als Veränderliche auftritt. 

Durch Umkehrung ergibt sich dann die Winkelgschwindigkeit als elliptische Funktion der 
Zeit. Mit Hilfe der Ausgangsgleichungen ergeben sich daraus die Komponenten der Geschwin- 
digkeit. 

3 Die ArbeitvonCarath&odory kannübrigens auch nach der Richtung ergänzt werden, 
daß statt der seitlichen Geschwindigkeit v = 0 diese konstant, aber von Null verschieden ange- 
nommen wird. Die Integration wird ähnlich wiebeiCarath&odory elementar, und es lassen 
sich ähnliche Bahnen für die Bewegung des Schlittens konstruieren wie in der erwähnten Arbeit 
angegeben. 

Verf. hofft bei einer späteren Gelegenheit ausführlicher auf die hier genannten Fälle einer 
Erweiterung der Carath&odory schen Arbeit zurückzukommen. 

[1] C©. Garath&odory, Z. angew. Math. Mech. 13, 71—76 (1933). Siehe auch A. Sommerfeld ‚Vor- 


lesungen über theoretische Physik, Bd. I Mechanik, Leipzig 1947, 3. Auflage, S. 220— 229. 
[2] E.Jahnke u, F.Emde ,Funktionentafeln mit Kurven und Tafeln. Leipzig u. Berlin 1938, S. 35 
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Ein neues Strömunsssgleichnis der Torsion 
Von E. Pestel in Hannover 


Im Gegensatz zu der bekannten hydrodynamischen Analogie von Thomson stellt das 
hier mitgeteilte Strömungsgleichnis eine mathematische Beziehung zwischen dem Verhalten 
einer viskosen Flüssigkeit und dem Torsions- 
verhalten eines prismatischen Stabes her. Dabei 
wird folgender Strömungsvorgang betrachtet. 

Eine Scheibe beliebigen Umrisses sei ge- 
lagert auf einer Flüssigkeitsschicht, deren Träger 
eine unendlich große Ebene ist (vgl. Bild ]). 
Scheibe und Ebene sind parallel zueinander. Wird 
bei einer Relativbewegung längs ihrer Normalen 
die Spalthöhe A zwischen Scheibe und Ebene 
z.B. verringert, so wird Flüssigkeit aus der 
Schicht zwischen Scheibe und Ebene herausge- 
drückt. Dadurch entsteht in dem Spalt ein 
i.a. von Ort zu Ort verschiedener Überdruck 9 
gegenüber dem Druck p, in der Flüssigkeit außer- 
halb der Scheibe. Es läßt sich nun u.a. zeigen, 
daß die Kraft, welche die beschriebene Bewe- 
gung der Scheibe hervorbringt, dem Drillmoment 
zugeordnet werden kann, das einen prismatischen 
Stab von einem Querschnitt gleich dem Umriß der Scheibe verdrillt. 

Unter den Voraussetzungen der St. Venantschen Theorie wird das Torsionsverhalten 
eines prismatischen Stabes, dessen Querschnitt im allgemeinen Fall mehrfach zusammenhängend 
ist, durch folgende Gleichungen hinreichend beschrieben. 


Bild1 


° ö 
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Die darin verwendeten Bezeichnungen bedürfen wohl im Hinblick auf Bild 1 keiner weiteren 
Erläuterung. Diesen Gleichungen des Torsionsproblemes können wir unter den Voraussetzungen 
für die sogenannte schleichende Strömung folgende Gleichungen für das soeben beschriebene 
Strömungsproblem gegenüberstellen. Hierin bedeuten q die in der Zeiteinheit strömende Flüssig- 
keitsmenge, deren Strömungsrichtung jeweils durch den Index gekennzeichnet ist, u die absolute 
Zähigkeit der Flüssigkeit, während die Bedeutung der übrigen Bezeichnungen dem Bild1 zu 
entnehmen ist. 
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Es ist beabsichtigt, dieses Strömungsgleichnis experimentell zu verwirklichen, zumal die Aus- Dre 
sicht besteht, auf diese Weise die Drillsteifigkeit von Stäben mit mehrfach zusammenhängendem 
Querschnitt auf dem Versuchswege besonders einfach zu bestimmen. Darüber hinaus bietet aber 
das Strömungsgleichnis wegen seiner Anschaulichkeit die gleichen Vorteile in der gedanklichen 
Bewältigung gewisser Torsionsprobleme wie das bekannte Prandtlsche Seifenhautgleichnis 
und die eingangs erwähnte hydrodynamische Analogie vonThomson. _ 


Das d’Alembertsche Prinzip in Motorsymbolik und seine 
Anwendung auf Stoßprobleme © 

Von W. Raher in Wien 2 E 

Im Rahmen der Motorrechnungt!) geht die reaktionskräftefreie Behandlung eines nicht- Bene: 


dissipativen Systems von n starren Körpern, das bestimmten Bedingungen unterliegt, vom 
d’Alembertschen PrinzipinLagrangescher Fassung aus?) 


I e- On 4 


Darin bedeutet $%x den Impulsmotor des k-ten Körpers, $tx die Schraube der auf ihn wirkenden 
eingeprägten Kräfte und 6x eine virtuelle Lagenänderung dieses Körpers. Allgemein ist der 
Motor der Lagenvariation 6x = (W;%) aus dem Geschwindigkeitsmotor & = (W;) durch 
Ersetzen der Differentialquotienten nach der Zeit durch solche nach einem Variationsparameter & 
entstanden zu denken; die Zahl der in den@x von (1) enthaltenen willkürlichen Variationsgrößen 
ist den Systembindungen entsprechend eingeschränkt. 

Die Spezialisierung von (1) auf Momentankräfte führt auf 
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% bezeichnet den Impulsmotor nach dem Stoß und © die Stoßschraube. Aus (2) sind leicht 
allgemeine Sätze aus der Theorie des Stoßes herleitbar (Carnotsche Theoreme für Systeme 
starrer Körper usw.), hier sei aber der praktische Ansatz von Stoßproblemen an Hand einfacher 
Beispiele gezeigt. 

Für das in Bild 1 dargestellte holonome System, bei dem ein Körper gelenkig gelagert, der 
zweite mit diesem durch ein Kugelgelenk verbunden ist, folgt für die Lagenvariationen bezogen 
auf O, bzw. O0, 


6, — (Bi; 0), 6, = (BB; DR, x. 


Wird auf den oberen Körper ein Stoß ausgeführt, so gilt: ©; = 0 und &, ist durch die Stoßkraft 
und ihr Moment bezüglich O, bestimmt: ©, = (6, My,). Nach Einsetzen in (2) und Ausführung 
der Produkte sind die Koeffizienten der beiden willkürlichen Vektorvariationsgrößen %, und ®, 
Null zu setzen. Die so erhaltenen 6 Skalargleichungen bestimmen den Geschwindigkeitszustand 
nach dem Stoß. 

Völlig analog gestaltet sich der Ansatz bei nichtholonomen Systemen, sofern nur die Bin- 
dungen linear in den Geschwindigkeiten sind. Bei nebenstehend skizziertem System, einem 
Räderpaar, das auf horizontaler Ebene reine Rollbewegungen ausführen soll, haben wir bei 5 Frei- 


I) R.v. Mises, Motorrechnung, ein neues Hilfsmittel der Mechanik. Z. angew. Math. Mech. 4, 1924. 
2)W. Raher und F. Selig, Die Verwendung der Motorsymbolik in der theoretischen Mechanik. Er- 
scheint demnächstin den Sitzungsberichten der Österreichischen Akademie der Wissenschaften. 
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vorher in Ruhe, so folgen nach Bildung der Produkte 

aus (2) durch Nullsetzen der Koeffizienten der beiden 
el Winkelvariationen die Gleichungen, 
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die gemeinsam mit (3) den Geschwindigkeitszustand nach dem Stoß festlegen. Wenn die Ebene 
nicht unendlich rauh ist, hat man nachträglich zu überprüfen, ob die nichtholonomen Bedingungen 
überhaupt physikalisch aufrechtzuerhalten sind. Aber auch dann wird man dieses und ähnliche | 
Probleme zuerst einmal so ansetzen. x 

Falls die Stoßkräfte von vornherein nicht gegeben, sondern durch stoßende Körper auf das 
System übertragen werden, unterteilt man in Kompressions- und Restitutionsperiode und setzt 
die Innenprodukte der Schrauben der Ausgleichsgeschwindigkeiten in den Motor der gemeinsamen 
Stoßnormalen gleich. Auf diese Art gelangt man beim einpunktigen Stoß zweier Körper rascher 
zum Resultat, als bei Bestimmung des Betrages der Stoßkraft aus energetischen Überlegungen. 

Auch die plötzliche Einführung skleronomer oder rheonomer Bindungen in ein System ist 
unmittelbar im Anschluß an (2) zu behandeln. 

Somit sind durch (2) verschiedenste Stoßprobleme auf einen gemeinsamen Ansatz zurück- 
geführt, von dem man, besonders bei räumlichen Problemen, vorteilhaft ausgeht. 


Über automatisch gesteuerte Bewegungen 
Von Stefan Schottlaender in Würzburg 


In verschiedenen Zweigen der Technik treten bei Störung einer Hauptbewegung Schwin- 
gungen auf, deren Eigendämpfung den Forderungen nach schnellem Abklingen nicht genügt 
oder sogar fehlt. Um in diesen Fällen dennoch rasch auf kleine Amplituden zu kommen, hat man 
automatische Regelungen oder Steuerungen geschaffen, die die Störbewegung möglichst schnell 
zum Verschwinden bringen sollen. Solche Steuerungen können entweder stetig oder unstetig 
arbeiten. Die unstetig arbeitenden haben hierbei den Vorteil apparativer Einfachheit, größerer 
Robustheit und geringerer Herstellungskosten, ihr Nachteil dagegen ist, daß sie von sich aus neue 
unerwünschte Störungen am schwingenden System hervorrufen können, wenn man die Steuer- 
parameter unzweckmäßig auswählt. Aber noch eine dritte Möglichkeit hat sich in neuester Zeit 
in der Praxis als außerordentlich erfolgreich erwiesen: die „gemischte‘‘ Steuerung!). In hin- 


A,W. H.Bihr Mi ps ‚, Graphical Solution of some Automatic-Control Problems involving Saturation Effects 
with Application to Yaw Dampers for Aircraft. NACA, TN 3034, Washington, USA, 1953. 
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mit der " Steuerfunktion 


gewst+ße. g=0 ee MR (2). T 


Dresehretet aber ll oder |5| einen gewissen Maximalbetrag, so verbleibt (z.B. als Folge von 
„Sättigungseffekten‘‘) die steuernde Kraft bzw. ihre Änderungsgeschwindigkeit auf einem Höchst- 


G=+G% (Stellungszuordnung) ey ver (3a) 
oder 


G=K,+tVt (Laufgeschwindigkeitszuordnung) ee OD 


also eine nichtlineare unstetige Steuerung. Für in dieser Art gesteuerte Bewegungen kann es 


keinen geschlossenen analytischen Ausdruck (im Großen) geben, mit Ausnahme der periodischen 
Bewegungen. Diese lassen sich durch Aufstellung eines einfachen linearen Gleichungssystems 
vollständig berechnen und von ihnen ausgehend können dann auch Aussagen über allgemeine 
Bewegungen gemacht werden, insbesondere die für die Technik so wichtigen Stabilitätsaussagen. 
Die Berechnung der periodischen Bewegungen einer ‚gemischten‘ Steuerung nach (2) und 
(3b) erfolgt nun nach einer Methode vonBilharz?). Die Steuerung verlaufe schon „gesättigt‘“, 
an die Stelle von g trete die tatsächliche periodische Steuerkraft (Periode 2a), (# = 0, 1,2, - - -) 


TE in I,:2va<t<(2v+1)a & 
a NE EDER ea i 
Lösungen von (1) mit der rechten Seite — @(t) sind dann 
_ fAı,e"+ B,, ea! —G(t) +22V, in En er 
Ar a ee u 
Hierbei sind A,, A, die Lösungen der charakteristischen Gleichung. Durch die Forderungen 


z(2» +)a—0) =z2(2» +N)a+0) 
&(2» +1) a0) =& (2»+1)a+0) 
x(2va +0) 2 (2» +29)a—0) 
&(2va +0) =: (2, +2)a—0) 
erhält man vier lineare Gleichungen, aus denen sich Ay, Aa, Bi,, Ba, eindeutig bestimmen lassen. 
Zu jedem «> 0 gibt es daher eine und Aur eine periodische Lösung der Periode 2a, insbesondere 
gibt es zu jedem a> O0 ein und nur ein Anfangswerttripel (&,, &%g &,) und zu verschiedenen Wer- 
ten a,, a, auch verschiedene Anfangswerttripel. 
Die Steuerung nach (4) muß nun aber noch der Steuervorschrift angepaßt werden, die be- 
sagt, daß die Umschaltpunkte t = n a der „gesättigten‘“ Steuerung auf der Steuerkurve (2) liegen 
müssen, d.h.k, =a2, + ß&, und ku, + V,a =ax(a) + P& (a), woraus folgt: 


|Stetigkeit in den Schaltpunkten = (2»+1)a, 
. h 
(6). 


I 


|Perioaizität in den Schaltpunkten t=2»a 


u + =— 5 V, (1. Steuerungsbedingung) nae1 ar 2% (7). 


Dies stellt eine einparametrige Geradenschar in der Ebene der Steuerparameter &, ß dar, be- 
stimmte Gebiete dieser Ebene sind hiernach für periodische Lösungen ausgeschlossen. Weiterhin 
ist zu fordern, daß |9] > V, ist zur Zeit 2?=0 (d.h. für den ersten Schaltpunkt, sonst würde ja 
keine „‚gesättigte‘‘, sondern eine lineare Steuerung beginnen), das bedeutet: 


la&o +Bßäol> Vo (2. Steuerungsbedingung)... . .... (8). 


Schließlich dürfen zwischen zwei aufeinanderfolgenden Schaltpunkten na und (n + 1)a keine 
weiteren Schnittpunkte von@ mit der linearen Steuerkurve g liegen (diese müßten sonst ja auch 
Schaltpunkte sein); dies führt im Falle E 2 <1 auf die Einschränkung 


0<a< er 2 (3. Steuerungsbedingung) .  .. .... (9). 
Mit (7) genügen die periodischen Lösungen von (5) noch der Funktionalgleichung 
x(t +na) = (— 1)* x(t) RE rn (10); 


derselben Funktionalgleichung genügen auch & (t) und & (t). 


2) H. Bilharz, Über eine gesteuerte eindimensionale Bewegung. Z.angew.Math. Mech. 22(1942), S. 206 
bis 215. 
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"Mit Benutzung einer Stabilitätsdefinition von Tre 
schen Lösungen stabil sind und alle anderen stabilen I 
hörigen) periodischen Lösungen konvergieren; die periodi 
Gleichgewichtslagen des schwingenden Systems. Besonder 
schen Lösungen mit möglichst kleiner Maximalamplitude, eine Betrachtung der durch ( 
in Bereiche unterteilten Ebene der Steuerparameter «, ß führt daher zu wichtigen Hinwe 
die Wahl dieser Steuergrößen®). Die Untersuchungen sollen noch fortgeführt und insbesondere auf 


Grundzüge eines Leitfadens zur praktischen Berechnung 
von Kreiszylinderschalen 
Von M. J. De Schwarz!) in Rom 
Für eine Kreiszylinderschale von gleichmäßiger Dicke seien Spannungen, Momente (pro 


Längeneinheit) und die Deformationen der Mittelfläche unter Eigengewicht und geeigneten ’ Se 


Randbedingungen zu bestimmen. 
Zu berechnen sind also Größen der Form 


a 


sin um ((<z<L) 
a) = 2.09) cos PTR VE Are ‚ 
L SYIY) > “ 


d.h. in erster Linie @,(y) und @,(y), wofür im Folgenden kurz @ gesagt sei. 

Es ist @ =@® +@’, wobei im Rahmen einer vereinfachten Theorie @° die sog. Membran- 
lösung und @’ eine zusätzliche Lösung ist, die gestattet die an den Längsrändern vorgeschriebenen 
Bedingungen zu erfüllen. 

Trotz zahlreicher Vereinfachungsvorschläge?) bleibt die praktische Berechnung der @' noch 
ziemlich umständlich, so daß der Gedanke gerechtfertigt erscheint, eine größere Anzahl auf 
typische Fälle bezüglicher numerischer Resultate mit den nötigen Gebrauchsanweisungen in 
einem Leitfaden zusammenfassen, welcher ih Istituto Nazionale per le Applicazioni del Calcolo 
in Rom unter Beratung durch Prof. Ing. G. Krallin Vorbereitung ist. 

Aus dem Gedankengang dieses Leitfadens sei erwähnt, daß die Typeneinteilung nach Aas- 
Jakobsen?) beibehalten ist. Alle Resultate beziehen sich auf gewisse Durchschnittstypen von 
Quertonnen (I) und Längstonnen (II), bzw. auf die Typen ],, I,, II,, II,, je nachdem ob es sich 
um die Berechnung von @, (y) oder @, (y) handelt. Durch die Typen wird das gegenseitige Ver- 
hältnis von Länge, Radius und Dicke der Schalen festgelegt, nicht aber ihr Öffnungswinkel %,. 

Dieser geht in @’, das für genügend großes %, stets auf die Form 


e= 7a [cos ßy(c4A + c5B) + sin By(— cp A + c4B)] 
+? [eos Yyp(coO + 6pD) + siny y(— cp + coD)] 
gebracht werden kann, nur vermittels der A, B, C, D ein, die durch die Bedingungen am Längs- 
randy = (0) zu bestimmen sind und darin bestehen, daß vier von den Größen @, bzw. vier Linear- 
kombinationen derselben gleich Null gesetzt werden. 
Für weniger großes y, ist G’ durch Zusatzglieder ergänzt zu denken, die y, enthalten. 
Durch passende Definition modifizierter Koeffizienten € ,, ©, ©, Cp ist aber — für ein 
symmetrisches Problem — immer zu erreichen, daß am Rand 


@ c„A +czB +c,0 +cpD 
ist, was das Aufstellen der Randbedingungen erleichtert. 


In dem Leitfaden sind für sämtliche ‚Typen‘ und alle interessierenden G@ — die nach dem 
Verfahren von Aas-J akobsen berechneten — 4, CB, Ca; Cp tabelliert. Tabelliert sind 


ferner außer der Membranlösung und den Faktoren e"Ycosß y --- efvsinyy auch die A, B 
C,.D, d.h. die Lösungen von je einem System von vier linearen Gleichungen, für folgende Rand- 


3)E.Trefftz, Zu den Grundlagen der Schwingungstheorie. Math. Ann. 95 (1925), S. 307—312 
. „.) Für eine ausführliche Darstellung sei auf die Arbeit vonH. Bilharz I St. R chottlaender 
Periodische Lösungen einer geregelten Bewegung, Archiv der Math. 5 (1954), S.479 ff. verwiesen. j 


ı) Diese Arbeit wurde im Istituto Nazionale per le Applicazioni del Caleoloin Rom hergestellt 
2) Vgl. A. Kuhelj, Transactions Technical Universityin Brno, XVII, 1948, S. 123—156. —W. Zerna 
Ingenieur-Archiv XX, 1952, S. 357—362. — J.M oe, Assoe. Internationale, Ponts et Charpentes, 1953 
®)A.Aas-Jakobsen,Der Bauingenieur 20 (1939), S. 394—405. 4 
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nd 50°, en a Tabellen m ntlichd edscrte ind igurn! für 
en Größen @ längs der ganzen Schale gegeben. 


Längsrichtung nahezu proportional verläuft. Bemerkenswert sind ferner Beispiele von zwei 


Schalen mit verschiedenem Öffnungswinkel, die am Längsrand zusammengefügt sind, und von 


einer Schale mit y, = 50° und Randglied, wobei klar ersichtlich wird, wie Schale und Randglied 
sich gegenseitig, hinaufziehen. Die G, (y) sind meist klein uler den G, (y), doch kommen 
für den Typ II auch Fälle von en a vor. 


Der Potentialbegriff i in Ger Motorrechnehe und seine Anwendung 
in der Theorie dünner Stäbe 


Von F. Selig in Wien 


Eine er Kraft ist definiert durch die Forderung, daß ihre Arbeit bei Er mit den. 


Bindungen verträglichen, endlichen Verschiebung nur von der Anfangs- und Endlage ihres 
Trägers abhängt. Die Arbeit muß dann eine eindeutige Funktion der Lagekoordinaten 4; sein, 
aus der durch geeignete Gradientenbildung der gegebene Kraftmotor wieder ableitbar sein muß. 

Bezeichnet man den Kraftmotor mit 8 = (8;M,), (8 bedeutet den Kraftvektor und M, 
den auf einen beliebigen Bezugspunkt O0 bezogenen Momentenvektor) und den Geschwindigkeits- 
motor mit 6 = (®;%,) -(& B; = RX 8), (W bedeutet den Winkelgeschwindigkeits- 
vektor, 3, die Bezugspunktsgeschwindigkeit, N den Ortsvektor von 0 zu einem beliebigen körper- 
festen Punkt 0’) dann ist das Differential der Arbeit gegeben durch das Produkt 


IA SEO ae di). 


eher von der formalen Delle der Winkelgeschwindigkeit als Zeitdifferential von 
Quasikoordinaten x; 


dr ; - 
Eee ER Pr ea es de ren (2) 
di 1 j=1. 
läßt sich auch der Geschwindigkeitsmotor anschreiben als 
6 nz mit de (die RER Fe (3). 


Die Koeffizienten des nichtvollständigen Differentials dx; sind Funktionen der. Lagekoordinaten 
Ay Ay As. 
Aus (1) folgt somit 
AZ HRS HaN I WMidnr. een Sr A) 


oder ausführlich geschrieben 
S | 
= NINE NUR re ne een aD) 
il 
Durch Untersuchung der Integrabilitätsbedingungen kann aus (5) stets die Konservativität über- 


prüft werden. 

Zu der der Motorrechnung angepaßten Gradientanbildung gelangt man im Anschluß an (4) 
einfach dadurch, daß man das negative Differential eines Potentials U an Stelle des Arbeits- 
differentials einführt, woraus sich sofort 
ou © 2 
ze, N EN 2 (6 
OR % 0 (6) 


ergibt. Um anzudeuten, daß es sich um eine symbolische Differentiation nach einer Quasikoor- 
dinate handelt, wurde in (6) die Klammer Er Man hat darunter zu verstehen 


u) = > Nuz7, . mit (N)=(Ay. 


# 
Unter den numerischen Resultaten besonders zu erwähnen ist, daß am Rand it 4% 
i positive, zum Teil sehr starke Längsspannungen auftreten, mit denen eine Dilatation in der 
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Zusammenfassend erhält man 

oU [U oU dU 7 

2%, +) Su 

Durch Vergleich von (3) und (7) erkennt man, daß durch Vertauschung von dr und 9/0 sowie 
AN und (9/07) die Lagenänderung d& in die Lagenableitung d/d& übergeht. 

Als Beispiele dienen gewisse Belastungsfälle aus der Theorie dünner Stäbe. An Stelle der 
Komponenten der Winkelgeschwindigkeit treten die beiden Krümmungen und die Torsion, an 
Stelle von dt treten die Verschiebungen, an Stelle der Zeit £ tritt die Bogenlänge s der Zentral- 
linie. Durch Einführung vonRodriguesschen Koordinaten 

j 2 


wurde der Fall näher untersucht, daß zwischen dem Momentenvektor My (0° Endpunkt der Zen- 
trallinie) und dem Tangentenvektor e - e5 5 1) ein linearer Zusammenhang:besteht. 


Bei konkret gestellten konservativen Problemen wird bei Verwendung des Dirichlet- 
schen Prinzips zur Untersuchung der Stabilität die explieite Berechnung des Potentials der Kraft- 
und Momentenbelastung nicht notwendig, da die nach dem Variationsparameter e abgeleitete 
Energie der Belastung wieder durch die Lagenableitung und die Lagenvariation ausgedrückt 
werden kann 

dU dU d& 
de, d& de’ 
Somit treten in den Randgliedern wieder nur die Bestimmungsstücke des Kraftmotors auf. Ähn- 
liches gilt auch für die zweite Variation der Formänderungsarbeit. Hierauf soll später ausführ- 
licher eingegangen werden. 


Optimale Schwingunsstilgung durch Fliehkraftpendel 
Von A. Slıbar in Wien 


Fliehkraftpendel werden nach den einführenden Arbeiten verschiedener Autoren zur Tilgung 
gefährlicher Drehschwingungen und zur Erzielung konstanter Drehgeschwindigkeit der Kurbel- 
welleim modernen Motorenbau häufig verwendet. 

Betrachtet man ein System nach Bild 1, so liefert die linearisierte Rechnung bekanntlich 
die als erste Abstimmbedingung bezeichnete Beziehung 


L 
D=Tan: 0. We. ee VER 


wenn 2 die Kreisfrequenz der erregenden Momente, &, die mittlere Winkelgeschwindigkeit der 
rotierenden Masse ist und @ =» w, gilt. Eine Aussage über die an der Trägermasse vom Träg- 


Abb. 1. Schema des als mathematisches Pendel aufgefaßten Tilgers Abb. 2. Schema des als physikalisches Pendel 
aufgefaßten Tilgers 


heitsmoment J günstigst anzubringende Pendelmasse m folgt aus der linearisierten Behandlung 
nicht. Schreibt man die Grundgleichungen für den als mathematisches Pendel !) betrachteten 
Tilger an und führt die Abkürzungen 

Eee mL AM 


Tor leiten EEE BED) 


EN o Ä > a a 8 (a; =ay vs 49 2a.) gas ir ; 
BE RUN ns OR ET LE NER NEN EB] re RL EN 
en: & z S ; Narr Be v o 4 (q, — 2A ı)* RK: Tr y2 41) K cos vT } 


dc re [ — (2 »)® Cı] % na 

f Teen) a] (a, — »* FR 
ee —-Bral-I; <a +) ara |. 

| Teen BRernee Sn a 


Die Rechtfertigung der Streichung der Glieder mit höheren als der een Potenz in K erweist 
sofort eine Größenabschätzung dieses Parameters nach ausgeführten Maschinen. Darnach 
ergibt sich K je nach Tangentialdruckdiagramm, Kurbelwellenkonstruktion und Drehzahl der 
_ Maschine in der Größenordnung von 105 —10-3. Das Glied 1. Größenordnung in K verschwindet 


cos2vr 


br “ mit Einhaltung von (1). Der für die Amplitude der bezogenen Winkelgeschwindigkeitsschwan- ; 
3 kung entscheidende Term 2. Größenordnung in K läßt sich als Funktion des Parameters O zu 
F einem Minimum machen. Man erhält dann für die optimale Pendelmasse den Ausdruck 
F. | An DRS 
= 5 Got Fr By; E;3 3 (v2 +1): . (5). 
Ferner zeigt sich, daß die jeweils doppelten Werte der optimalen Pendelmassen zu Resonanzen 
führen. 


Wenn die nicht linearisierte Behandlung von Pendelschwingungstilgern auch Aussagen über 
die zur Erzielung geringster Winkelgeschwindigkeitsschwankung günstigst anzubringende 
Pendeimasse liefert, so legt die Betrachtung der ausgeführten Tilgerformen doch deren Behand- 
lung als physikalische Pendel nahe ?). Außerdem ist zu beachten, daß, während das mathe- 
matische Pendel durch die beiden Größen l und m festgelegt ist, bei der Behandlung als physi- 
kalisches Pendel der Trägheitsradius 3; des Pendelkörpers um seinen Schwerpunkt als neue : 
Kenngröße auftritt. Dieser dritte Parameter gewährt aber zur Erzielung optimaler Verhältnisse 
zweifellos wesentlich größere Freizügigkeit in der Konstruktion. 

Legen wir eine Anordnung nach Bild 2 mit den dort angegebenen Bezeichnungen zugrunde 
und verwenden die Abkürzungen 

s? +33 a en, m. L? AM “ 

LT. ir, 12% 7 =@G, ae &t=Ten anby 
so erhalten wir nach verschiedenen Umformungen eine Lösung analog Gl. (4), wenn die Größen 
y- a;(A, ©) durch die Größen a} (A*, ©, o) und die Größen c; (A) durch die Werte c; (4*) ersetzen. 

Bei Einhaltung einer, der ersten Abstimmbedingung des mathematischen Pendels entsprechenden 
} Beziehung 3 
E s? +25 
3 er SD BET TER N 
verbleibt als ausschlaggebendes Glied für die Amplitude der bezogenen Winkelgeschwindigkeits- 
| schwankung Aw/w, der Koeffizient des mit K® multiplizierten Gliedes 


2 ERTL BREREE el 


290.0 +12 |3c Bee nel 


y? 
der sich nach einfacher Umformung für alle Werte von » darstellen läßt. 


1)K.DesoyerundA.Slibar, Zur Berechnung von Pendel-Schwingungstilgern. Ing.-Arch. 21. Band, 
1953, S. 208. 
2) A. Slibar und K. Desoyer, Zur Erzielung optimaler Wirkung bei Föndel- Schwingungstilgern. 
Ing.-Arch. 22. Band, 1954, S. 36. 
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Pa» Py. Wenn in diese Vollscheibe ei 
in Richtung der Hauptspannungen 
‚ursprüngliche gleichmäßige 
 vorgerufene Störfunktion en 


der elliptischen Lochkontur durch einfache Kräfte, deren 
Spannungen mit 1/r vom Lastangriffspunkt abklingen, 
läßt sich das Ziel nicht erreichen. Die erforderlichen 
Grundtypen der Belastung ergeben sich wie folgt: 


2. achsensymmetrische Druckpole 


ws 2 antimetrische Pole 
z. —: B: + fürk— co, Ax—0 
un er 
\ Bildi Bild 2 


Nähern sich zwei entgegengerichtete Kräfte auf ihrer Wirkungslinie und nehmen dabei in 
einer Weise zu, daß sich bei unendlich kleinem Abstand noch ein endlicher resultierender Wert 
ergibt, dann erhält man eine gerichtete Belastung deren Spannungen mit 1/r® vom Angriffspunkt 
abklingen. Durch Überlagerung der gleichen Belastung, jedoch um 90° verdreht, resultiert der 
achsensymmetrische Druck- (oder Zug-)Pol, dessen Spannungsfunktion F, lautet?): 


M=-=0,.lhr: 
Hat die um 90° verdreht überlagerte Belastung entgegengesetztes Vorzeichens, dann resultiert ein 
antimetrischer Pol F,;: 
I, =0,008520, 


Nähern sich zwei achsensymmetrische Pole F, entgegengesetzten Vorzeichen in der oben für zwei 


Kräfte beschriebenen Weise, dann folgt die Spannungsfunktion F, der Belastung deren Span- 
nungen mit 1/r? abklingen: | 


cos 
F,=—(, 


Es läßt sich nun nachweisen, daß der Spannungszustand, der durch das elliptische Loch erzeugt 


wird, ebenso entsteht, durch eine Belastung der Vollscheibe längs der Verbindungslinie der 
Brennpunkte des gedachten elliptischen Loches, wie folgt (s. Bild 2): 


1. Belastung durch achsensymmetrische Druckpole F, der Intensität 
b2 a? 
=(p +9) Vor du. 


‚L.Föppl,,,Drang und Zwang“ 3. Bd., S. 124 ff.. R. Oldenbourg, München 1947. 
pl,,,‚Drang und Zwang“ Bd. I, 3. Aufl. (1941) S. 275 ff. Verlag R. Oldenbourg, München. 
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3. Belastung durch Pole F, der Intensität 


Ban (m Zr) WR du. | a: 


- 
- 


Br Durch Integration und Berechnung der Spannungen wurde im Vergleich mit der bekannten Lö- | 
sung für das elliptische Loch die Richtigkeit dieser Ersatzbelastung der Vollscheibe nachgewiesen. 
h "Beachtung verdient der Sonderfall, in dem sich p,/p, wie a/b verhält. In diesem Fall ver- 
schwinden 0, und O, und es verbleibt nur die Belastung durch achsensymmetrische Druckpole. 


z K Für das Kreisloch vom Radius a lautet die Spannungsfunktion der Störung, wenn p, =0 

2 | LEN Px a? . a? } 3 
S$ 7 =Pe (mr + 0029 — 200529). | 
32 Fa In diesem Ausdruck ist keine dem Pol F, entsprechende Spannungsfunktion, deren Spannungen 
Br mit 2 abklingen, enthalten, dafür klingen die Spannungen des dritten Gliedes in der Klammer 


Pa’ gm 00529 = k, mit 1/r? ab. 


Male N 


DerPolF, entsteht aus Polen F, entgegengesetzten Vorzeichens bei dem oben beschriebenen 
Grenzübergang. e 
3 Diese Untersuchung stellt einen Anfang dar; es wäre die weitere Aufgabe Zusammenhänge 
e- aufzuspüren, die einer Verallgemeinerung den Weg zeigen, um damit noch unbekannte Lösungen 
zu finden. i 


Die Plattentheorie als asymptotisches Integrationsproblem 
der Elastomechanik 


Von E. Weinel in Jena 


Die Kirchhoffsche Plattentheorie, auch in der vonE.Reißner!) undM.Schä- 
f er?) erweiterten Gestalt, ist theoretisch immer noch unbefriedigend, weil in ihr mit Spannungs- 
verteilungen und Verzerrungszuständen gearbeitet wird, die mit den Verträglichkeitsgleichungen 
im Widerspruch stehen und somit in einem ideal-elastischen Körper eigentlich gar nicht möglich 
sind. Die unbestreitbare Bedeutung der Reißner-Schäferschen Entwicklung besteht 
darin, daß drei Randbedingungen erfüllt werden können, wie es der Natur des elastizitäts- 
theoretischen Problems entspricht. Die strengen Lösungen des Plattenproblems vonMichelPP), 
Love®), Neuber?) genügen zwar den Grundgleichungen der Elastizitätstheorie, sind aber so f 
speziell, daß nurz wei Randbedingungen vorgeschrieben werden können. 
j Das im Vortrag entwickelte asymptotische Integrationsverfahren mit Hilfe des Euler- 
F Maclaurinschen Entwicklungssatzes hat sich bei näherem Zusehen als ein Weg zur stren- 
gen Lösung des dreidimensionalen Plattenproblems zu erkennen gegeben. Aus den im all- 
gemeinen divergenten Entwicklungen nach Bernoullischen Polynomen kann leicht die 
allgemeine und konvergente Lösung des dreidimensionalen Plattenproblems gewonnen werden. 

Als ein erstes einfaches Resultat ergibt sich unmittelbar eine Einsicht in den Zusammenhang 
zwischen den Lösungen von Michell-Love-Neuber einerseits und dem Vorgehen von 
Reißner-Schäfer andrerseits. 

Ist nämlich das Verschiebungssystem 


Ur, 92), Vol® Y 2), Welt, Y, 2) 


an 


U a a Ed A A A al a Ed Hm: 


1) E. Reissner, On bending of elastic Plates. Quart. Math. 7. (1949) S. 223—228. 

2) M. Schäfer, Über eine Verfeinerung der klassischen Theorie dünner schwach gebogener Platten. 
Z. angew. Math. Mech. 32. (1952) S. 161—171. 

3) Michell, On the direct Determination of Stress in a elastie solid. London Math. Soc. Proc. 31. (1899) 
S. 100. 

4) Love, Lehrbuch der Elastizität. (1906) $ 299 und folgende. 

5) Neuber, Über das Kerbproblem in der Plattentheorie. Z. angew. Math. Mech. 20. (1940) S. 199 
bis 209. 
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eine Lösung des Plattenbiegungsproblems nach Michell, Love oder Neuber, so erhält 
man in der Gestalt 


u— ul, y,2) + 59@ y) sin 5 | 
0 nz ee BET 
= u(2, %, 2) — 7, y) sin 57 
w= wy(X, %, 2) ) 
eine neue Lösung der Grundgleichungen, wenn y(z, y) der Differentialgleichung 
RYPy—la2?y=0.....-. EU EN 


enügt. 
i ED ist aber gerade die beiSchäferfür die Momentenfunktion y angegebene Differential- 
gleichung, wenn man von dem numerisch unwesentlichen Unterschied zwischen (2/4) und (5/2) 
absieht. 

Mit (1) lassen sich dann drei Randbedingungen erfüllen und man überzeugt sich leicht, daß 
die Näherungstheorie vonReißnerund Schäfer die aus (1) gewonnenen elastisch-strengen 
Ergebnisse numerisch gut wiedergibt. 

Eine ausführliche Darstellung des dreidimensionalen Plattenproblems mit allgemeinen 
Randbedingungen, insbesondere auch die Behandlung von Schwingungsproblemen folgt in 
einer späteren Arbeit. 


Referate folgender Vorträge sind nicht eingegangen: 


N.Forbat, Mons: Vibrations de flexion d’un arbrecylindrique en rotation rapide, 


K. Magnus, Frankfurt/Main: Die Stabilität des kräftefreien unsymmetrischen 
Kreisels in kardanischer Lagerung. (Eine ausführliche Arbeit über diesen Gegenstand wird 
in einem der nächsten Hefte dieser Zeitschrift erscheinen.) i 


H. Schwieger, Halle/Saale und G. Haberland, Wittenberg: Die voll- 
ständige Bestimmung des Biegespannungszustandes elastischer quadratischer Platten 
mittels des optischen Zweischichtverfahrens und einer neuen Integrationstheorie. 


D. Strömungslehre 


Zur Differentialgeometrie der ebenen Strömung von Gasen 
Von A. Basch in Wien 


Bei der stationären ebenen Strömung von Gasen treten an Stelle der für unzusammen- 
drückbare Flüssigkeiten geltenden Cauchy-Riemannschen Gleichungen die Beziehungen 


po Typ Leo BR ee 
0% 0% 


in denen ® das Geschwindigkeitspotential, Y die Stromfunktion, o die mit dem Ort veränderliche 
Dichte, go, deren Wert in dem als Koordinatenursprung gewählten betrachteten Feldpunkt O be- 
deutet. DieBernoulische Gleichung 


2 


112 1 
rad — = — —- grad Pe N . RRRSBERER 
grad 7 nr 2) 


gibt einen Zusammenhang der Veränderlichkeiten des Betrages W der Strömungsgeschwindig- 
keit % und des Druckes p. Wird die Zustandsgleichung der Gase 
p=RoT wz.n... Se ee), 


(R Gaskonstante, 7 absolute Temperatur) herangezogen und adiabatischer Charakter der Zu- 
standsänderungen vorausgesetzt, etwa in der Form 


aT d 
Tun 


!)Vgl.R.Sauer, Theoretische Einführungin die Gasdynamik. S. 20. Berlin 1943. 


La alu Do la 2 ur 


ur a 


” 


Z. angew. Math. Mech. h el 
Ba. 34 Nr. 8/9 Aug./Sept. 1954 D. Strömungslehre 333 


(x Verhältnis der spezifischen Wärme bei konstantem Druck und bei konstantem Volumen), so 


besteht zwischen den Änderungen des Geschwindigkeitsbetrages und der Dichte der Zusammen- 
hang 


Narr IE: z 
gradp = P grad u ee are (5), 
wobei 
* Po 
Eu N er | A ae TE er 6 
= (6) 


die adiabatische Schallgeschwindigkeit im Koordinatenursprung 0 bedeutet. . 
Geschwindigkeitspotential und Stromfunktion seien in der Umgebung des betrachteten 
Feldpunktes O (0, 0) bis auf Glieder zweiter Ordnung entwickelt 


1 
® =6d, +9,27 + 0d,y - 5 @u® +20.2Y + ©?) 


1 
„=-W%+9Mr+ P,y eo; Fu +2 P.2Yy+ Pay?) 


und das Koordinatensystem derart orientiert, daß die x-Achse Tangente an die Niveaulinie, die 
y-Achse Tangente an die Stromlinie ist. Dann ergibt sich 


(D,,)o == o, = 0, (D,)o == ®, —— W, a OEL AO De (8) 
und die Komponenten der Strömungsgeschwindigkeit 8 sind 
u=®, =D,T + DY v= Pd, =Dd+ 9:7 + Day Be Rn (9), 
somit bis auf Glieder erster Ordnung 
W"=w+v”=03+20, 0:2 +20,0,Y u Or ec (10) 


und die Komponenten des auf der rechten Seite von (5) auftretenden Gradienten 


W°® o /W° 
2 — ®, D;s = W, Ds 9 5) = DB, DB. = W, D; A een (11) 
und damit zufolge (5) ! 
ae ei 
0 Co 
Die hieraus folgenden Werte von o, und o, in (1) eingesetzt ergeben die Beziehungen 
Y=--9%,=-—MW %=-9d,=-0, 
2 ra il EN 7 (13). 
Yı re 0.1 = =) ’ Pie = —d, (1 ar en ’ Pa=b.| 
2: 69 


Approximiert man Niveaulinie und Stromlinie durch den Feldpunkt 0, (= 0, Y = Y,) 
durch die Parabeln 


1 
= er a Du N a N: (14), 


er 
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wobei k, und ky die ersten Krümmungen von Niveaulinie und Stromlinie bedeuten, so gibt der 


Vergleich dieser Gleichung mit (7) 


® a ee 
ee: Fee ge a | 


und hieraus die von jeder Einheitenwahl unabhängige Vektorgröße 


= 1 ’ : 
grad In W = — grad n m (Disi+ Dj) 
1 ee ko... 3:23 up er al ee 
2 Pat pl TE Mar (1) BIN. i 
er. c? ® ca 


Hierbei bedeuten i und j die Grundvektoren, r, und r„ die Krümmungsradien der Niveaulinie 
und der Stromlinie, £ und n die x- und y-Komponente des durch (16) bestimmten Gradienten, die 
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Bild 2 


zu den Krümmungsradien der beiden Feldlinien in den Beziehungen stehen 


1 
AN NIS ee ee 


Wir ziehen nun in den Bildern 1 und 2 die Verbindungsgerade der Krümmungsmittelpunkte 
von Stromlinie Öy und Niveaulinie 0, (Ce sär 0o-Gerade)?). Der Endpunkt @ des von O aus ge- 


zogenen Vektors OQ = grad In W ist der Pol der Cesär o-Geraden und zwar im Fall von Unter- 
schallgeschwindigkeit (W, < c,) bezüglich der Ellipse 


2 
+1 -)y=1. N RER 
0 
(Bild 1) im Fall von Überschallgeschwindigkeit (W, > c,) bezüglich der Hyperbel 
2 
= (R-1)y - Se Ze ae nz 
Co 


(Bild 2). In einfachem Zusammenhang mit den durch (16) gegebenen Gradienten stehen die 
Gradienten für Dichte und Druck. 


?) Bezüglich ihrer Bedeutung in der Differentialgeometrie ebener Kurvenscharen vgl. E.Cesäro ‚ Vor- 
lesungen über natürliche Geometrie. Deutsche Ausgabe von G.Kowalewski. $.148. Leipzig 1901. —R.v. 
Lilienthal, Vorlesungen über Differentialgeometrie. $. 168. Leipzig 1908. 

®) Beiunzusammendrückbaren Flüssigkeiten tritt an Stelle dieser Ellipse der Einheitskreis. Vgl. A. Basch, 
Mh.Math. Phys. 41, 300—321(1934). — Proc. Second InternationalConference of SoilMecanics and Foundation 
Engineering, Rotterdam, 1948, Bd. 5, Sect. 10 u. 5, 8. 280-285. 
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en von der linearen Potentialgleichung für das unsymmetrische pl Fi in = EL 


© Zilinderkoosoinaten 
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Shen Entlang einer nach außen een Charakteristik Hi = const) gilt das positive und 


längs der nach innen gerichteten das negative Vorzeichen. 


Die neuen Variablen fru und vo + = ändern sich in stetig veränderlichen Strömungs-. 


feldern relativ wenig, wodurch es möglich ne mit verhältnismäßig großen Schritten zu arbeiten. 
Für den Grenzfall v—+ &- = const erhält man mit dem Potentialansatz #=r-+ —Z 2 ee 
schlanken Körpern der er = 7, direkt die bekannte erste Näherung von Tsien: 


dh 
uh) = 27,- 


Die Berechnung erfolgt leicht in Form einer Konstruktion in der Geschwindigkeitsebene 
Di = ß rw für feldinnere Punkte und erfordert etwas numerische Rechnung am Rand. 


Es ist vorgesehen, daß diese Arbeit in der ZFW erscheinen soll. 


Erweiterte Theorie der Grundwasserströmung 
Von @. Heinrich in Wien !) 


Für eine einwandfreie Darstellung der Grundwasserströmung ist es wesentlich, zwischen 
Mikro- und Makroströmung zu unterscheiden. Unter Mikroströmung sei die tatsächlich auf- 
tretendeBewegung der Flüssigkeit in dem unregelmäßigen Porenraum des Festkörpers verstanden. 
Sie ist nur insoweit von Interesse, als sie durch geeignete Mittelwertsbildungen durch eine Makro- 
strömung ersetzt werden kann. 

Es werden folgende Annahmen getroffen: 1) die Körner des Festkörpers berühren einander 
punktförmig, 2) der Festkörper befindet sich dauernd in Ruhe, 3) in der Mikroströmung erfährt 
die Flüssigkeit nur an den Körneroberflächen eine Hautreibung (im Innern heben sich die Schub- 
spannungen im Mittel auf), 4) im Innern der Strömung soll die Flüssigkeit den Porenraum des 
Festkörpers vollständig erfüllen, 5) die Filterwirkung des Festkörpers sei so groß, daß bereits bei 
sehr kleinen Strömungsgeschwindigkeiten die Trägheitskräfte in der Flüssigkeit verschwindend 
klein werden gegen die übrigen Kräfte, 6) Kapillarkräfte und Verdunstung werden nicht be- 
rücksichtigt. 

Aus dem System Festkörper-Flüssigkeit werde ein beliebiges Volumen herausgeschnitten, 


“so daß kein Korn von der Schnittfläche zerteilt wird. Wir teilen die wirkenden Kräfte in Ober- 


flächen-, Volums-, und Quasivolumskräfte ein. Die letzteren erhält man, indem man die in der 
Mikroströmung zwischen Körnern und der Flüssigkeit wirkenden Kräfte in der Makroströmung 
auf das Volumselement gleichmäßig verteilt. Man kann sie in einen Reibungs- und einen Nor- 
maldruck-Anteil zerlegen. Für den Flüssigkeitsanteil des herausgeschnittenen Volumens gilt bei 
quasistationärer Flüssigkeitsbewegung die Vektorgleichung 


— [pa [rnpvar + [avar + [a-mrpar =0. ER RSER Ey: 


1) Der vorliegende Bericht bezieht sich auf eine Arbeit, die gemeinsam mit K.Desoyer, Wien, durchge- 
führt wurde. 


ern sich mit Fa M’,—1,F,= uundF, = ventlang der Charakteristiken die Proportiona- 
_ litätsbeziehungen 


- £ ” 
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Hierin bedeuten: n das auf die Volumseinheit bezogene Porenvolumen, p den Druck in der 


Flüssigkeit, do einen nach außen positiv gezählter Vektor des Oberflächenelementes, yy das spez. 
Gewicht der Flüssigkeit, U das Potential der Volumskraft (Schwerkraft), dV das Volumselement, 
N bzw. (1—n) Yp den Reibungs- bzw. Normaldruckanteil der Quasivolumskraft. Da die 
obige Beziehung für ein beliebiges Volumen gelten muß, folgt — bei Anwendung des Gau ßB- 
schen Satzes auf das Oberflächenintegral- — aus dem Verschwinden des Integranden die Be- 
ziehung 


Dale) Fa a 

Analog erhält ER für den Festkörperanteil: 
[&x: o- | Pr vuar— | mar | a—mrrar=0.. . 8). 
©x bedeutet den eisch Spannungsaffinor der Korn-zu-Korn-Kräfte, die man in der 


Makrobetrachtung gleichmäßig über das Oberflächenelement gemittelt denkt, yz ist das spez. 
Gew. des Festkörpermaterials. Analog wie oben folgt aus (3) bei Verwendung von (2): 


7-& = u z a er > (4). 


Die Gl. (4) ist für die Beurteilung der Rutschgefahr von Dämmen von Bedeutung. 

Wir definieren einen Filtergeschwindigkeitsvektor b so, daß dg = d : do das durch do pro 
Zeiteinheit hindurchtretende Flüssigkeitsvolumen bedeutet. Unter Beschränkung auf jene Fälle, 
bei denen ein linearer Zusammenhang zwischen v und Aybesteht, setzen wir, um auch anisotrope 
Festkörper (Schichtungen) zu erfassen: 


NRreslat uu yIN re 
worin X ein beliebiger i. a. unsymmetrischer, ortsabhängiger Affinor ist. Aus (5) und (2) folgt: 
E 
= er + 2) 8:92 mit =) Se 
BE. N Yw 


X kann als Durchlässigkeitsaffinor bezeichnet werden, h ist die Standrohrspiegelhöhe. Für 
sotropen Festkörper erhält man daraus das Darcy sche Gesetz 


be—krh. 
Wegen der Quellfreiheit der Filtergeschwindigkeit (7 - d = 0) ergibt sich aus (6) die Beziehung: 
ve r)=- VD: WM FE 
Für homogenen, aber anisotropen Festkörper geht (7) über in: 
8.4, P)h = ee 


Bei symmetrischem St kann man in letzterem Fall $ auf das generelle Hauptachsensystem &, n,& 
transformieren und erhält aus (7a) die Differentialgleichung: 


.0®h 0°h o®h 
gr 0&® im Kun On? I ker 0&° 


al 


(7b) kann durch eine volumsgetreue Affinverzerrungin die Laplace sche Differentialgleichung 
übergeführt werden. Damit ist der anisotrope Fall auf den isotropen zurückgeführt. Die Diffe- 
rentialgleichung (7), (7a) oder (7b) für die Größe Ah ist nun für die bekannten Arten von Rand- 
bedingungen, wie sie bei Grundwasserströmungen vorkommen, zu lösen. 

Um auch instationäre Probleme zu erfassen, benötigt man noch eine Differentialgleichung 
für die Bewegung des Grundwasserspiegels. Sie kann durch eine volumetrische Bilanz ge- 
funden werden. Wir betrachten einen prismenförmigen Körper von infinitesimalem Quer- 
schnitt und bringen zum Ausdruck, daß der Überschuß der durch die Seitenflächen eintretenden 
gegenüber der austretenden Flüssigkeitsmenge zu einer äquivalenten Vergrößerung der von 


?).. bedeutet das doppelt skalare Produkt, das Zeichen, das dyadische Produkt. 


, bzw. 2, die Höhe der undurchlässiger 


v 


de 


2 y% ()de Pe 


02 Be; 
= n(z, Y, %,) eu dx dy. 


Daraus folgt, wenn man berücksichtigt, daß die undurchlässige Schicht eine Stromlinienfläche ist: 


920, 1 092, (2, y,t) 


W 02, (2, y Ü) | (8). 
ot N(T, %, 2) 0% oYy. 1 


PB (%, Y%, 2, 1) — v,(%, %, 2a l) —U,(8, Y,2, 2) 


Im stationären Fall wird der rechte Klammerausdruck Null und dies bedeutet, daß die Spiegel- 


* 


fläche ebenfalls zur Stromlinienfläche geworden ist. In der Veränderlichen h geschrieben lautet (8) 
für homogenen Festkörper: 0 


(1 Re, Abt 
ot N (x, Y; 2.) r 05 Fe 0£ 0% a on 0Y längs der Spiegelfläche ; 


Mit Hilfe dieser Gleichungen kann man, am vorteilhaftesten mittels der Relaxations- 
methode, instationäre und stationäre Probleme der Grundwasserströmungen untersuchen. Auf 
diese Weise wurden bereits die instationäre Durchströmung eines Dammes und der stationäre Zu- 
fluß zu einem rotationssymmetrischen Brunnen behandelt, wobei sich insbesondere die genaue 
Lage der Hangquelle ergibt, die in vielen Fällen nach den bisherigen Methoden nicht einwandfrei 
ermittelt werden konnte °). 


Überschallströmungen und langsam pendelnde Drehkörper 
Von (. Heinz in Lörach 


Das bekannte Singularitätenverfahren von v. Kärmän,Mooreund Tsien zur Be- 
rechnung der stationären Überschallströmung um schlanke, schwach angestellte Drehkörper 
wurde von Sauer auf die Strömung um langsam und mit geringer Amplitude pendelnde Dreh- 
körper ausgedehnt. Hierbei tritt zu den beiden die axiale und schiefe stationäre Strömung be- 
herrschenden, im allgemeinen numerisch zu lösenden, Integralgleichungen noch eine weitere, die 


- zur Berechnung der nichtstationären Strömung zu lösen ist. 


Die numerische Lösung dieser Integralgleichungen, die schwierig und mit erheblichen 
Rechenfehlern behaftet ist, wurde vom Verf.!) dadurch umgangen, daß sie näherungsweise durch 
lineare gewöhnliche Differentialgleichungen ersetzt wurden. Dieses Verfahren wird auf die nicht- 
stationäre Strömung ausgedehnt. Zu den beiden die stationäre Strömung bestimmenden Diffe- 
rentialgleichungen tritt noch eine dritte, deren Koeffizienten mit der der schiefen stationären 
Strömung übereinstimmen. Die in den Differentialgleichungen auftretenden unbekannten 


- Funktionen sind die A belschen Transformierten der Belegungsfunktionen der die Strömung 


bestimmenden Quellen- und Dipolverteilungen auf der Achse des Drehkörpers. Diese zeigen bei 
etwaigen Unstetigkeiten der Meridiankurve dasselbe Unstetigkeitsverhalten wie die Meridian- 
kurve selbst. Dadurch wird die Behandlung von Drehkörpern, deren Meridiankurven Knicke 
aufweisen, auch gegenüber dem Verfahren von Lighthill wesentlich erleichtert. Da die 
numerische Lösung einer gewöhnlichen Differentialgleichung der numerischen Quadratur einer 
Funktion entspricht, ist das Verfahren sehr genau; der Rechenaufwand wächst andererseits nur 
linear mit der Anzahl der Rechenschritte. Die Koeffizienten der Differentialgleichungen sind 
Funktionen, die weitgehend ein für allemal tabelliert werden können. 


3) Eine ausführliche Darstellung soll demnächst im Ing.-Arch. erscheinen. 
1) Tagung GAMM 1953 (Aachen). 


a stationär ist, d. h. damit das Geschwindigkeitsfeld zeitunabhängig b 


‚ gen ug 
en ne ale zur 7 Zeit 1a. ein vorge: Jene: 
_ durch ein entsprechendes Wirbelfeld erzeugt denken kann, 
Zeit i, das Geschwindigkeitsfeld beschaffen sein, damit 


Wir setzen voraus, daß zur Zeit , die Wirbellinien Wirbelröhre 
bilden, die von doppelpunktfreien, stetig gekrümmten Flächen ‚begrenzt wer 
- weder ins Unendliche erstrecken oder vom topologischen Typus des Torus sind. Die $ roml; 2 
diesen Wirbelröhren seien einfach geschlossene stetig gekrümmte Kurven; etwaige 
(d.h. Punkte mit verschwindender Geschwindigkeit) sollen nur isoliert vorkommen : za 
isolierten Staulinien zusammenschließen. 2 

Wird allgemein die Flüssigkeitsbewegung durch die Gleichungen 


z—= %a,b,6,0, y= yYa,b,cl, 2= :Ha,b,ct) 
beschrieben, [d. h. also, das zur Zeit t, in (a, b, c) gelegene Teilchen befindet sich zur Zeit t im er 
Punkte (x, Y, 2)], so bestehen zwischen den zu £, bzw. £ gehörigen Wirbelkomponenten &,, 270» &o € 
bzw. S n,& ein und desselben Teilchens die Cauchy schen Relationen h1 a 


0 B 
a rn er 


Soll also die durch ein gegebenes Geschwindigkeitsfeld u,(a, b, c), E- b, c), w(a, b, c) festgelegte 
Flüssigkeitsbewegung stationär sein, so muß, da ja die Bahnkurven mit den Stromlinien zusam- 
menfallen und der lokale Wirbelvektor zeitunabhängig ist, auf den in den Wirbelbereichen zur 
Zeit i, gelegenen Stromlinien die Wirbelverteilung den Beziehungen (1) entsprechen. 

Da diese Beziehungen zunächst nur für Bewegungen einen Sinn haben, hat man in unserem 
Falle, d.h. bei einem Stromlinienfeld, darunter folgendes zu verstehen: Wir erklären eine Ab- 
bildung der gegebenen Stromlinienkonfiguration in den Wirbelbereichen auf sich selbst derart, 
daß die Stromlinien in sich übergehen. Bezeichnen etwa (a, b, c) bzw. (x, y, 2) die hierbei einander 
entsprechenden Punkte einer Stromlinie und sind diese durch die Bogenlängen s, bzw. s festgelegt, 
so erscheint die Abbildung in der Gestalt 


ve 9ob,c;9), Y = v(a,b;c}H), 2 (a, b,c039)>) aa. (2); 


Staupunkte sollen in sich selbst übergehen. Ersetzen wir schließlich auf den Stromlinien den 
Parameter s durch 


°) 
=-4h+ [| — (3) 
: ara = 
‚So 
so erhält man die Gleichungen (2) in der Gestalt 
x = .o0(0,5,c;0); %= Ylasb.0,0), 2 .larbee: oO) . (4). 


0% 0% 


Bildet man hieraus die Größen en ° und setzt diese in (1) ein, so bekommt (1) 


es 0 
”0b’ "dc 
jetzt auch für das Stromlinienbild zur Zeit it, einen Sinn, und wir erkennen, da während der 
stationären Bewegung { die Rolle von o übernimmt, daß das Bestehen der im genannten Sinne 
interpretierten Ca uch y schen Relationen auf den Stromlinien eines einer stationären Bewegung 
entsprechenden Stromliniensystems notwendig ist. 

Wir behaupten nun: Das Bestehen der Relationen (l) auf den 
Stromlinien des gegebenen Geschwindigkeitsfeldes ist auch 
hinreichend für eine stationäre Bewegung. 

Für das weitere bieten die bei den eingangs erwähnten Existenzsätzen entwickelten 
Methoden die erforderlichen Grundlagen. Lichtenstein geht dort (vgl. c.!)) nicht von 


RL Vgl.L. Liechtenstein, Grundlagen der Hydromechanik, Berlin 1929, S. 415 if. 
°) Man beachte, daß in den Relationen (2), die man auch als die Gleichungen des Stromliniensystems auf- 
fassen kann, a, b, c nicht voneinander unabhängig sind; eine Stromlinienschar ist bekanntlich zweiparametrig. 


3) oist bei der stationären Bewegung mit der Zeit t, die das Teilchen (a, b, c) braucht, um von s,nach s 
zu gelangen, identisch. 


la dan un 


ger 
in von nn ser die rn Funktionen 
z=x(a,b,c, t)..., das den Druck nicht mehr enthält. Der Beweis der Existenz genau eines 
Lösungssystems erfolgt dann dort für hinreichend kleine it — bo ‚durch esenet konstruierte 
sukzessive Näherungen. 3 
Zum Beweise unserer Behauptung an man das Lic h tensteinsche en | 
und benutzt abweichend als erste Näherung die Verschiebung der Teilchen längs ihrer Strom- 


linie. Unter der Voraussetzung (1) zeigt sich, daß diese Näherung bereits die Lösung darstellt. 
Man erkennt leicht, daß in jedem Raumpunkt dauernd dieselbe Geschwindigkeit herrscht und 


also auch die Betrachtungen „im Großen“, d. h. für beliebiges t — t, gelten. Es ergibt sich weiter, 
daß für Bewegungen, die bezüglich eines gleichförmig rotierenden oder sich parallel zu sich selbst 


verschiebenden Bezugssystems stationär sind, entsprechende Betrachtungen zum Ziel führen. 

Die zunächst recht umständlich erscheinende Bedingung (1) vereinfacht sich wesentlich, 
wenn wir speziell ebene oder rotationssymmetrische Ausgangsfelder (d. h. Wirbelringe) zugrunde 
N Im ersten Fall reduziert sich (1) auf 

Bene, a 
falls man die Vorgänge in der x, y-Ebene betrachtet. 

Im zweiten Fall erhält man 


= m. u a a a Be a 5 


mit der x, z-Ebene als Meridianebene; a bzw. x sind der Anfangs- bzw. Endabstand eines Teil- 
chens von der Symmetrieachse. Da zudem in den genannten Spezialfällen Stromfunktionen 
existieren, erleichtert sich die Konstruktion von Stromliniensystemen, auf denen die Bedingungen 
(5) bzw. (6) gelten. Auf diesem Wege finden Untersuchungen von Lichtenstein über 
Systeme geradliniger Wirbelröhren?®) ihr notwendiges theoretisches Fundament. Die entsprechen- 
den Betrachtungen vonL. über die Existenz gewisser rotationssymmetrischer Wirbelbewegungen?) 
bedürfen allerdings noch einiger Modifikation, da dort die stationäre Bewegung eines Ringes mit 
überall konstantem Wirbelvektor gefordert wird, den es wegen (6) nicht geben kann. 


Über die Anwendung der Impulsmethode bei Messungen 
an Schaufelgittern 


Von N. Scholz in Braunschweig 


Für die Ermittlung des Strömungswiderstandes von Tragflügelprofilen und anderen um- 
strömten Einzelkörpern nimmt seit längerer Zeit die sogenannte Impulsmethode eine hervor- 
ragende Stellung in der Meßtechnik ein. Diese Methode gestattet mit sehr großer Genauigkeit die 
Bestimmung des Profilwiderstandes durch punktweise Ausmessung des hinter dem umströmten 
Körper infolge der Strömungsverluste auftretenden Nachlaufes. Für die Auswertung derartiger 
Nachlaufmessungen sind dieB etz sche und Jonessche Formel in Gebrauch (vgl. [1], S. 462 ff.). 
Da die Ermittlung der Strömungsverluste durch Integration des örtlichen Gesamtdruckverlustes 
über den Nachlauf des Körpers besonders bequem ist, sind diese Formeln auch so umgeschrie- 
ben worden, daß neben einem Integral über die Gesamtdruckverlustdelle als Hauptglied ein 
additives Korrekturglied auftritt, das von der Form der Nachlaufdelle abhängt. Durch Annahme 
einer universellen Form für das Nachlaufprofil ist von A. D. Young [2] dieses Korrekturglied 
universell berechnet worden. 

Bei der Anwendung der Impulsmethode für Nachlaufmessungen an ebenen Schaufelgittern 
sind in vielen bisher bekannt gewordenen Untersuchungen lediglich der Integralwert des Gesamt- 
druckverlustes sowie Mittelwerte des statischen Druckes und des Abströmwinkels benutzt worden. 
Ebenso wie beim Einzelprofil treten hierbei jedoch Korrekturglieder auf, die nicht nur bei 
der Ermittlung des Strömungswiderstandes, sondern auch bei der Berechnung des statischen 
Druckes und des Abströmwinkels der Strömung hinter dem Schaufelgitter eingehen. Die Ab- 
weichungen gegenüber den Mittelwerten der Meßebene können je nach Abstand der Meßebene 
hinter dem Schaufelgitter bis zu 20% betragen. Auch hierbei liegt es nahe, diese Korrekturen 
durch Zugrundelegung einer universellen Form des Nachlaufes ein für alle Mal zu berechnen und 
damit die Auswertearbeit erheblich abzukürzen. 

Betrachten wir eine Kontrollfläche X, (Bild 1), die ein Schaufelprofil eines unendlich langen 
ebenen Schaufelgitters so umschließt, daß in der Eintrittsebene (Ebene 1) und in der Austritts- 


4) Val.1.c. 450 ff. 
5) Vgl.1.c. 443 ff. 
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ebene (Ebene 2) der Strömung eine homogene Strömung herrscht, d.h. alle Zustandsgrößen der 
Strömung nicht mehr von der Koordinate y längs der Gitterfront abhängen, so läßt sich mit Hilfe 
des Impulssatzes sofort die resultierende Schaufelkraft angeben, wenn die Strömung in den 
Ebenen 1 und 2bekannt ist. Da infolge 
des Nachlaufes hinter jedem Schaufel- 


Ebene 2 erst sehr weit hinter dem 
Schaufelgitter erreicht wird, muß die 
Messung in einer Ebene 2, erfolgen, in 
der die Strömungsgrößen, insbesondere 
aber die Geschwindigkeit noch von der 
Koordinate y abhängen. Durch noch- 
maliges Anwenden des Impulssatzes 
auf eine Kontrollfläche K, (Bild 1) 
kann jedoch aus dieser Messung in 
der Ebene 2y der Strömungszustand 
der homogenen Strömung in der 
Ebene 2 ebenfalls ermittelt werden, 
so daß damit die Bestimmung der 
Schaufelkräfte auf eine Ausmessung 
der Strömung in der Meßebene 2% zu- 
rückgeführt ist. 


Kontrollflache K, 


Bild1: Berücksichtigt man bei der Durch- 
Zur Anwendung des Impulssatzes bei Messungen an ebenen Schaufelgittern. führung der Messung die Veränderlich- 
Ebene ®: sehr weit vor dem Gitter homogene E >” 2 e 
Ebene (): sehr weit hinter dem Gitter f Strömung. keit des Strömungswinkels 1 yın der 
Ebene 2y) : Meßebene, inhomogene Strömung. Meßebene 2y, so ergeben sich ziemlich 
Kontrollfläche X,: begrenzt durch die Ebenen @) u. @). umständliche Auswerteformeln. Durch 
Kontrollfläche X,: begrenzt durch die Ebenen (29) u. @ Einführung eines konstanten Mittel- 


wertes für den Strömungswinkel ß3, 
läßt sich jedoch die gesamte Auswertung auf die Berechnung folgender drei Integralwerte 
zurückführen: 


1. mittlerer Gesamtdruckverlust 


Yot 
a (a) 
t PA Qı ' 
2. mittlere statische Druckdifferenz 
I ( ) 
je Pay\Y) — Pı 
f # a a an 
% 
3. Korrekturglied 
ER nn: y) 3 i 
a aut) a ar Wa y\Y 
K an u dy E EI RT RI EEE NN 
Yo Y 
Hierbei bedeutet gq, = 2 w; den Staudruck der Zuströmgeschwindigkeit und g den Gesamt- 


[\ 


druck der Strömung. 


Diese Integrale sind in der Ebene 2y über eine Teilung zu nehmen. Die gesuchten Strö- 
mungsgrößen der Ebene 2 erhält man dann zu: 


91 — 9 
— 2 GR . . . . . . . . + . . . . . . 
N - I 
OP ; 
m Pan Re 
K 
eig ß = eig Pry|1 | —). TE IERTR 0): 


Wird nun für das Geschwindigkeitsprofil des Nachlaufes längs einer Schaufelteilung ein uni- 
verseller Ansatz in Form einer Ga u ßschen Glockenkurve gemacht, dann läßt sich das Korrek- 


profil die homogene Strömung der _ 


r des Geschwintigkeitsprohle, tritt ale, der ee Hr 
kgeführt auf die Größen @ und P) und die Dellen- Ei n 
den maximalen örtlichen Gesamtdruckverlust in der 2% auf. Die fol- 


thält die berechneten Werte des Korrekturgliedes: 


P nach Gl. (2), 


ERS g ee 


1—-P 


. 9 — Iaymin _ 


u Tabollel: : Bereehnete Wertedes Korrekturgliedes Knach Gl. @), @nach Gl. (), 


a ri a Aa 


2 0,2431 
3 0,1902 
4 ‚0,1491 
„5 0,1156 
8 0,0869 
7 0,0616 
2 0,0389 


Die hiermit berechneten Werte der homogenen Abströmung zeigen vorzügliche Überein- 
stimmung mit Auswerteergebnissen unter Berücksichtigung der Veränderlichkeit des Strömungs- 
winkels ß,,in der Meßebene 2y. Die Abweichungen der homogenen Strömung der Ebene 2 gegen- 
über den Mittelwerten der Meßebene 2, sind vor allem beim Gesamtdruck und statischen Druck 
beträchtlich. Die hier mitgeteilte Auswertemethode hat sich bei den im Institut für Strömungs- 
mechanik der Technischen Hochschule Braunschweig durchgeführten Gittermessungen gut be- 


währt (vgl. [3]). 


MB: Schlichting: 


[2] A. D. Young: Note on a Method of Measuring Profile Drag by Means of Integrating Comb, ARC Rep. and 


Mem. No. 2257, (1948). 


[3] H. Schlichting: Ergebnisse und Probleme von Gitteruntersuchungen. Zeitschr. f. Flugwissenschaften 
Bd. 1, S. 109—122, (1953). Vgl.auch: H.Schlichting: Problems and Results of Cascade Investigations. 


0,2353 
0,1832 
0,1425 
- 0,1093 


0,0809 


0,0559 
0,0335 


0,2276 


0,1761 
0,1359 


0,1031 
0,0751 


0,0503 


0,0281 


0,2199 
0,1690 
0,1293 
0,0969 
0,0691 
0,0447- 
0,0227 
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0,2122 
0,1619 
. 0,1227 
0,0907 
0,0632 
0,0390 
0,0173 


0,2044 
0,1548 
0.1161 


0,0845 


0,0573 


0,0334 


0,0119 


Grenzschichttheorie. Verlag Braun, Karlsruhe, (1951). 


Journ. Aeron. Sciences, Vol. 21, p. 163—178, (1954). 


Zur Ermittlung von Unterschallströmungen mit der Transformations- 
methode bei quadratischer Approximation der Adiabate 


Durch Übergang zur Hodographenebene wird die Kontinuitätsgleichung für die Strom- 
funktion einer kompressiblen Unterschallströmung um ein gegebenes Profil nicht wie beiTsch.a- 
pligin und bei Sauer auf dieLaplacesche Potentialgleichung sondern auf die Schwin- 
gungsgleichung transformiert und gezeigt, daß die zugehörigen Druck-Dichte-Beziehungen eine 
quadratische Approximation der Adiabatenkurve gestatten. Ferner wird in der Hodographen- 
ebene die Singularität ermittelt, die die Stromfunktion im Bildpunkt der Grundströmung besitzt. 


Von H. Schubert und E. Schincke in Halle (Saale) 


0,1967 


0,1477 _ 


0,1095 
0,0782 
0,0514 
0,0277 
0, 0065 


0,1812 


0,1335 


0,0963 
0,0658 


0,0396 


0,0165 


Eine ausführliche Darstellung wird in einem späteren Heft dieser Zeitschrift erscheinen. 


Kegelige Überschallströmung in Schallnähe 


Die Überschallströmung an Kreiskegeln ist vielfach berechnet worden. In diesem Zusam- 
menhang braucht nur auf die in allen einschlägigen Lehrbüchern zitierten Arbeiten von G. J. 
Taylor—J. W.Macoll, A. Busemann undauf die Tabellen von Z. Kopal hinge- 
wiesen werden. Dabei handelt es sich stets um die numerische Integration einer gewöhnlichen, 
nicht linearen Differentialgleichung, die keinerlei grundsätzliche Schwierigkeiten bietet, wenn 
auch zwischen Kopfwelle und Kegelmantel in Schallnähe Übergänge von Überschall- zu Unter- 
schallströmungen auftreten. Wenn das Problem der Überschallanströmung von Kreiskegeln für 
beliebig dünne Körper und für Machzahlen der Anströmung M„ beliebig nahe an Eins 
(M„-— 1) demnach grundsätzlich gelöst erscheint, solange nur die Kopfwelle anliegt, so besteht 
dennoch ein Bedarf für die Lösungen bei schlanken Kegeln in Schallnähe. 


Von Klaus Oswatitsch und Lars Sjödin 


reichen die vorhandenen Tabellen und Diagramme in keiner Weise aus. 


In diesem Gebiet 
Aus diesen lassen sich 
keine allgemein gültigen Formeln und Ähnlichkeitsgesetze für Schallnähe ableiten und es fehlt 
der Anschluß an die üblichen Methoden für schlanke Körper in Schallnähe. 


_ ausgeführten, aber nur in Berichtform 


zu machen. 


Dabei unterschei 


y ist es dann ni 


Ausgehend von der Kreiskegelströmung ist es da 
valenzsatzes von Osw.atitsch Aussagen über beli 
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Zur Wechselwirkung zwischen Rumpf und Tragflügel Bel Br 
- Überschallgeschwindigkeit De, : > 


Von A. J. Stetter in München a, . 
Zur Bestimmung der Überschallströmung um eine Rumpf-Tragflügelkombination in lineari- Wz 7a © 


sierter Näherung hat man einerseits die Methoden der „slender body theory“ und andrerseits 
einige strenge Verfahren von Ferrari,Morikawa,Nielsen u.a. Erstere liefern sehr 
schlechte Ergebnisse, letztere sind für allgemeinere Fälle durchweg recht kompliziert. 

In Verallgemeinerung einer Überlegung dieser Autoren wurde im Referat gezeigt, wie man 
ohne große rechnerische Schwierigkeiten zu praktisch recht brauchbaren Ergebnissen kommt. 
Die zunächst vorliegende allgemeine Anfangs- und Randwertaufgabe für die Wellengleichung 
wird dadurch vereinfacht, daß man zwar die Einwirkung des Rumpfes auf den Flügel voll be- 
rücksichtigt, die Rückwirkung des Flügels auf den Rumpf jedoch weitgehend vernachlässigt. 
Die vereinfachte Aufgabe läßt sich mit Hilfe des Kalküls der Distributionstheorie allgemein lösen 
und liefert in einem Teil des Flügels die strengen Strömungswerte, im Restbereich eine obere 
Grenze für die wahren Werte (Trennungslinie zwischen den beiden Bereichen: die vom Flügel- 
ansatzpunkt ausgehende Mac h linie). 


Für den Fall eines aus einem unendlich-langen Zylinder bestehenden Rumpfes wurden die 
Ergebnisse der Berechnungen an zwei Beispielen gezeigt. Der von der Rumpfeinwirkung her- 
rührende Zusatzdruck verläuft von seinem Maximum am vorderen Flügelansatzpunkt längs der 
von diesem ausgehenden Mach linie in einem ‚Grat‘. Ferner zeigt sich, daß für eine grobe 
Abschätzung der Wechselwirkung bei verschiedenen Flügeln mit Überschallvorderkante in 
einem größeren M a c h zahlenbereich eine einzige Tabelle genügt. 


Lignes de courant limites a la paroi pour l’aile droite 
d’envergure infinie derapee 
de Adalbert Oudart ä Paris 


Il s’agit du comportement des lignes de courant parietales de couche limite aux abords de la 
generatrice de separation du bord d’attaque et aux abords de la generatrice de decollement pour 
l’aile droite d’envergure infinie derap£e, en fluide incompressible. 


Que l!’on considere le fluide parfait ou le fluide visqueux, il existe une surface I qui delimite 
en amont de l’aile les deux nappes fluides enveloppant l’aile. Sur cette surface , les particules 
fluides decrivent des lignes de courant qui se deduisent toutes de l’une d’elles par une translation 
parallele aux gen£ratrices du cylindre constituant la surface extefieure de l’aile. 


On demontre dans les deux cas que ces lignes de courant sont pratiquement tangentes ä la 
generatrice de separation et s’y dedoublent en fourche en conservant cette generatrice pour tan- 
gente. En fait, theoriquement cette generatrice de separation est en general une asymptote. On 
ri aussi que la fourche, est 2,18 fois plus ouverte en &coulement visqueux qu’en &coulement 
parfait. 


Mutatis mutandis, la ligne de courant limite est tangente ou asymptote ä la generatrice de 
decollement: ce point s’Eclaircit egalement pour le cas du decollement laminaire, ce qui n’a pas 
echappe &Searsetä Wild. 


9 Lo A D . D . 
L’experience, m&me en fluide compressible, confirme les comportements theoriques precises 
par la presente communication avec laquelle les figures donnees par Searset H ayes sont 


grosso-modo en accord. A noter que la figure obtenue par Wild!) exige une modification au 
bord d’attaque. 


CE. H.Schlichting, Grenzschichttheorie p 164. Bild 108. 


RENNEN 


a die ee re Einflusses‘ des Br dien 
a inter Zugrundelegung des Prandtlschen Schubspannungsansatzes wird 
gem rung des logarithmischen Geschwindigkeitsgesetzes ein Wandgesetz abgeleitet, da 2 
fluß in Wandnähe Rechnung trägt. Die empirischen Koeffizienten des Wandgesetzes Mer 
eiden empirischen Koeffizienten des logarithmischen Gesetzes und vom Druckgradienten 
ler Strömung unabhängig. Das Wandgesetz kann unter Umständen dazu dienen, die Wandschub- 
spannung aus der TR Geschwindigkeitsverteilung i in Wandnähe zu bestimmen?). | 


Berechnungsmethoden von Lavaldüsen auf Grund eines 
Formelsystems für die Strömung in der Düsenkehle 
Von W. Haack und J. Zierep in Berlin 


w, d seien die Polarkoordinaten des Geschwindigkeitsvektors einer Gasströmung durch eine 
Lavaldüse, deren Meridian f(x) ist. In der Umgebung der Düsenkehle, wo der Übergang von 
Unter- zu Überschallgeschwindigkeit erfolgt, werden w, ® durch Reihen nach Potenzen von z, r/f 


‚dargestellt. Die so gewonnenen Formeln dienen als Ausgang für das Charakteristikenverfahren 


im Überschallgebiet. Als Anwendung wird die Methode zur Berechnung einer rotationssymme- 
trischen Parallelstrahldüse und einer Kegeldüse mit idealer Quellströmung und schließlich eines 
ebenen Krümmers ausführlich erläutert; dabei kann stets erreicht werden, daß die Beschleunigung 


j steuig ist. Ausführliche Veröffentlichung in Z. f. Flugwiss. 1954, Heft 2. 


1) Publikation i im Ing.-Archiv, Bd. 22 (1954). 


Referate folgender Vorträge sind nieht eingegangen: 


W.Bader, Berlin: Die Düsenströmung in vereinfachter Betrachtung. 


U. Domm, Berlin: Die Stabilität der Kärmänschen Wirbelstraße unter Berück- 
sichtigung realer Geschwindigkeitsprofile der Einzelwirbel. 


BrEichelbrenner‘, Paris: Eine Verallgemeinerung der Theorie der schlanken 
Körper. 


O0. Emersleben, Berlin u. Greifswald: Über eine doppeltperiodische Parallel- 


strömung zäher Flüssigkeiten (die Arbeit wird in dieser Zeitschrift erscheinen.) 


E. Eujen, Braunschweig: Bemerkungen zur Durchflußmessung mit Drossel- 
meßgeräten bei kleinen Re-Zahlen. 


K. Mangler, Farnborough: Berechnung der Druckverteilung für einen langsam 
schwingenden Flügel in der Nähe der Schallgeschwindigkeit. 


M.Schäfer, Göttingen: Über die stetige Rückkehr von Überschallströmungen in 
den Unterschallbereich bei gemischten Strömungsfeldern (eine ausführliche Arbeit über 
dieses Thema wird im ‚Journal of Rational Mechanics and Analysis‘‘ erscheinen). 


W. Wuest, Göttingen: Asymptotische Absäuge-Grenz- un. an längs- 
angeströmten zylindrischen Körpern. 


| Ferner ist das Referat des Vortrages von 


H. Gebelein, Bamberg: Vielfachkorrelation und Hauptachsentransformation 


nicht eingegangen. 
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Über das Spektrum bei Eigenschwingungen ebener 
Laminarströmungen 
Von D. Grohne in Göttingen ° 


(Aus dem Max-Planck- Institut für Strömungsforschung) 


Die Stabilitätsuntersuchung einer Laminarströmung nach der Methode der kleinen Schwingungen führt = 
einer gewöhnlichen Differentialgleichung vierter Ordnung auf eine nicht selbstadjungierte u e. 
Von den Eigenschwingungen dieser Randwertaufgabe hat man bisher vorwiegend die neutralen und angefae ge 
betrachtet. In der vorliegenden Arbeit soll nun dasgesa mte Spektrum der Eigenwerte c sowie die Gestalt der 
Eigenschwingungen untersucht werden. ; 


If the stability of a laminar flow is investigated by the method of small oscillations, a non-self-adjoint boun- 
dary problem must be solfed following from am ordinary differential equation of the fourth order. Till now, from 
the natural oscillations of the problem, preferably the neutral and the incited were sonsidered. The paper deals 
with the whole spectrum of the caracteristic values c and the shape of the natural oscillations. 


L’examen de la stabilite d’un courant laminaire suivant la methodes des petites oscillations conduüt a l’aide 
d’une equation differentielle ordinaire de quatrieme ordre a un probleme non adjoint de valeur ‚marginale. Des 
oscillations propres de ce probleme de valeur marginale on a consider jusqu’a present par pröference de neutres 
et siimulees. Dans l’article en question le spectre entier desvaleurs propres c ainsi quela forme des oscilliations 
propres va Etre examinee. 


MconexoBaHune YCTOÜYHBOCTH JNAMHHAPHOTO TeyeHUA IIO METOAY MA,JIEIX KolreÖaHnufi IPH- 
BOAHT, IPH OÖBIEHOBEHHOM AbHhPepeHuHansHOM YPaBHeHHH YeTBePToro IIOPAMKA, K HECaMO- 
cONPAKeHHON RpaeBoä sanaye. Ms yncıa BCex COOCTBEHHBIX KOJeÖAHHH 9TOH kpaeBoä sayaun 
AO CHX IICp IIPeuMyIecTB2HHO PACCMATPHBANHCh HEÄTPANBHBIe HU BOsÖy:KNeHHBIe. Mer HacTos- 
mei Pa6OTEI COCTOUT B UCCJIENOBAHHH BCETO CHEKTPA COÖCTBeHHEIX 3HAYEHHÄ c, a TAKıke 
hoPpMEL COÖCTBEHHBIX KONeaHHH. 


1. Einführung und Problemstellung 
Bei der Untersuchung einer ebenen Laminarströmung auf Stabilität gegenüber kleinen Stö- 
rungen beschreibt man eine stromabwärts (in Richtung der x-Achse) wandernde Störung der 
Stromfunktion durch den ‚Partialwellenansatz‘‘ 


u=o lg) N IT u 


und erhält alsdann für die Verteilung der Störamplitude @ (y) senkrecht zur Hauptströmung die 
sog. Stabilitätsdifferentialgleichung 4. Ordnung 


u 1 " 
W-).W- RU Lug +aig) ..... (2, 


worin U(y) das Geschwindigkeitsprofil der in z-Richtung fließenden Laminarströmung bezeichnet. 
Hinzukommen die hier auf den Parallelkanal bezüglichen Randbedingungen 


(1) 20, EDER 


welche das Haften auch der gestörten Strömung an den begrenzenden Wänden zum Ausdruck 
bringen. In diesen Gleichungen sind die Geschwindigkeiten U und c auf eine Bezugsgeschwindig- 
keit U, bezogen, ferner die Längen x, y und 1/« auf die halbe Kanalbreite 5, und schließlich die 
Zeit t auf die Zeiteinheit b/U,. DieReynolds sche Zahl R ist definiert durch - 
perl 
v 

Die aus Differentialgleichung und Randbedingungen bestehende Randwertaufgabe bestimmt 
zu jedem Parameterpaara und Rein Spektrum unendlich vieler Eigenwerte c,. Die zugehörigen 
Störungen (1.1) verlaufen gedämpft, wenn Im (c,) < 0, dagegen angefacht, wenn Im (c,) > 0 
ist; « wird als positiv reell angenommen. Eine Grundströmung heißt für einen R-Wert stabil, 
wenn das gesamte Eigenwertspektrum c, für alle möglichen Werte von« nur gedämpfte Stö- 
rungen enthält. Durch diese Begriffsbildung wird der Bereich der Reynoldsschen Zahl R in 
einen Stabilitätsbereich 0 < R < R* und in einen Instabilitätsbereich R> R* aufgeteilt, die 
beide durch die Stabilitätsgrenze R* voneinander getrennt werden. 

Da in der bisherigen Literatur fast ausschließlich neutrale, allenfalls angefachte Schwin- 
gungen untersucht worden sind, soll in der vorliegenden Arbeit auf Anregung von Herrn Prof. 
Dr. W.Tollmiendasgesamte Spektrum der Eigenwerte c, in seiner Abhängigkeit von «& 
und R untersucht werden, wobei zur Vereinfachung der Nachdruck auf die Abhängigkeit von «R 
gelegt sein möge. Eine geschlossene Lösung dieser Aufgabe ist in den beiden folgenden Sonderfäl- 
len möglich: 1. in dem mit U = 0 gleichwertigen Falle U = const. Es handelt sich um die bereits 
vonLordRayleigh behandelten „Schwingungen auf einer ruhenden Flüssigkeit‘. Die Lösung 


“are 


er 

ß ‘ L 

a 

ia a Ze 


Ge g | n v (y) anten de en „in den len 
analytisch-approximativ bestimmt werden: 


Tr im Grenzfall «R— 0 bei beliebiger Ordnung. n der Bigenwerte Om 
_ II. im Grenzfal n— o bei festgehaltenem «R. 
m. im Grenzfall «R— oo bei beschränkter Ordnung n. 


ae stetige Überführung der drei Fälle ineinander bei gleichbleibender Indizierung istin 
den obengenannten Sonderfällen U = 0 und U = y möglich. Die Indizierung der Eigenwerte u, 
kannin EN Fällen I und II nach wachsender Dämpfung ee d. h. nach der Vorschritt nee 


ge RE ER 


c Diese Vorschrift ist jedoch nicht immer auf den Fall DEE anwendbar, wenn "die Indizierung bei 
stetiger Änderung von a und R gleichbleiben soll. 


Die in (1.2), (1.3) formulierte Randwertaufgabe ist im allgemeinen nicht selbstadjungiert, 

- womit der Rückgriff auf die bekannten Aussagen und Abschätzungen der Sturm-Liou- 
villeschen Theorie entfällt. Die Eigenfunktionen bilden im allgemeinen kein Orthogonal- 
system. Sie bilden jedoch, wie 0. Haupt [3] gezeigt hat, unter gewissen Voraussetzungen ein 
Funktionensystem, welches bezüglich jeder den Randbedingungen (1.3) genügenden und viermal 
stetig differenzierbaren Funktion were ist, und in ein Orthogonalsystem umgewandelt wer- 

den kann. 


2. Die Grenzfälle « R— 0 bei beliebiger Ordnung rn der Eigenwerte c, sowie n— oo 
j bei beschränktem «R 
Wie bereits Lord Rayleigh [8] fand, läßt sich das gesamte System der Eigenfunktionen 
und Eigenwerte im Falle der Grundströmung U = 0, d.h. bei einem ruhenden Medium geschlos- 


sen angeben. Da diese Eigenwerte auch bei allgemeineren Grundströmungen zu approximativen 2 4 

Darstellungen geeignet sind, seien sie hier kurz hergeleitet. Im Falle der Grundströmung U = 0 Bi 

vereinfacht sich die Stabilitätsdifferentialgleichung (1.2) auf i A 
8W —- 2028’ +nt® +ioakR-0(8’—0)=0........ (2.1), * 


wobei die Eigenwerte zur Unterscheidung von den Eigenwerten c der allgemeinen Stabilitätsdiffe- _ 
rentialgleichung mit C bezeichnet seien. Die Gleichung wird gelöst durch jede der Funktionen f Be. 


E Dily) = Sofoy—Cofaygen Pu(y) = Sin oy— Sinay Zn ® u} E: 
wenn . 
2 r Pe KB a a ru a (2.3) 


gesetzt wird. Der Anteil ® (+ 1) = O der Randbedingungen wird identisch erfüllt. Die restlichen 
Randbedingungen ®’(+ 1) = 0 führen auf die diesbezüglichen Zweige der Eigenwertgleichung: 
oITgda=maTgw im Falle I, (24). 

a&otga =wlotgw im Falle II 


Die Gleichungen (2.4) haben zu positivem & keine Wurzeln » außerhalb der imaginären Achse der 
komplexen w-Ebene. Mit & ist auch— w ein Eigenwert zur gleichen Eigenfunktion. Es genügt 
deswegen, nur die positiv-imaginären Eigenwerte zu betrachten. Bezeichnet man gerade Eigen- 
funktionen mit geradem Index und ungerade Eigenfunktionen mit ungeradem Index, so mögen 
die GlIn. (2.4) die Lösungen 

Oy @y @g '- imFalleI und w, @, @, '-- im Falle II 


besitzen. Es läßt sich stets erreichen, daß die Be eine monotone Folge 
0< <= =. =<er 


bilden. Die zugehörigen Eigenwerte C,, lauten nach ES 


a 


a > nn in mern ALTE 
On t«R A 
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sie sind also im Sinne wachsender Dämpfung angeordnet. Der n-te Eigenwert läßt sich nach oben 
und unten abschätzen durch 


a2 +m(% +3) <iaR-<ar+m (a+ 1) ne (2.6). 


Aus der Darstellung (2.5) geht hervor, daß die Eigenwerte C sowohl für « R—- 0 als auch fürn — oo 
sehr groß werden. Dasselbe Verhalten findet nun auch bei allgemeineren Geschwindigkeitspro- 
filen U(y) statt, denn die Hauptanteile in der Stabilitätsdifferentialgleichung (1.2) werden dann 
durch die Gl. (2.1) dargestellt. Wir wollen diesen Gedankengang jetzt dadurch präzisieren, daß 
wir die Differenz c — © der Eigenwerte c bei allgemeineren Profilen gegenüber dem Eigenwert C 
der Gl. (2.1) einer genaueren Abschätzung unterwerfen. Unter Einführung der Differentialope- 
ratoren 


L[pl = U — a2 9)—U’p, Mlpl=p"— ap, Nlpl=p®—202p"+a!p. . (27) 
läßt sich die Stabilitätsdifferentialgleichung (1.2) in der Form schreiben 


©: Mg] = Lfpl— op: Nlgl. 


Entsprechend lautet 2.1) 0-M[d]= _n -N[Ö] 


Unter Ausnutzung der Selbstadjungiertheit der Operatoren M und N erhält man aus diesen 
beiden Gleichungen die Beziehung 


+1 +1 
(ce — YfpM [Bd] dy = /®Llpl 0 (2.9). 
Wird die für p noch freie Normierung durch die Vorschrift 
+1 4 
Sale -dy = /[PMIe] Ar, aule See a EE (2.10) 
festgelegt, so folgt nach Einführung der Hilfsgrößen 
+ +1 
_J BL19]y JOLlp— 2] dy 
= 7 0, g= = 7 ——.......(@1) 
JPMI®l dy [PM[9] dy 
E= -ı 
aus (2.9) die Darstellung 
a En a N Er 1 ul. (2.12). 


Hierin können CO und Q als vermöge der in (2.2) dargestellten Funktionen ® bekannt angesehen 


werden. Die Eigenwerte 0, wurden bereits in (2.6) eingegrenzt. Für Q erhält man direkt die Ab- 
schätzung 


ri 
0-3: |v-a+old), | ee re er (2.13). 


Für g erhält man in Verbindung mit einer simultanen Abschätzung der Funktion (o — ©)” 
— a? (p — ®) die Abschätzung 


Setzt man beides in (2.12) ein, so erhält man unter Beachtung von (2.5), (2.6) die beiden Teilaus- 
sagen 


Vz 


(2.152), 


(2.15b). 


Die letztere Abschätzung besagt, daß die Eigenwerte c, der Stabilitätsdifferentialgleichung bei 
genügend hoher Ordnung n schließlich gegen die (im Realteil um die mittlere Geschwindigkeit der 
Grundströmung vermehrten) Eigenwerte C,, der „Strömung Null‘ [vgl. (2.1)] streben. Eine 


a ee it 3 


ih daß bereits F. N oether 1) auf S. 239, Formel (28) eine ‚asymptotische Dar 


le — ©,| L A: (a RU, 


_ wobei ne ein Näherungseigenwert ist, der (in unserer Bezeichnungsweise) durch die Glerchung 


+ 
NT 


Im [VT=a): PUR 


Bene ist und etwa unserm in Gl. (2.1) eingeführten Näherungseigenwert C,, entspricht. In 
der vorstehenden Abschätzung vonMorawetz dürfen aber weder« Rnochn beliebig groß wer- 
den; denn im ersten Falle würden die Eigenwerte in die ausgeschlossene Nachbarschaft von 
c=w(y,) (w = U; y, bezeichnet die Wand) rücken, und im anderen Falle würde die Abschätzung 
inhaltslos, da das Verhalten der Größe A für unbegrenzt wachsendes » nicht mitgeteilt wird. 


3. Geradlinige Couetteströmung. Der Grenzfall « R> oo bei endlicher Ordnung n 

Im Sonderfall der ruhenden Grundströmung U = 0 wird das Verhalten der Eigenwerte c, 
für unbegrenzt wachsendes « R durch die Formel (2.5) beschrieben, in welcher die Größen w, 
nicht mehr von x R abhängen. Abweichend von diesem Gesetz ergibt sich bei allgemeineren Ge- 
schwindigkeitsprofilen ein Verhalten wie 


—U(-)$; für & ko | EEE (3.1), 


5 


worin die komplexwertigen Größen 7’, nicht mehr von & R abhängen. Werden diese Eigenwerte 
c„ durch stetigen Übergang von « und a R an die Eigenwerte (2.15a) angeschlossen, so geht bereits 
in Sonderfällen wie der geradlinigen Couetteströmung das Ordnungsprinzip (1.4) nach wachsen- 
der Dämpfung verloren. Will man also eine gleichbleibende Indizierung der Eigenwerte c, bei 
stetiger Überführung der Grenzfälle «R—0 und «R— o ineinander erzielen, so muß man 
diesen Vorgang auch tatsächlich durchführen, was eine geschlossene oder nur soweit approximative 
Lösung des Eigenwertproblems voraussetzt, als die einzelnen Eigenwerte noch voneinander unter- 
scheidbar bleiben. Eine Lösung in diesem Sinne ist uns nur im Sonderfall der geradlinigen 
Couetteströmung gelungen. Sie möge deswegen den Gegenstand des folgenden Abschnittes 
bilden. 

Nach Einsetzen des Geschwindigkeitsprofils U = y der geradlinigen Couetteströmungin 
die Stabilitätsdifferentialgleichung (1.2) ist diese mittels der Substitution 


a a ER (3.2) 


auf die Besselsche Differentialgleichung in der Nebenform 
y’—feaRy—o)+ly=0....... N 


reduzierbar. Um über die Randbedingungen zu den Eigenwerten zu gelangen,. werde (3.2) in der 


Form 


y r 
an au 2 nee a 
il 


umgekehrt. Die Randbedingungen  (— 1) =p’ (— 1) = 0 sind dann identisch erfüllt; die rest- 
lichen Randbedingungen @ (+ 1) =’ (+ 1) = O verlangen das Bestehen der beiden Gleichungen 


+1 Fi 
IE (y): Eng =  [wweotayay =. ER. 
—öt —ı1 


24* 


‚was anders definierte Eigenwerte bei unbegrenzt wachsender Ordnung angegeben hat, P 
| gs nur unter Andeutung eines Beweisganges, gs beträchtliche Schwierigkeiten yollzu 


- Erwähnt sei weiter eine von C. S.Morawetz [6], S. 580 angegebene Abschätzung für die 7 
Bigenwerte c ! 
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Die Differentialgleichung (3.3) ist mittels der Substitution 


y—y=en mit e=lah? ud wech: ee (3.6) 
in die Differentialgleichung er 
ar Er 347 
überführbar. 
Seien 
y(n) und. a Gm) ae RER 3), 
zwei geeignete Fundamentallösungen dieser Gleichung, und bezeichne 7-1, 7+1 die vermöge (3.6) 
den Wänden y = — 1 und y = + 1 zugeordneten Werte 
1 v0? -ı( 39 
=), HZ c fe ); 
so folgt aus (3.5) die Eigenwertgleichung 
71 n+1 N+1 Si "+1 
| in oe 
| Yı) m | Ynlm)Cof aen dn — |y6n) nn dy OS #endy=0 (3.10). 
n—1 Hl DT 71 


Zur weiteren Behandlung dieser Gleichung empfiehlt sich die Einführung einer Folge von 
Funktionen A,(n) durch das Laplace-Integral 


1 nztHi — 
x A,„(n) ==] e 2 E zn—1 , dz RE (3-11), 
o U 
in welchem der Integrationsweg X in derin Bild 1 gezeichneten Weise vom 
Unendlichen ins Unendliche läuft. Die Funktionen A,„(n) genügen der 
Differentialformel 
Bar AH...) 


und der Rekursionsformel 
vi Ans + N Antı +n: 4, =0 ER EI ER ROLE: 8:13); 


mittels derer die Gesamtheit der Funktionen A,„(n) und ihrer Integrale und Ableitungen aus den 
drei Grund funktionen 


An)» A,(n); AD 
rekursiv aufgebaut werden kann. 


Die Bedeutung dieser Funktionen A,(n) für das Stabilitätsproblem liegt darin, daß die 


beiden in der Eigenwertgleichung (3.10) benötigten Partikularlösungen der Differentialgleichung 
(3.7) darstellbar sind in der Form 


2ni 2ri 
vo) = AM), Yu) =e? A >) a HS 


Das Erfülltsein der Differentialgleichung (3.7) folgt aus den Formeln (3.12) und (3.13). Die 
lineare Unabhängigkeit der beiden Funktionen folgt aus dem Nichtverschwinden der ja kon- 
stanten Wronskischen Determinante an der Stellen =0. 

Die Grundfunktionen A,(n)und A,(n) = Aı(n) sind (in etwas anderer Bezeichnungsweise) 
vonH.H.Aiken [1] für ein quadratisches Punktgitter mit der Maschenweite 0,1 innerhalb 
des Kreises |n| <6 der komplexen n-Ebene numerisch vertafelt worden. Die Grundfunktion 
AH) = / A,(n) dn kann daraus durch eine numerische Integration beschafft werden. 


Außerhalb dieser Tafel kann das Verhalten der Funktionen A,(n) aus der asymptotischen 
Reihendarstellung 


a 
A,)zEn € Ge EEE REN GCHNN, 


v—=0 


erschlossen werden, welche für m — © im Winkelraum — = oeaınz u mit be- 


> Eig ellung Nullstellen 7y der Funktion A,(n) benötig 
tische Berechnung dieser Nullstellen für |y|> 1 ist mittels (3.16) nicht ur 
' möglich, da die Nullstellen aus dem Gültigkeitsbereich dieser Darstellung herausrücke 
Man umgeht diese Schwierigkeit durch Anwendung der zweiten der aus dem Integral (3.11 
 hältlichen Relationen | Par hl IRRE) 


“ 


RE | | | (, BREISRE er en. 
Braten Au\ne s) ze. ei .e ) RT EN BER (En PR 
FR gr angt | (* = konjugiert komplexe Werte) Se 


e und 
Er _- RE A se  ;ten ER ee 
B* a A,ı)+te ° le ae 3 er ;) Pl MRS (3.18), Zi kn. 
: worin die Polynome P,(n) vom Grade —n derselben Rekursionsformel (3.13) Be R 
= | nn RP = 0 I ee 3.19) A 


mit den Anfangselementen 
Be PENRANER 
genügen. I 
* Kombination der Formeln (3.16) und (3.18) liefert alsdann die für 2] > oo im Winkelraum 
 jargz| <r — ö gültige asymptotische Darstellung : 


‚5 cos E 282 1 2 
PR ( ei z1— = 
Yn: zes 


Hieraus folgt für die Nullstellen N „ von A, (nv) = 0, die nach (3.17) paarweise symmetrisch zur 


E Graden arg = 5/6 liegen, die asymptotische Darstellung 

Er | 77 = (Bmn)B.ert 0 9e® N 1 I ee 3.21). 

- Das niederste Nullstellenpaar wurde nach der Tafel von Aiken [1] zum Werte 

r Js + 0,2708] | r 
nn = 4,257 e Ba Reale) 2 


berechnet. 
Zur Weiterbehandlung der Eigenwertgleichung (3.10) empfiehlt es sich, die Funktionen 
Sinsenund@ojaeninihre Taylorreihen zu entwickeln und alsdann die Summation mit 
der Integration zu vertauschen, was durch einen Satz von Bromwich (vgl. [2], S. 398) ge- 
rechtfertigt wird. Die erhaltenen Reihen konvergieren nach Sätzen der Laplace- Transfor- 

mation für jeden Wert von ae. Werden diese Reihen nach den ersten Gliedern abgebrochen, - 
mit Restgliedern versehen und in die Eigenwertgleichung (3.10) eingesetzt, so vereinfacht sich 

diese deswegen unter Beachtung von (3.17) auf | 

An) — 40-1) » A) Al mi) 

4) — 4: 9-ı) » 3) Ani) 

Was geschieht jetzt, wenn « R über alle Grenzen wächst, d. h. wenn e— O strebt? Wegen 

der aus (3.9) folgenden Beziehung 7, — n-ı = 2/e muß mindestens eine der beiden Größen n,, 7-ı 

für e— 0 auf einer Parallelen zur reellen Achse ins Unendliche streben. Es genügt, dies bezüglich 

7, anzunehmen, denn im anderen Falle würde alles spiegelbildlich zur imaginären Achse der 

n-Ebene verlaufen (wie ja grundsätzlich bei der Couetteströmung mit C auch — C* ein Eigen- 

wert ist). Unter Beachtung des asymptotischen Verhaltens (3.16) vereinfacht sich dann die 

Eigenwertgleichung (3.24) auf 


=0(8 2) für |mjsconst (3.24). 


‚4 (N-1) € . 2 
i een iin. Ber. 2. . (3.25). 
4.n-1) 2 e) 
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An dieser Formel ist erkennbar, daß n-ı für &— 0 gegen die in (3.21) abgeschätzten Null- 
stellen 7y der Funktion A, (n) streben muß. Man erhält so für die Eigenwerte c unter Beachtung 


von (3.9) die asymptotische Darstellung 


a a ; RE 
mitnyaus Aulnn) = 0 


Hiermit haben wir die angekündigte Eigenwertformel (3.1) am Sonderfall der Couette- 
strömung nachgewiesen. 

Um den Verlauf der Eigenwerte c über den ganzen Bereich 0 <aR <w zu verfolgen, 
muß man auf die Eigenwertgleichung (3.10) oder ihre Näherungsform (3.24) zurückgreifen, wobei 

die Funktionen A,(n), Aı(n) As(n) als bekannt 
MT Iyorauszusetzen sind«. Wir Haben iermsch ale 
12 niedersten Eigenwerte c in ihrer Abhängig- 
keit von aR für einen festen Werta = 1 be- 
rechnet und in Bild 2 dargestellt. Die Verän- 
.derlichkeit der Eigenwertkurven mit «& ist nur 
gering und wird z.B. für« R> & mit. klein 
von der Ordnung 0(a? e2). 

Das Eigenwertspektrum der geradlinigen 
Couetteströmung wurde bereits von L.Hopf 
[5] diskutiert. Hopf ersetzte etwa auf der 
Stufe unserer Eigenwertgleichung (3.10) die 
von ihm durch Hankelfunktionen der Ord- 
nung 1/3 dargestellten Lösungen %,, %7ı durch 
die ersten Glieder ihrer asymptotischen Reihen 
(3.16), wodurch sich die Eigenwertgleichung 
auf eine algebraische Gleichung aus Hilfs- 
argumenten sowie Kreis- und Hyperbelfunk- 
tionen vereinfachte. Da Hopf jedoch bei 
den asymptotischen Darstellungen der Han- 

Bud2. Geradlinlas Con ettenkrömung. Die rwält niedsfsten ' relfunktionen. gewisse, Versehensurker un 
sind, bedürfen seine Ergebnisse teilweiser Kor- 
rekturen. Fallen diese Änderungen bei kleinen Werten von « R kaum ins Gewicht, so ändern 
sich z.B. die Werte von ny in einer (3.26) entsprechenden Formel schon erheblich. Auch der 
topologische Zusammenhang der Eigenwertkurven c=c(a,aR) sieht uns anders aus als bei 
Hopf. Erhalten bleibt jedoch das qualitative Bild der Eigenfunktiönen, die daraus gezogenen 
De len Folgerungen und das Hauptergebnis, daß nämlich alle Schwingungen gedämpft 
verlaufen. 


4. Der Grenzfall a R—> & bei endlieher Ordnung r für symmetrische Grundströmungen 
Bei einer Grundströmung mit symmetrischem-Geschwindigkeitsprofil 


UW)=-Uemir er er (4.1) 
besitzt die Stabilitätsdifferentialgleichung (1.2) stets ein Fundamentalsystem von vier Lösungen 


91 9 Pz 94 so daß r 2 
®(Y)»  B4(y) gerade Funktionen von y | 


und %(Y); ®(y) ungerade Funktionen von y| ee (42) 


sind. Soll nun eine Linearkombination dieser Lösungen den Randbedingungen (1.3) in der Reihen- 

folge g— 1) =0,9(—1) =0,9(+1) =0,9(+1)=0 genüge leisten, so muß die folgende 
- unter Beachtung der Symmetrien (4.2) vereinfachte Determinante verschwinden: 

el Er Et 

Da Pı (—1) 9; (1) P, (1) 9, =D) 

AAN HN) Ba) A 

AA) Bar BA) — 5) 

x; diese Determinante als Produkt der beiden zweireihigen Determinanten geschrieben werden 

ann: E 


=), 


D 19, (=) D OR (—1) ; 


iR BCM] 
91), (I) 


P, ei) ’ ®, (—1) Be 


4 


TEEN OP TERIUHER 


+ Ableitungen nach y gegen eine Grenzfunktion 


= streben, führt zu der sogen. „reibungslosen Differentialgleichung‘“ 


(ER ED TE BR ne 3 


?, 50 bestimmen, daß sie geeignete 

lich sind. Als dem lem angepaßt erweisen sich die vo: 
ptotische Integration der Stabilitätsdifferentialgleichung‘“) beschrie 
en, deren asymptotische Darstellungen mit Restgliedabschätzungen versehen 


m d © > 


po 


3. Um den Anschluß dieser Fundamentallösungen an die unsrigen herzustellen, ist es une- 
lich, zuvor auf den Begriff der „reibungslosen Approximation‘‘ einzugehen. Die Frage nach 


ösungen der vollständigen Stabilitätsdifferentialgleichung (1.2), die füra R> © mitsamt ihren 


der solche Grenzfunktionen notwendigerweise genügen müssen. Will man die Lösungen dieser 
reibungslosen Differentialgleichung zur Approximation der Lösungen der vollständigen Differen- 
tialgleichung für & R— oo verwenden, so darf man den Geltungsbereich der Grenzwertaussagen 
(4.4) in der komplexen y-Ebene nicht außer Acht lassen. Nach W. Wasow [13] gilt darüber 
der folgende Satz: | ET | 
„Von den vier Fundamentallösungen (4.2) sind je eine gerade und eine ungerade so be- 
stimmbar, daß mit zwei geeigneten reibungslosen Lösungen 7,(y) und 7,(%) die Approximationen 


& 1 
AN) euy)H0 = » %ı(y) = ungerade Fkt. von y, (4.6) 


= 1 - E 
9%(y) =x(y) +0 nn » %(y) = gerade Fkt. a Y 


in jedem festen Innern eines Doppelbereiches (I -+ II) oder (IT + III) oder (III u D zu Recht 
bestehen und jeweils im komplementären dritten Bereich III bzw. I bzw. II ungültig werden. 
Das Gleiche gilt für die Ableitungen nach y.‘‘ (Vgl. Bild 3.) 


Die jeweilige Grenze zwischen den Bereichen I, II und III 


genügt der Gleichung 


Re I% J Yu=eay = Ba 


wenn %; die durch Noll 2 
nee, AT) 


definierte ‚‚kritische Stelle‘ bezeichnet. Näheres zu den Ge- pila 3. Die Gebiete I, IL und IILin der 
bieten I, IIu. IIIbei Wasow [13]. komplexen y-Ebene 


Die reibungslose Differentialgleichung (4.5) besitzt an der kritischen Stelle U(yx) =c 
eine singuläre Stelle der Bestimmtheit. Zwei Fundamentallösungen verhalten sich dort wie 


re D% 
4ly) = (y— u) Bıly u), a d=RWY— Yu) + v Bı(y — Y) : (y — Yı) In (y—yr), (4.9), 
wenn ®, und ®, Potenzreihen mit dem Anfang 


BO=142 4,400. Pe) =1+0(2) TEN BER) 


bezeichnen (vgl. W. Tollmien [12], S.35). Der gemeinsame Konvergenzradius dieser 
Potenzreihen wird beschränkt entweder durch den Konvergenzradius einer entsprechenden Reihe 
für U—c oder durch die nächstbenachbarte Nullstelle von U — c als einer singulären Stelle der 


Differentialgleichung. 


ay.ln - a 
Li} Kar 
» Wi 


| „im oa R)=x(). EURE (4.4) 


% u, 
3 A 
iv 


Pr 
A 


b Ber 
N 
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Für den weiteren Verlauf empfiehlt sich die Einführung einer Folge von Funktionen B,„ (n) 
durch das mit (3.11) vergleichbare Laplace-Integral 


a 2 
B,0) = ne Eee ee (4.10). 


Darin bezeichnet © = 0,5772... die Eulersche Konstante. Der Integrationsweg ® verläuft 

in der in Bild 4 angedeuteten Weise in der komplexen längs (0 , — i 00) aufgeschnittenen 2-Ebene 

vom Unendlichen ins Unendliche. =D 
Die Funktionen B,(n) genügen der Differentialformel 


dB 
u Bel) a San (4.11) 
N 
a Euranse und der Rekursionsformel 
$% i- Batzs+n' Bor tn Bn=P, - N . (4.12), 
| in welcher P,(n) die in (3.19) erklärten Polynome sind. 
Bild 4. Integrationsweg ® Mittels dieser beiden Formeln kann die Gesamtheit der Funktionen 


in der komplexen lä : : ! : 
(0° fe) nilfgeschnittenen B,(n) und ihrer. Ableitungen und Integrale rekursiv aus den drei Grund- 


ann funktionen 
EROBERN De, en ee ee . . (4.13) 
aufgebaut werden. Die Grundfunktionen B, und, sind mittels der Darstellung 
A), A n®), 
B,o) =2ri|, ; nr 4.13a 
dr a, ram) a 


(*—=konjugiert-komplexer Wert) auf die Funktionen A, zurückführbar (vgl. W.Tollmien[10], 
S.27). 
Die Bedeutung der Funktionen B,(n) für unser Stabilitätseigenwertproblem liegt darin, daß 
die Funktion B,(n) zufolge (4.11) und (4.12) der Differentialgleichung 
i Sn ren 1 
di N a Sg A er 


genügt, die unter der Bezeichnung ‚Differentialgleichung für die Reibungskorrektur‘‘ von 
W.Tollmien [10], [12] als wesentlieher Bestandteil in die asymptotische Integration der Sta- 
bilitätsdifferentialgleichung eingeführt worden ist. 

Nach diesen Vorbereitungen wenden wir uns den vier von Tollmien im Hinblick auf 
asymptotische Entwickelbarkeit konstruierten Fundamentallösungen 97, 97 Pr Yıy der voll- 
ständigen Stabilitätsdifferentialgleichung zu. Diese vier Lösungen lassen sich nach W. Toll- 
mien[12] S. 77 unter Verwendung der Substitution 


yzy=en mit e=(aRU,)® N und y aus Uly)=c zeezatıe 


so bestimmen, daß sie sowohl in einem festen der (aus Gründen der analytischen Fortsetzung) 
komplexen »-Ebene als auch in jedem fest en Innern des Bereiches II der reibungslosen Approxi- 
mation [vgl. (4.6a)] die folgenden asymptotischen Darstellungen besitzen: 


EA EI en 

Ur \ 

YuYy) =Rlen)+ U; Bılen) e[Bıam)+n-+nine] + 0(eIne); 

verschärfend gilt in jedem festen Innern von II | .  . (4.17). 
Pı(ly) = X2(y) +0) 


Des weiteren gilt 


Pr (y) = A-ı (m) +0), "undzwarin n]= constr. 2.2 2er 


as, gilt nach W. Wasow [13] quotientasymptotisch in jedem festen Innern von II 
vgl. (4. 


y 
5/4 PEREN  . 
Yırı(Y), Pıv (y) 7 const - en - exp Nor. | (7=e lc): 
%k 


Entsprechende Formeln gelten für die Ableitungen. 


| 
Da il a a a 


Ay au Ss den reibungslosen Lösungen %, und 72 vgl. @. A Des weiteren eilt 2 


en a 2 n ; 4; | = 2 are Eee: 20, 
Br ro Sala (Opel SE Er 


er #s, u) = ke + oles Se el: 1 A-ı(n) a. 00 a nee: 

_ Entsprechende Formeln gelten für die Ableitungen. 
Schreiben wir jetzt die beiden in (4.3) erhaltenen Eigenwertgleichungen als Produkt ne 

Faktoren, z.B. Dur 


(abe N) Br ; 
0=D=[9,—1 ee 1 Be 15 Bee. 
= beliebige Konstante), so liefern die Nullstellen 'der beiden ersten Faktoren keinen Beitrag. Be 


zum Eigenwertgebilde, da sie durch entsprechende Pole des dritten Faktors kompensiert werden, es 
sei denn, daß gleichzeitig die Ableitung #3(— 1) verschwände. Es genügt deswegen, nur die Null- 
stellen des 3. Faktors rare Nach Einsetzen der Approzimabionen (4.21) erhält man so 


Be a inet 142 + But) 
eine Pur — | ———— +0 (en) (n=n-) 
Yıı 14:7 Elndne+14 819) + B- Ve 


oder 


me +14 S1a+ Bon) 
0 es 2 4 ee Holen) | (29) 
— 1+eglnlne+1+S12+ Bm] 


mitn =n-ı =(—1-— yr)/e und 8,, 8, gemäß (4.21a) aus den reibungslosen Lösungen. Die aus 
einer nächsthöheren Näherung in (4.21b) Rap asende Funktion y(n) lautet 


oe 0 ne] | Be 
TH dn\ A-ı(n) 
und kann mittels der Formeln (3.12), (3.13) auf die drei vertafelten Grundfunktionen A,, A,, A, 
zurückgeführt werden. 


Wie verhalten sich nun die Eigenwerte c, wenn wir in der Eigenwertgleichung (4.22) R> 
d. h.e— 0 streben lassen ? Offenbar strebt 7_ı dann gegen die Nullstellen 7y der Funktion A,(n).' 


_ Durch Taylorentwicklung um diese eg folgt genauer 


FE URD) 
N— nt ZA > -elne+0(e). 


Die Eigenwerte c verhalten sich dann asymptotisch wie 


c— U(—1) = Urena—gi Ü un 2) aine +06) (4.23). 
mit n„ aus A,(nn) = 0 
In Ergänzung zu (3.21) seien hier einige Nullstellen y und Werte s A,/A, angegeben: 
PATER NN | v A,nnlAdı(nn) . ä s 
ee lee 
24952 09.3.037.|.1,902. 7. 2.. 0,851 ‚dem — 2 4.24) 
BR ar Ei; Anm) i pe 


Er 


Bet, 
- 


gestattet. Hierzu muß a 


1 


schaften „‚gerade‘‘ oder „ungerade‘‘ der Eigen! 

den Konstanten $, und 8, berücksichtigt werden. 

Auch zur numerischen Berechnung der Eigenw 

wir die Eigenwertgleichung in der Form (4.22) herangezogen ] 
Poiseuilleströmung und eine Strömung mit Wendepunktprc 


Bild 6. Wendepunktprofil 


v=(V2-1ı) +(2—V2)cosy3 n/a. 
Die vier niedersten Eigenwerte cin Abhängigkeit von R 


Bild 5. Ebene Poiseuilleströmung. Die vier niedersten heit 05 
=0,5. 


Eigenwerte cin Abhängigkeit von R bei& = 0,87 


vier niedersten Eigenwerte in Abhängigkeit von R für einen festen Wert von & in den Bildern 5 
und 6 dargestellt. Die numerische Rechnung selbst ist nach Zurückführung der in der Eigenwert- 
gleichung enthaltenen nichtalgebraischen Elemente auf die drei vertafelten Grundfunktionen 4,, 
A,, A, und auf die reibungslosen Lösungen eine Aufgabe der numerischen Praxis, auf deren Ein- 
zelheiten hier nicht eingegangen werden kann. Erwähnt sei nur die folgende Näherungsdarstellung 
der reibungslosen Konstanten $, 


U; (U 7 
! 0 

01 U — 0)? dy 

1 


{vgl. W.Tollmien [11], S. 100), welche mit Vorteil bei kleinen Werten von a anzuwenden ist. 

Die Indizierung der aus (4.22) gewonnenen Eigenwerte c im Sinne eines stetigen Anschlusses 
an den Grenzfalla R— OÖ bleibt hier ein offenes Problem. Auch im Geltungsbereich von (4.22) 
allein läßt sich eine durchgängige Indizierung nach der Vorschrift Im (cn+1ı) < Im (c,), d.h. 
nach wachsender Dämpfung, grundsätzlich nicht durchführen. Zwar lassen sich die Nullstellen in 
(4.23) nach wachsendem Imaginärteil ordnen, aber die Kurven Im (c,) können einander bei 
Anderung von e durchdringen. 


-8,=4:0-+0() für «<1 mit 4= en wln 


5. Die reibungslosen Eigenwerte innerhalb des Grenzfalles « R> «. 
Die Ermittlung der angefachten Eigenwerte 


Die Approximation (4.6) mittels der reibungslosen Lösungen sei entweder auf den Doppel- 
bereich I+ II (vgl. Bild 3) oder auf den Doppelbereich II+ III zugeschnitten, wodurch in den 
reibungslosen Lösungen jeweils das logarithmische Glied eindeutig festgelegt wird. Unter Ver- 
wendung der Approximationen (4.17) erhält man dann über die Eigenwertgleichungen (4.3) 
Eigenwerte c, die für R— oo gegen die sog. „reibungslosen Eigenwerte‘‘ c() (x) streben, welche 
durch die Randbedingung 4,(— 1) = 0 bzw. %,(— 1) = 0 der ungeraden bzw. geraden reibungs- 
losen Lösungen %,, %, definiert sind. Über diese durch den Geltungsbereich der Grenzwertaus- 
sagen (4.6) eingeschränkten reibungslosen Eigenwerte lassen sich nach W. Tollmien[11]z.T. 


oc a ü, 


Be Wendepunktprofi we ne er Ha) 2 “ Y. Die nur zu ae Eigen- Er 


funktion Enden reibungslosen Eigenwerte c sind in den Bildern 7 und 8 dargestellt. Der ” 5 


Bild 7. Ebene Poiseuilleströmung. Der f ; ' 
reibungsiose Eigenwert c zu gerader i Bu BU Wendepunktpzofil 
Eigenfunktion in Abhängigkeii von & u =(V2-1ı) +(2— Y2) cos y3 n/a. 


Der reibungslose Eigenwert c zu gerader Eigen- - 
- funktion in Abhängigkeit von & 


Existenzbereich dieser Eigenwerte ist jeweils durch ein IntervallO <« < const gegeben. Für den 


reibungslosen Eigenwert c hat Tollmien [11], S. 100 die er Näherungsformeln auf- 
gestellt: 


U (— U7 e 
(EA CENE — a) i seele) a mit Ulyr)=0 BE (5-1). 


r (U’—c,)2 dy 


Wir fragen nun nach dem Zusammenhang, in dem die reibungslosen Eigenwerte mit den bisher 


- diskutierten stehen. Die geschlossene Lösung bei der Couetteströmung vermag hierzu keine 


Antwort zu liefern, weil die betreffenden reibungslosen Eigenwerte dort garnicht existieren. Es 
ist aber möglich, die reibungslosen Eigenwerte in die Gl. (4.22) einzufügen, und sie so als Grenzfall 
innerhalb der Eigenwerte (4.23) zu verstehen. | 

Es werde also an der Eigenwertgleichung (4.22) der Grenzübergang x R> »,d.h.e—0 
bei festgehaltenem (— 1 — Yı) = en-ı vollzogen, wobei die Berechtigung dieses Vorganges aus 
der (4.22) vorangehenden Gleichung zu nehmen ist. Mittels der asymptotischen Formeln (3.16) 
folgt dann für« R> © 


BEL NT Ihe, 


Das sind aber genau die ersten Taylor glieder der reibungslosen Eigenwertgleichung y,— 1)=0 


bzw»x,(— 1) =, dieman gemäß der Bedeutung (4.21a) von S,, 8, bei Taylorentwicklungim 


Sinne der Reihen (4.9) erhalten würde. 

Auf Grund dieser Erkenntnis kann die Ermittlung der angefachten Eigenwerte (Im c > 0) 
vereinfacht werden. Da angefachte Eigenwerte nach den Ergebnissen des 2. Abschnittes nur 
innerhalb der ersten endlich vielen Eigenwerte und für genügend große Werte von R auftreten 
können, genügt es, die Gl. (4.22) auf angefachte Eigenwerte hin zu untersuchen. Und hierbei gilt 
dann die Alternative: Angefachte Eigenwerte sind entweder durch die reibungslose Randwertauf- 
gabe in Verbindung mit genügend großen Werten von R (approximativ) bestimmbar, oder sie 
liegen in einer Umgebung von c = U(— 1) und sind mittels einer der Gln. (4.22) oder (4.23) als 
zu endlichen Werten von n_ı gehörig bestimmbar. 

Wie hieraus für genügend große Werte von R folgt, aber in den Beispielen auch für die 
kleineren Werte von R bestätigt wurde, ist bei Wendepunktprofilen die stärkste Anfachung stets 
mit dem reibungslosen Eigenwert bzw. dessen Fortsetzung nach kleineren. R-Werten verbunden. 


L) 
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Hieraus folgt die wohlbekannte Tatsache, daß bei Wendepunktprofilen das Stabilitätsverhalten 
bereits aus der reibungslosen Differentialgleichung allein zu erschließen ist. Man vergleiche hierzu 
die Berechnung der reibungslosen Eigenwerte von G. Rosenbrook [9] für ein Wendpunkt- 
profil, welches dieser in einem divergenten Kanal gemessen hatte. 


6. Gestalt der Eigenfunktion. Die innere Reibungsschicht. 
Der zeitliche Verlauf einer Störung 
Um den zeitlichen Verlauf einer Störung zu beurteilen, wird man diese in Partialwellen der 
Art (1.1) zerlegen. Es ist dann erforderlich den Verlauf der Amplitude g(y) über die Kanalbreite 
hinweg zu kennen. Wir betrachten hier nur den Fall sehr großer Werte von«a R. 
Bei der Couetteströmung folgt über (3.4) in der Bezeichnungsweise (3.6), (3.9) und (3.15) 
füre<1,d.h.a R> 1 als Näherungsausdruck für die Eigenfunktion 
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Die Randbedingungen 9(+ 1) = 0 sind identisch erfüllt; die restlichen Randbedingungen sind 
mit der Eigenwertgleichung (3.24) identisch. Wie bereits aus der Differentialgleichung folgt, ist 
9*(— y) eine Eigenfunktion zum Eigenwert — c*. Wir haben hiernach für« =1 und R = 10° 
die Eigenfunktion zum Eigenwert c = —0,7 —i 0,3 berechnet und in Bild 9 dargestellt. 

In diesem Bilde fällt auf, daß sich die wesentlichen Änderungen der Eigenfunktion in einer 
Schicht —1 < y < y, abspielen, die etwa durch den Winkelraum arg <x/6 des starken An- 
wachsens von A_,(n) zu definieren wäre. In der Veränderlichen y lautet diese ‚innere Reibungs- 
schicht‘“ zufolge (3.6), (3.26) näherungsweise 


bey il tel Bens  Jalmn): ei me a 


Aus dieser Darstellung geht hervor, daß die Breite der Schicht mit wachsender Ordnung n der 
Eigenfunktionen zunimmt, wobei gleichzeitig die Dämpfungsgröße anwächst. Die Geschwindig- 
keit der zugeordneten Störwelle ist etwa gleich der 
Geschwindigkeit der Grundströmung in der Mitte 
der Schicht. Die Dicke der Schicht geht ferner mit e 
gegen Null. Die physikalische Interpretation dieses 
Sachverhaltes besagt nach Hopf [5], S.57, „daß 


Bild 9. Geradlinige Couetteströmung,. Eigenfunktion Bild 10. Ebene Poiseuilleströmung. Eigenfunktion „9 (%) 
o'(y)bix=1l1undR= 10° zum Eigenwert bei @=1 und R=7,7.10° zum Eigenwert 
ce = — 0,70 —i 0,30 c = 0,178 —i » 0,049 


jede beliebige Störung für große Werte von R so abklingt, daß schließlich nur von den Wänden, 
und zwar ohne gegenseitige Beeinflussung Störungen auszugehen scheinen, ein Verhalten, wel- 
ches an reibungslose Flüssigkeiten gemahnt‘‘. 

Bei allgemeineren Grundströmungen hat Tollmien[12]einen Näherungsausdruck für die 
Eigenfunktion aufgestellt, an dem sich im Dämpfungsfalle wieder eine ‚innere Reibungsschicht‘“ 
erkennen läßt, die durch den Winkelraumr/6 < argn = 5/6 der starken Wachstumsänderungen 
a bzw. A-ı(n) zu definieren wäre. In der Veränderlichen y lautet diese Schicht zufolge 

+Y3c _ —Y3e; 
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gehörigen reibungslosen Eigenwert ab.. ö 

‚In Gegenüberstellung der inneren Reine schnehe zur Grenzschicht kann man sagen, daß > 
_ die Grenzschicht denjenigen Strömungsbereich darstellt, in dem das Verhalten der laminaren 
Grundströmung maßgeblich durch die innere Reibung beeinflußt wird, während dieinnere 


eispiel der Pois euille Strömung zu 1 und R = = 2 7: 10 beree 


und i in Bild 10 dargestellt. Der Eigenwert ce = 


Reibungsschicht den Bereich angibt, in dem die Störung maßgeblich unter dem Einfluß der 
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BUCHBESPRECHUNGEN 


M.N. Balschin) Reibungsuntersuchun- 
gen an geschmierten Oberflächen 
bei tiefen Temperaturen (mit einem 
Vorwort von Prof.Dr.-Ing. Dr.-Ing. E.h. E.Heidebrock, 
Dresden). VIII + 64 S. m. 46 Abb. u. 3 Zahlentafeln. 
Berlin 1953. Akademie-Verlag. Preis geb. 7,50 DM. 

Das Verhalten von Schmiermitteln bei tiefen Tem- 
peraturen bestimmt nicht nur die Anlaßmomente 
beim Starten von Verbrennungsmotoren, sondern 
auch die Funktion von Meßinstrumenten im Flug- 
und Verkehrswesen in der Kälte. Planmäßige Unter- 
suchungen nach dieser Richtung lagen noch nicht vor; 
jetzt wird erstmalig von M. N. Butschin dieser Auf- 
gabe eine eingehende Experimentalar beit gewidmet. 

Bei einleitenden Versuchen an ungeschmierten 
Flächen zeigte sich eine Abhängigkeit der Reibung 
vom Vorhandensein adsorbierter Schichten. Bei ge- 
schmierten Flächen ergaben sich für jedes Schmier- 
mittel bestimmte kritische Temperaturen, unterhalb 
deren die Reibung sehr große Werte annahm, größer 
als bei den Versuchen ohne Schmiermittel. Diese 
kritische Temperatur gilt daher als untere Grenze für 
die Brauchbarkeit, sie steht nicht in engerer Verbin- 
dung mit dem Stockpunkt, den man bisher als Kenn- 
zeichen für die Verwendung bei tiefen Temperaturen 
ansah. 

Für die Versuche wurde eine Apparatur mit schiefer 
Ebene verwendet, also der Reibungswinkel gemessen, 
bei dem die Bewegung eintritt. Leider fehlen alle 
Angaben über die Abmessungen der Versuchsappa- 


ratur. Die Ergebnisse sind dennoch äußerst aufschlu ß- 
reich. Ferner wurden kleine Wälzlager untersucht, 
wie sie in Meßinstrumenten Verwendung finden. 
Zwischen den Gleitversuchen und denen mit Kugel- 
lagern ergeben sich in bezug auf die Wirkung der 
Schmiermittel enge Zusammenhänge, 

Die sorgfältig durchgeführten Versuche verdienen 
besondere Beachtung und geben hoffentlich Anreiz 
zur Erforschung vieler noch offener Fragen über das 
Reibungsverhalten von Maschinen und Apparaturen 
bei großer Kälte. 


Göttingen. G.Vogelpohl. 


Dr.-Ing. A. Nasvytis, Die Gesetzmäßig- 
keiten kombinatorischer Technik. 
(Wissenschaftliche Normung. Schriftreihe, heraus- 
gegeben mit dem Seminar für Technische Normung 
an der Technischen Hochschule Hannover. Heft 3) 
VIII +103S. m. 29 Abb. Berlin/Göttingen/Heidel- 
berg 1953. Springer-Verlag. Preis geh. 9,— DM. 

Der Verfasser hat es sich zum Ziel gesetzt, die 
mathematischen Grundsätze herauszuarbeiten, nach 
denen die verschiedenen Einzelfragen der Normung 
möglichst einheitlich zu behandeln sind. Dabei soll 
der Begriff ‚‚Normung‘‘im weitesten Sinne verstanden 
werden; er soll z. B. auch die Gesetzmäßigkeiten einer 
Sprache, den Aufbau einer Schrift und die Ordnung 
des Zahlensystemes mit umfassen. Es geht immer 
darum, aus gewissen Elementen — z. B. Grundmaßen, 


BEN NE RER 
El Getriebeteilen, 


und Ziffern 


mühen, die Zahl der zum Aufbau der gewünschten 


Gebilde benötigten Elemente nicht zu groß werden 


zu lassen. Im allgemeinen widersprechen diese beiden 
Forderungen einander, so daß man nach einer opti- 
malen Kompromißlösung zu suchen hat. 

Da hierbei in der Hauptsache kombinatorische 
Fragen auftreten, werden in Abschn. 21 zunächst die 


 Grundtatsachen der Kombinatorik zusammengestellt. 


Ausführlich werden dann in 22 bis 24 Kombina- 
tionen betrachtet, die in der zeitlichen Aufeinander- 
folge einzelner Elementarvorgänge bestehen — der 
Verfasser nennt sie ‚kinetische Kombinationen“. 
Zwar wird durch die Betrachtung solcher zeitlicher 
Folgen kein grundsätzlich neues Element in die 
Kombinatorik eingeführt, es treten aber einige bisher 
weniger beachtete kombinatorische Aufgaben mehr 
in den Vordergrund. 

Auf andersartige kombinatorische Fragen führt in 
25 die Untersuchung räumlicher Kombinationen. 


Die allgemeinen Grundsätze für die Behandlung von 


Normungsaufgaben werden nun in Abschn. 3 zuerst 
auf das Zahlensystem angewendet. Insbesondere 
werden Systeme mit verschiedenen Grundzahlen ver- 
glichen, und dabei zeigt sich, daß das Zehnersystem 
einen recht günstigen Kompromiß zwischen den beiden 
Forderungen liefert, einerseits die Ziffern-Anzahl und 
andererseits die. Stellenzahl klein zu halten. 

Hieran schließt in Abschn. 4 unmittelbar eine all- 
gemeine Untersuchung ‚‚additiver‘“ und ‚‚multipli- 
kativer‘“‘ Größensysteme an, die für die Schaffung 
genormter und möglichst vielseitig verwendbarer 
Elemente grundlegend sind. 

Die Abschnitte 5 bis 7 behandeln spezielle An- 
wendungsgebiete: Schriftsysteme, Baukastensysteme 
und Kombinationsmaschinen. 

Bemerkenswert ist der in dieser Schrift zutage 
tretende Umstand, daß technische Fragestellungen 
in Teilgebiete der Mathematik hineinführen können, 
die den Anwendungen zunächst scheinbar fernliegen. 


Hannover H. Epheser. 


Paul F. Byrd and Morris D. Friedman, Hand- 
bookofElliptie Integrals for Engi- 
neers and Physicists. (Die Grundlagen 
der mathematischen Wissenschaften in Einzeldar- 
stellungen. Bd. LXVII.) XIII + 355 S. mit 22 Abb. 
Berlin / Göttingen / Heidelberg 1954. Preis geh. 
36,— DM, geb. 39,60 DM. 

Allen theoretisch arbeitenden Ingenieuren, nicht 
nur denen, die wie die Verf. sich mit Aerodynamik 
beschäftigen, begegnen heute in zunehmendem Maße 
nicht-elementare Integrale und höhere transzendente 
Funktionen, insbesondere elliptische Integrale und 
Funktionen. Ihnen will das vorliegende Buch helfen, 
ohne daß seine Benutzung eine eingehende Beschäfti- 
gung mit diesen Integralen und Funktionen voraus- 
setzt. Es gibt eine Zusammenstellung von etwa 
3000 Integralen und Formeln, die nach praktischen 
Gesichtspunkten ausgewählt sind. Dabei sind in der 
Hauptsache die Formeln und Bezeichnungen nach 
Jacobi benutzt, während die von Weierstraß in 
einem Anhang verwiesen sind. Nach einer kurzen 
Einführung, in der an fünf Beispielen der Gebrauch 
der Tafeln gezeigt wird, werden zunächst die Defi- 
nitionen und Fundamentalrelationen zusammen- 
gestellt. Den größten Teil des Buches füllen dann 
Tafeln für die Reduktion algebraischer, trigono- 
metrischer und hyperbolischer Integranden auf Jaco- 
bische elliptische Funktionen. Es folgen u.a. Ab- 
schnitte mit Tafeln von Integralen Jacobischer 
Funktionen, mit elliptischen Integralen dritter Art, 


grale und Funktionen, insbes 
funktionen, gebracht. 
tafeln der wichtigste 


zeichnete Dienste leisten. 


Dresden. _ Willers. 


J.C. Belgrano unter Mitarbeit von A.Lopez Nieto 
und J.M.Urcelay, Tratado de nomografia 


mit einem Vorwort von J. Rey Pastor. XII + 387 S. 
mit etwa 250 Abb. Madrid 1953. Insitutio tecenico 
de la construction y del cemento. 
125, — ptas. k 

Nach Schilderung der grundlegenden Gedanken, 
die zur nomographischen Darstellung führen, werden 
Kurventafeln mit festen und beweglichen Ebenen er- 
örtert, anschließend die Fluchtlinientafeln mit einer 
Geraden als Ablesemittel und in Erweiterung mit 
einem Winkel, einem Parallel-Lineal oder einem 
Storchschnabel als Ablesemittel. Die breite und 
in viele Einzelheiten gehende Darstellung enthält 
zahlreiche, den verschiedensten Gebieten entnommene 


Beispiele und schließt mit einem Anhang, der neben 


einem Hinweis auf die vielfältigen Möglichkeiten der 
Nomographie einen englischen und einen französischen 
Auszug enthält sowie ein Verzeichnis der wichtigsten 
nomographischen Werke (einzelne Literaturhinweise 
finden sich im Hauptteil an verschiedenen Stellen) 
und eine Zusammenstellung der Anwendungsbeispiele 
geordnet nach Sachgebieten. Das vorliegende Werk 
wird zweifellos der Anwendung der Nomographie sehr 
förderlich sein. 


Aachen. W.Meyer zur Capellen. 


Dr. W.Reichardt, (Prof. a. d. Techn. Hochschule 
Dresden. Grundlagen der Elektro- 
akustik. VIII + 464 S. m. 311 Abb. u. mehreren 


Tabellen. Leipzig 1952. Akademische Verlags- 
gesellschaft Geest u. Portig K. G. Preis geb. 
32,— DM. 


Das folgende Lehrbuch, dessen Kapitel-Über- 
schriften lauten: Physikalische Grundbegriffe, das 
Schallempfinden des Menschen, die Schallwandler, 
Mikrophon-Einsatz, Beschallungstechnik, Elektro- 
Schall-Übertragung, Schall-Aufzeichnung stellt eine 
umfassende Einführung in die technische Akustik 
überhaupt dar. Wer sich darüber hinaus für ein 
Spezialgebiet besonders interessiert, findet alle er- 
forderlichen Hinweise in dem 14 Seiten umfassenden 
Literaturverzeichnis. 

Das spezifisch ‚Elektrische‘ an diesem Werk liegt 
eigentlich nur darin, daß mehr Kenntnisse der 
Elektrotechnik als der Mechanik vorausgesetzt 
werden und daher dem Leser das mechanische Ver- 
halten anhand von Analogien zur Elektrotechnik 
nahegebracht wird. Dabei kann bekanntlich mit 
gleichem Recht der Schalldruck mit der Spannung 
wie mit dem Strom verglichen werden. Der Verfasser 
gibt der zweiten Art den Vorzug, weil sie leichter ge- 
stättet, mechanische Systeme in Schaltungen zu 
übersetzen. Er wehrtsich daher (mit Recht) dagegen, 
den Quotienten aus Druck und Schnelle als ‚„‚Wider- 
stand‘ zu bezeichnen. Er schlägt statt dessen den 
Ausdruck ‚„Standwert‘“ vor und für den reziproken 
„Leitwert‘“ den Ausdruck „Mitgang“. (Doch klingt 
auch hierbei eine größere Zusammengehörigkeit der 
zuletzt genannten Ausdrücke an, die eben doch 
tiefere Gründe hat). 


Preis brosch. 
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Physiker und dem forschenden Ingenieur ausge- Be: 
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ch etaker der ZaMM 
ur wegen seines Gegen- 
‚der Mechanik greift, 
| mehr wegen seiner grundsätzlichen 
 Auseinandersetzungen über die zweckmäßigste Dar- 
stellung von Schwingungs- und Wellenvorgängen. 
Wie sehr es dem Verfasser dabei um Einheitlichkeit 
und Konsequenz geht, zeigen die Zusammenstel- 
lungen über die verwandten mathematischen Schreib- 
weisen, die verwandten Buchstaben und das ver- 
wandte Maß-System, die noch vor der Einleitung 
‚stehen. Die Elektro-Akustik hat es hierin besonders 
schwer, denn in sie gehen mechanische und elek- 
trische Größen ein sowie physikalische Vektoren, die 
gleichzeitig Zeiger der komplexen Ebene bedeuten. 
Als umfassendes Maß-System wird daher das Volt- 
Ampöre-Meter-Sekunde-System nach Giorgi ge- 
wählt, wobei als zusammengesetzte Einheit nur noch 
bei mechanischen Problemen das Newton für die 
Kraft zugelassen wird. Der Versuch des Verfassers, 
jede Verwechslung in der Schreibweise auszu- 
schließen, führt freilich zu einer Zuhilfenahme von 
Neu-Einführungen, die möglicherweise stets die 
Eigenheit dieses Buches bleiben werden. (Sehr an- 
.schaulich erscheint dem Referenten dabei die Kenn- 

zeichnung des reellen Teiles einer komplexen Größe 
durch einen horizontalen und des imaginären Teiles 
durch einen vertikalen Strich.) 

Der Anklang, den das Buch gefunden hat, zeigt 

sich darin, daßin Kürze die zweite Auflage erscheinen 
wird. Bei dieser Überarbeitung hat der Verfasser die 
Mühe nicht gescheut, die Schreibweise entweder den 
inzwischen erfolgten Beschlüssen des AEF anzu- 
passen oder dort, wo solche noch ausstehen, auch 
seine bisherige Schreibweise zu verlassen, und durch 
eine ihm geeigneterer scheinende zu ersetzen, wo- 
möglich zu vereinfachen. 
. Wir wünschen dem Werk noch viele Auflagen mit 
asymptotischer Annäherung an das von ihm erst- 
malig angestrebte Ziel einer allen Bedürfnissen der 
Elektro-Akustik gerecht werdenden Schreibweise, 
die dann allgemein anerkannt und benutzt werden 
möge. 


München. 


L. Cremer. 


Professor Dr. Alfred Lottermoser, (weil. Direktor 
des Kolloidchemischen Instituts an der Technischen 
Hochschule Dresden. Kurze Einführung 
in die Kolloidchemie unter besonderer 


‘ Berücksichtigung der anorganischen Kolloide. Dritte, 


ergänzte Auflage von Dr.-Ing. Carl Kalauch (For- 
schungsinstitut für chemische Technologie Meinsberg 
bei Waldheim). VIII + 260 S.m. 111 Abb. Dresden 
und Leipzig 1954. Verlag Theodor Steinkopff. Preis 
geb. 10,— DM. 

Die bewährte Einführung von A. Lottermoser in 
das Studium der Kolloidchemie ist in der Neubear- 
beitung von C. Kalauch durch eine Reihe von Er- 
gänzungen bereichert worden. Die Darstellung des 
gesamten Stoffes, die sehr elementar und anspruchs- 
los ist, wurde unverändert beibehalten. Neben der 
Ostwald’schen Einteilung der Kolloide werden neuere 
Einteilungsprinzipien, z.B. von Staudinger und 
Kruyt, erwähnt. Bei der Behandlung der Micell- 
bildung wird die Wasserstoffbrückenbindung als 
häufig wesentlicher Faktor eingeführt. Die modernen 
Apparaturen für die Elektrophorese werden ab- 
gebildet, insbesondere auch die Grassmann’schen 
einfachen Vorrichtungen zur Papierelektrophorese. 
Ein kurzer Abschnitt über die Chromatographie ist 
auch eingefügt worden, allerdings ohne die wissen- 
schaftlichen Grundlagen näher zu erörtern. Weiter 
wird das Verfahren der Diasolyse sowie die Messung 
von Dialysepotentialen und ihre Theorie behandelt. 
Die Erscheinungen der Thixotropie werden anhand 
einer Reihe von Abbildungen erläutert. Auch die 


d lebendig ge- neuen 
wie das ronenmikroskop, werden kur 


und technisch wichtige kolloidchemische 
wie der Waschvorgang und der Flotationspr 
werden — wenn auch nur andeutungsweise — 
Grund der neueren Literatur besprochen. BE 
Im Ganzen ist das gerade für den Anfänger sehr 


leicht lesbare Buch von Lottermoser durch die in 


der Neuauflage vorgenommenen Ergänzungen wesent- 


lich bereichert worden und wird unter den Studie- ? 


renden sehr viel Beifall finden, 
Dresden. 


“ Dr, W. Großmann (0, Prof. a.d. T. H. Hannover) 
Grundzüge der Ausgleichsrechnung 
nach .der Methode der kleinsten 
Quadrate nebst Anwendungin der 
Geodäsie. VIII +261 S.m.54 Abb. Berlin/ 
Göttingen/Heidelberg 1953. Springer-Verlag. Preis 
geb. 19,80 DM. 

Wie der Titel sagt, ist das vorliegende Buch ins- 
besondere für Geodäten geschrieben. Es kann als 
Hilfe bei Vorlesungen gebraucht werden, eignet sich 
aber auch als Nachschlagewerk für den Praktiker. 
Entsprechend dem vorgesehenen Benutzerkreis ist die 
Darstellungin der heute meist üblichen Form gehalten, 
Von der Benutzung der Matrizenrechnung wird abge- 
sehen. Daher wird z. B. auch nicht auf die von Bana- 
chiewicz angegebene Umformung der Gauß’schen 
Methode zur Auflösung linearer Gleichungssysteme 
eingegangen, diein der ihr von Zurmühl gegebenen 
Form zur Vereinfachung der Rechenarbeit wesentlich 
beiträgt. Nacheinander werden folgende Abschnitte 
behandelt: Grundzüge der Fehlerlehre, Ausgleich von 
direkten Beobachtungen, Ausgleich von vermittelnden 
Beobachtungen, Ausgleich von bedingten Beob- 
achtungen. Ein Anhang, der sich mit der Ausgleichung 
durch schrittweise Annäherung, dem Glätten von 
Beobachtungsreihen und dem mittleren Fehler von 
Genauigkeitsmassen befaßt, dürfte allgemeineres 
Interesse haben. An zahlreichen Beispielen wird die 
Anwendung der Fehlerrechnung insbesondere auf 
Fragen des Vermessungswesens behandelt. Dabei 
wird Wert auf die praktische Anordnung der Rechnung 
gelegt. Das Buch, das die Mitte zwischen größeren 
Darstellungen, wie es z. B. die von Helmert ist, und 
den vielen heute vorhandenen kleineren hält, kann 
allen Geodäten warm empfohlen werden. 


Dresden.’ Willers. 


Joseph Peres (Prof. & la facalti des sciences de 
Paris, M&canique gönörale. IV+A408S. 
m. 64 Abb. Paris 1953. Masson et Cie. Preis brosch. 
2110 frs., geb. 2545 frs. 

Vorliegendes Buch entspricht im wesentlichen den 
Pariser Vorlesungen des Verfassers und soll vor- 
wiegend dem Gebrauch der Studierenden dienen, Der 
Inhalt umschließt etwa die sonst als analytische 
Mechanik an den Universitäten übliche Vorlesung, 
jedoch mit verschiedenen Erweiterungen. Ebenso 
wie Lagrange in seiner analytischen Mechanik 
hat der Verfasser größtenteils auf Abbildungen ver- 
zichtet, so daß der Stoff in reichlich abstrakter Form 
dargeboten wird. Das Buch gliedert sich in folgende 
Abschnitte: Prinzipien der klassischen Mechanik; 
Mechanik starrer Körper; Methode der virtuellen 
Arbeit und Energiesatz; Untersuchungen einer 
fundamentellen Differentialgleichung der Dynamik; 
Mechanik des materiellen Punktes; klassische Fragen 
zur Mechanik des starren Körpers; die Gleichungen 
von Lagrange und Appell; Variationsprin- 
zipien und kanonische Gleichungen; Bewegungen 
eines Systems bei linearen Differentialgleichungen 
mit konstanten Koeffizienten; Stoßmechanik und 
Ergänzungen zu Systemen mit einsinnig veränder- 
lichen Freiheitsgraden; Mechanik der Kontinua mit 
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H.v.Jürgensonn, Elastizitätund Festig- 
keitim Rohrleitungsbau. Statische 
Berechnung der Rohrleitungen und 
ihrer Einzelteile. Zweite, neubsarbeitete 

_ Auflage. VIII + 379 S. m. 236 Abb., 38 Zahlentafeln 
und 31 Zahlenbeispielen. Berlin-Göttingen-Heidel- 
berg 1953. Springer-Verlag. Preis geb. 39,— DM. 

Die stürmische Entwicklung des Hochdruck-Rohr- 
leitungsbaus hatte den Verfasser zur Herausgabe der 
ersten Auflage des Buches im Jahre 1940 veranlaßt. 
Er kam damit dem Wunsche wsiter Kreise nach einer 
‘Zusammenstellung der Berechnungsgrundlagen für 

 Hochdruck-Rohrleitungen entgegen und brachte 
neben den notwendigsten theoretischen Grundlagen 
ausführliche Anleitungen für die in der Praxis be- 
sonders bewährten Berechnungsverfahren. Diese 
Auflage war seit Kriegsende vergriffen. In der nun- 
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Die besprochenen und angezeigten Bücher sind durch den Buchhandel zu beziehen. 


G.Bangen und Dr.R.Stender, Wahrsch 
lichkeitsrechnung sth 
tische Statistik. (Sch zur 
tung des mathematischen Unter 


) 8.- 2 
mit 28 Abb. und 2 Tafeln. Frankfurt am Main/Pinne- 


berg 1954, Otto Salle-Verlag. Preis 4,40 DM. e 
Die vorliegende Schrift gibt eine für die Hand de 
Lehrers bestimmte Skizze zur einführenden Behan! 
lung der Wahrscheinlichkeitsrechnung und Statist 


im Schulunterricht. Im 1. Teil werden methodische : 


Fragen besprochen. Der 2. gibt den eigentlichen 


Lehrgang. Hier werden zunächst Probleme mit dis- _ 


kreten, sodann mit kontinuierlichen Zufallsveränder- 
lichen, wesentlich in Anlehnung an Gedankengänge 
von v.Mises, behandelt. 
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EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


Prof. Dr.-Ing. €. Weber, Schwingungen im 
Maschinenbau. 151 S. m. 93 Abb. Düsseldorf 
1953. Deutscher Ingenieur-Verlag GmbH. Preis 
kart. 17,— DM (VDI-Mitglieder 15,30 DM.). 


E.Melan (o. Prof. a. d. T. H. Wien) und H. Parkus 
(Prof. am Michigan State College, East-Lansing, Mich. 
USA.), Wärmespannungeninfolge sta- 
tionärer Temperaturfelder 5-+11I4AS. 
m. 30 Abb. Wien 1953. Springer-Verlag. Geb. 

18,50 DM. 


Premier Colloque sur les &quations aux 
deriv6es partielles. (Tenu & Louvain du 17. au 
19. d&cembre 1953). Centre Belge de recherches math6- 
matiques. 1288. Liege, Paris 1954. Masson & Cie, 
G. Thone. Preis geh. 200 fr. belges. 1400 francs fran- 
gais. 


JenaerJahrbuch1953(Wissenschaft- 
liche Veröffentlichungen des Zeiss- 
werkes). Herausgegeben vom VEB Carl Zeiss, 
Jena. 269 S. mit zahlr. Abb. Jena 1953. Kommissions- 
verlag VEB Gustav Fischer. Preis geb. 18,— DM. 


Dr. W. Haack, (o. Prof. an der T. U. Berlin), Dar- 
stellende Geometrie I]. Die wichtigsten 
Darstellungsmethoden. Grund und Aufriß eben- 
flächiger Körper. (Sammlung Göschen Bd. 142) 
110 S.m. 117 Abb. Berlin 1954. Verlag Walter de 
Gruyter & Co. Preis geh. 2,40 DM. 


HigherTranscendentalFunktions, Vo- 
lume I. Based, in part, on notesleft by Harry Bate- 
man (Late Prof. at the California Institute of Tech- 
nology) and compiled by the Staff of the Bateman 
Manuscript Projekt. XXVI + 302 S. New York/Te- 
ronto/London 1953. Verlag MeGraw-Hill Book 
Company, Inc. Preis geb. 56/6s. . 


Dr.K.Bechert (o. Prof. an der Univ. Mainz) und 
Dr.Ch.Gerthsen (o. Prof. an der T.H. Karlsruhe). 
Atomphysik. Bd.Ill. Theorie des Atom- 
baus, I. Teil von Dr. K. Bechert. 3. umgearb. 
Aufl. 148 S. m. 16 Abb. Berlin 1954. Bd. IV. 
Theorie des Atomsbaus, II. Teil von Dr. K. 
Bechert. 3. umgearb. Aufl. 180 S. m. 14 Abb. 
Berlin 1954. (Sammlung Göschen Bd. 1123 bzw. 
1165). Verlag Walter de Gruyter & Co. Preis geh. 
je 4,80 DM. 


Dr.-Ing. W.Matz, Anwendung des Ähn- 
lichkeitsgrundsatzes in der Ver- 
fahrenstechnik. (Verfahrenstechnik in Ein- 
zeldarstellungen Bd. 3) VIII + 115 S.m.29 Abb. 
Berlin/Göttingen/Heidelberg. 1954. Springer-Verlag 
Preis geh. 13,50 DM. > 


Dr.-Ing. habil. W. Kaufmann (o. Prof. an der T.H. 
München). TechnischeHydro- undAero- 
mechanik. VIII-+352S. m.254 Abb. Berlin/ 
Göttingen/Heidelberg 1954. Springer-Verlag. Preis 
geb. 36,— DM. 


Herausgeber und Hauptschriftleiter: Prof. Dr. Dr.h.c. Fr. A. Willers, Dresden. Verlag: Akademie-Verlag GmbH., Berlin W8, Mohrenstraße 39. 


Fernsprecher Sammelnummer: 200386. Postscheckkonto: Berlin 35021. 


Bestell- und Verlagsnummer dieses Heftes: 1009/34/8—9. Die Zeit- 
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„ 


Am Abend des 4. Oktober 1954 entschlief in Landshut (Bayern) 
im 77. Lebensjahr der em. Professor der Technischen Universität 


Berlin-Charlottenburg 


Dr. Dr. E.h. 


Georg Hamel 


Mit ihm ist der letzte der Mathematiker dahin geschieden, die 
der Zeitschrift für angewandte Mathematik und Mechanik seit 
ihrer Gründung nahe standen. Bis zu seinem Tode war er Mit- 
glied des Wissenschaftlichen Beirates der Zeitschrift. Hamel war 
ein außerordentlich vielseitig interessierter Mathematiker, der 
zahlreiche Arbeiten aus den Gebieten der Mathematik wie der 
Mechanik geschrieben hat. Eine Reihe dieser wertvollen Aufsätze 
sind in unserer Zeitschrift veröffentlicht. Schriftleitung Re Ver- 


lag werden seiner Mitarbeit stets dankbar gedenken. 


Willers 
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